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CHAPTER 1 



Introduction 

1.1 Motivation 

The neutrino was proposed by Pauli in 1929 "....as a desperate remedy to save the 
principle of energy conservation in beta decat/^ . This particle which can take part only in 
weak interactions was discovered experimentally by Reines and Cowan in 1956. Neutrinos 
can come in different ffavors in analogy to the different ffavors of charged leptons. However, 
the decay width of Z° boson in the LEP2 experiment restricts the number of light active 
neutrinos to three [Q. We now have experimental evidence for the existence of all the 
three types of neutrinos, z/g, z/^ and Ur- 

The neutrino is known to be a neutral particle which carries a spin 1/2. But whether 
it has mass or not has been an intriguing issue ever since it was proposed. The direct 
upper limits on neutrino masses are quite poor. The limit on z/g mass is obtained from 
tritium beta decay experiments and the best bound is m^^ < 5 eV 0. The bounds on z/^ 
and Ut masses are much more weaker, z/^ < 190 keV and Ur < 18 MeV 0]. If the neutrinos 
are assumed to be Majorana particles then another experimental bound on neutrino mass 
comes from the neutrino-less double beta decay experiments (/3/?0z/ or PPOuJ, where J is 
a Majoron). This is a lepton number violating process which will be possible only if the 
neutrinos are massive Majorana particles. From the non-observation of this process the 
most stringent bounds on the Majorana neutrino mass is m^^ < 0.35 eV at 90% CL. 0]. 

Even though these upper bounds, particularly the ones coming from direct mass 
searches are still weak (these limits on neutrino masses far exceed the cosmological bound 
which we will discuss later), they clearly indicate that the neutrino masses are much 
smaller than the mass of the corresponding charged leptons. This leads to an intra- 
familial hierarchy problem, a challenge for any model which can predict neutrino mass. 
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In the Glashow-Weinberg-Salam standard model of particle physics, which is consistent 
with all known experimental data till date, the neutrinos are assumed to be massless. 
However there is no underlying gauge symmetry which forbids neutrino mass, unlike as 
in the case of the photons. In most extensions of the standard model, the Grand Unified 
theories and the supersymmetric models, the neutrino is massive 0. 

If one does allow for non-zero neutrino mass then the flavor eigenstates of the neutrino 
can be different from the mass eigenstates. One then encounters the quantum mechanical 
phenomenon of neutrino flavor mixing where one neutrino flavor oscillates into another 
flavor due to interference effects. The parameters involved in this process is the mass 
square difference between the two states which mix, Am^ and the mixing angle between 
them, 6. This mechanism can probe very small neutrino masses. Neutrino flavor oscilla- 
tions was conjectured long ago as a plausible solution to the solar neutrino problem |^ 
and the atmospheric neutrino anomaly p. 

The thermonuclear fusion reactions responsible for energy generation in the Sun release 
a huge flux of neutrinos. This flux of pure u^, arriving from the Sun have been measured for 
almost the last 40 years now by the Homestake, SAGE, GALLEX/GNO, Kamiokande, 
Super-Kamiokande and the SNO experiments. All these experiments have observed a 
deficit of the solar neutrino flux predicted by the "standard solar model" (SSM) 
This discrepancy between theory and experiment came to be known as the solar neutrino 
problem. Neutrino flavor mixing - either in vacuum or in solar matter |[TT[ 



can 



account for the solution to this apparent anomaly. The most favored solution to the 
global data demands a Am^ ~ 10~^ eV^ and large values for the mixing angle. 

The atmospheric neutrinos are produced due to collision of the cosmic ray particles 
with the air nuclei. These neutrinos were detected in large water Cerenkov experiments 
which reported a depletion of the z/^ flux compared to expectation. Solution to this 
atmospheric neutrino anomaly again called for neutrino flavor oscillations with Am^ ~ 
10~^ eV^ and maximal mixing to reconcile data with predictions. The results from the 
Super-Kamiokande (SK) atmospheric neutrino experiment in 1998 became a hallmark 
in the history of particle physics when it flnally conflrmed that the atmospheric neutrinos 
do oscillate and are hence indeed massive. The observation of neutrino mass in the SK 
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experiment (although indirect) is the first and till date the only evidence of physics beyond 
the standard model. 

There have been many terrestrial neutrino oscillation searches using both accelera- 



tors as well as reactors as neutrino sources [12, 13]. But all of them except the LSND 



experiment at Los Alamos, have yielded negative results for neutrino oscillations. The 
LSND experiment has continued to give positive signal for oscillations since 1996 with 
Am^ ~ eV^ WM. Thus we have three indications of neutrino oscillations. However since 
the three different hints demand three different values of Am^ , it is difficult to explain all 
the experimental data in the framework of three neutrinos. There have been quite a few 
attempts in the literature to explain all the three experiments with three flavors but it is 
largely believed that if the LSND results are correct then one has to introduce a fourth 
sterile neutrino. 

Neutrinos are known to play a pivotal role in the supernova dynamics and nucleosyn- 
thesis. A huge flux of neutrinos and antineutrinos are released during the thermal cooling 
phase of a core collapse supernova. These neutrinos can be detected in the terrestrial 
detectors. The detection of the neutrinos from SN1987A in the Kamiokande and 1MB 
p!5[] gave birth to the subject of neutrino astronomy and heralded the beginning of a new 
era in neutrino physics. A careful study of the resultant neutrino signal from a galactic 
supernova in the terrestrial detectors, can throw light on not just the type-II supernova 
mechanism but also on the neutrino mass and mixing parameters. 

The neutrinos are the most abundant entities in the Universe after radiation. Hence 
even a small mass for the neutrinos can make a huge difference to the total mass of the 
Universe. The energy density of the neutrinos cannot exceed the total energy density of 
the Universe and this puts a strong upper bound on the sum of all the light neutrino 
species |]16[ , J2i < 46 eV . Small neutrino masses can be an important component 



of the dark matter. Since neutrinos were relativistic at the time of structure formation 
they are called "hot" dark matter. It is believed that a combination of hot+cold dark 
matter is required for a correct explanation of this problem. Neutrino properties and 
number of neutrino generations are also severely constrained by primordial nucleosynthesis 
arguments. 



Neutrinos play a crucial role in the understanding of stellar evolution, type-II super- 
nova mechanism, nucleosynthesis and dark matter studies. A huge amount of theoretical 
and experimental effort has gone into the study of neutrino properties. Thus the im- 
portance of the subject warrants a detailed analysis of the neutrino mass and mixing 
parameters in the context of the current experimental data and a careful study of its 
implication for astrophysics. 

1.2 Plan of Thesis 

In this thesis we explore the implications of neutrino oscillations in the context of solar 
neutrino problem, atmospheric neutrino anomaly and the terrestrial accelerator /reactor 
neutrino oscillation experiments. We perform detailed statistical analyses of the solar and 
atmospheric neutrino data and map out the regions of the parameters space consistent 
with the experiments. We investigate the effect of neutrino mass and mixing on the pre- 
dicted neutrino signal from a galactic supernova in the current water Cerenkov detectors. 

We begin in chapter 2 with the presentation of the basic aspects of neutrino oscillations 
in vacuum and in matter. We discuss the adiabatic and non-adiabatic propagation of 
neutrinos in a medium with varying density and present the expressions for the survival 
probability. 

In chapter 3 we give a detailed description of the current status of the neutrino oscil- 
lation experiments in terms of two flavor oscillations. For the solar neutrinos we present a 
brief review of the SSM predictions, summarize the main experimental results, work out 
the unified formalism for the solar neutrino survival probability, describe our solar neu- 
trino code and perform a comprehensive analysis of the global solar data. Similarly for 
the atmospheric neutrinos we discuss the theoretical flux predictions, the SK experimental 
data, our atmospheric code and perform a fit with the z/^ — oscillation scenario. We 
briefly review the bounds from the most important accelerator/reactor experiments. 

We extend our study of the solar neutrinos in chapter 4 where we probe the potential 
of the energy independent solution in describing the global solar neutrino data. We 
investigate the signatures of this scenario in the future solar neutrino experiments. 
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In chapter 5 we analyse the SK atmospheric neutrino data and the accelerator/reactor 
data in a three- generation framework. We present the results of the analysis of (1) 
only the SK data and (2) SK+CHOOZ data and compare the allowed regions with those 
obtained from the accelerator/reactor experiments. For the only SK analysis we indicate 
some new allowed regions with very small Am^ which appear due to matter effects peculiar 
to the mass spectrum of the neutrinos that we have chosen here. We study the implications 
of this mass spectrum in the K2K experiment. 

In chapter 6 we work in a scheme where one of the components in the neutrino beam 
is unstable and explore the viability of this decay model as a solution to the atmospheric 
neutrino problem. 

In chapter 7 we make quantitative predictions for the number of events recorded in the 
current water Cerenkov detectors due a galactic supernova. We examine the signatures 
of neutrino mass and mixing which can show up in the detectors and suggest various 
variables that can be used to study the effect of oscillations in the resultant signal. 

We present our conclusions in chapter 8 with a few comments on the physics potential 
of the most exciting future detectors. 
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CHAPTER 2 



Neutrino Oscillations 



Neutrino oscillations in vacuum is analogous to — oscillations in its quantum 
mechanical nature. Neutrinos produced in their flavor eigenstates lack definite mass if 
they are massive and mixed. The states with definite mass are the mass eigenstates 
and neutrino propagation from the production to the detection point is governed by the 
equation of motion of the mass eigenstates. Finally the neutrinos are detected by weak 
interactions involving the fiavor states, leading to the phenomenon of neutrino oscillations, 
first suggested by Bruno Pontecorvo [|T| and later by Maki et a/ . Some very good reviews 
on neutrino oscillations include [^, H, ^ ^, |^ 

When neutrinos move through matter they interact with the ambient electrons, pro- 
tons and neutrons. As a result of this they pick up an effective mass much the same way 
as photons acquire mass on moving through a medium. This phenomenon has non-trivial 
impact on the mixing scenario of the neutrinos. In section 2.1 we first develop the for- 
malism for neutrino oscillations in vacuum. In section 2.2 we discuss in detail the effect 
of matter on the mass and mixing parameters of the neutrinos. 

2.1 Neutrino Oscillations in Vacuum 

The flavor eigenstate \va) created in a weak interaction process can be expressed as a 
linear superposition of the mass eigenstates |z/j) 



where U is the unitary mixing matrix analogous to the Cabibo-Kobayashi-Maskawa ma- 
trix in the quark sector. We consider the general case of neutrinos with flavors and 



N 




i=l 
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henceforth set ^ = c = 1. After time t, the initial \va) evolves to 



N 



(t))=^e-^^'U.k.) (2.2) 



1=1 



where Ei is the energy of the i mass eigenstate. For simplicity we assume that the 3- 
momentum p of the different components of the neutrino beam are the same. The energy 
Ei of the i*^^ component is given by 



= v/p2 + mf (2.3) 

where rrii is the mass of the i*^ mass eigenstate. However since the masses are non- 
degenerate, the Ei are different for the different components and eq. (|2.2|) is a different 
superposition of |z/j) compared to eq. (|2.1| ). Hence one expects the presence of other flavor 
states in the resultant beam in addition to the original flavor. The amplitude of finding 
a flavor z/^j in the original Ua beam is 

N 

{MMt)) = J2e-'^''U^.U;i (2.4) 

i=l 

Hence the corresponding probability is given by 

N 

P.^., = Y.\Up^\'\Ua^\' + ^eY.U;,UpjU^.U:^ eH(i^.-i^.)*] (2.5) 

i=l i^j 

For ultra relativistic neutrinos with a common definite momentum p, Ei ^ \p\ + and 
t can be replaced by the distance traveled L. Then we obtain 

Pu^., = - 4 J2U^,U;,U:^Ufs, sin' (^] (2.6) 

j>i \^ijJ 

is the oscillation wavelength which denotes the scale over which neutrino oscillation effects 
can be significant and Ajj = mf — m?-. The actual forms of the various survival and 
transition probabilities depend on the neutrino mass spectrum assumed and the choice of 
the mixing matrix U relating the flavor eigenstates to the mass eigenstates. The oscillatory 
character is embedded in the sin^ (f^) ^^^i^- Depending on the value of E and Ajj, if 
the oscillation wavelength is such that Xij ^ L, sin^ (7rL/Ajj) 0, the oscillations do 
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not get a chance to develop and the neutrino survival probability is ~ 1. On the other 
hand, Xij <^ L implies a large number of oscillations, so that once the averaging over 
energy and/or the distance traveled is done sin^ (vrL/Ajj) ^1/2 and one encounters what 
is called average oscillations. But when Ajj ~ L then one has full oscillation effects and 
constraints on Aij can be put from observations of the resultant neutrino beam. 

If we restrict ourselves to two-generations for simplicity then the mixing matrix takes 
a simple form 

U=( (2.8) 

V — sm u cos a J 

z/g = cos 9^1 + sin 9^2 (2.9) 
z/^ = — sin 9^1 + cos 9^2 (2. 10) 

where 9 is called the mixing angle. The expression of the transition probability from 

to a different flavor z/^ reduces to 



y) 



P,^,^ = sm^29sm^{^] (2.11) 



The survival probability of the original neutrino beam is simply 



P = I — P 



1 - sin^ 2^ sinM — ) (2.12) 



A 



2.2 Neutrino Oscillations in Matter 



So far we have discussed the neutrino transition and survival probabilities in vacuum. 
When neutrinos move through a medium they interact with the ambient matter and this 
interaction modifies their effective masses and mixing. It was pointed out in [P, |lOl that 



the patterns of neutrino oscillations might be significantly affected if the neutrinos travel 
through a material medium. This is because normal matter has only electrons and so the 
z/g experience both charged as well as neutral current interactions while the and z/,- can 
participate in only neutral current processes. 

Since the flavor eigenstates are involved in weak interactions we work in the flavor 
basis for the time being and will introduce neutrino mixing later. All the expression for 
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the interaction terms given below are for neutrinos. The corresponding expressions for 
antineutrinos are same with a -ve sign. The charged current scattering of with electrons 
gives the following contribution to the Lagrangian 

>Ce// = -^{e7'^(l - 75)^e}{^e7M(l " l^)e} (2.13) 

where G^r is the Fermi coupling constant. On rearranging the spinors using Fierz trans- 
formation one gets 

^eff = -^{i7el^{l - l,)ye}{ei,{l - l,)e} (2.14) 

For forward scattering of neutrinos off electrons, the neutrino momentum remains the 
same so that charged current contribution after averaging the electron field over the 

background becomes 

- ^{^e7'^(l - l>e}{ei,{l - 75)e) (2.15) 

The axial current part of (67^(1 — 75)6) gives the spin in the non-relativistic limit while the 
spatial part of the vector component gives the average velocity, both of which are negligible 
for a non-relativistic collection of electrons. So the only non-vanishing contribution comes 
from the 70 component which gives the electron density. 

(e7oe) = ne (2-16) 

where rig is the ambient electron density. The extra contribution to the Lagrangian 
then reduces to —\/2GF'nf.i'eLl^^eL- Hence for unpolarized electrons at rest, the forward 
charged current scattering of neutrinos off electron gives rise to an effective potential 

Fee = \/2Gi.ne (2.17) 

The effect of this term is to change the effective energy of the i^e in matter to 



Eeff = Jp^ + m'' + V2GFne 



2 

m 



IpI + 7^ + V2GFne (2.18) 
2|P| 

\p\ + :^{m^ + 2V2GFneE) 
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Thus the charged current interaction gives rise to an extra effective contribution to the 
mass square called A or the Wolfenstein term 

A = 2V2GFneE (2.19) 

The neutral current contribution to the effective neutrino energies can be calculated in 
an identical manner and comes out to be 

Vnc = V2Gf Y1 [4l - 2 sin^ 9wQ^ ] (2.20) 
/ 

where / stands for the electron, proton or neutron, is the density of / in the surrounding 
matter, is the charge of / and I^^^ is the third component of weak isospin of the left 
chiral projection of /. Since for the proton = 1 and I^j^ = 1/2, while for the electron 
Qf = -1 and //^ = -1/2 and since normal matter is charge neutral ensuring that He = n. 



pi 

the contributions due to electron and proton cancel each other exactly. Hence the only 
remaining contribution is due to the neutrons for which Qf = Q and j/j^ = —1/2 so that 

Vnc = -V2GFnn/2 (2.21) 

In order to see the effect of these extra contributions due to interactions of the neutrino 
beam with the ambient matter, on neutrino mass and mixing parameters, let us look at 
the equation of motion for the neutrino states. We restrict ourselves to a two-generation 
scheme for simplicity {u^ mixing with either z/^ or z/^) and first write down the evolution 
equation for the mass eigenstates in vacuum which can be subsequently extended to 
include matter effects. The equation of motion for the mass eigenstates in vacuum can 
be written as 

'i(:!)=Mi:!) (2.22) 



dt \l'2 J \l^2 

where is the mass matrix, diagonal in the mass basis 

As the interaction terms are defined in the flavor basis we have to look at the evolution 
equation in the flavor basis. Since the mass eigenstates are related to the flavor eigenstates 
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by the relation and since the mixing matrix in two-generations is given by (|2.8|) , the 
equation of motion in terms of the flavor states can be written as 

M] = U M (2.25) 

m[+m| A /-cos2^^ sin2^\ 
~ ^ 4E ^ 4E\ sin 2^ cos 2^ J ^^'^^^ 

where Mj is the vacuum mass matrix in the flavor basis and A = m| — m^. We next 
include in the mass matrix the interaction terms in the mass matrix for neutrinos moving 
through matter. Once these extra contributions to the Hamiltonian due to interactions of 
the neutrino with matter are taken into account, the equation of motion for the neutrino 
beam in the flavor basis is given by eq. (|]24D with the Mj in vacuum replaced with Mj^ 



in matter given by 

cos 26 A sin 26 
^Je\ a sin 2^ Acos2e 

If A- Acos26 Asin2^ \ , . 

^4E\ A sin 2^ -A + A cos 2^7 ^ 

where A = 2\p2G Fnf.E as defined before. The terms proportional to the identity matrix 
play absolutely no role in flavor mixing and can be safely dropped. Hence we see that the 
neutral current term which affect all the flavors equally, falls out of the oscillation analysis 
and the only non-trivial contribution comes from the charged current interaction. 

The energy eigenvalues in matter are obtained by diagonalising the mass matrix Mj^. 
We define E"^ to be the mass eigenvalues and Um fo be the mixing matrix in matter. In 
analogy with the vacuum case Vm can be parametrized as 

/COB 9 Bine X 

V-smftm CO?, 6m J 
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where 9m is the mixing angle in matter. Then since 



m fm 

the eigenvalues of the mass matrix are 



7-rt jj _( E"^ 



EJ: 



(2.30) 



1 



(2.31) 



where 



mi + m^2 + A)T V ("^ + ^ cos 26)^ + (A sin 2ey 



(2.32) 



Hence the mass squared difference A in vacuum is modified in presence of matter to 

1/2 



-A + A cos 29y + (A sin 29)' 



2.33) 



while the mixing angle in matter is given by 



tan 29r^ 



A sin 2^ 



or equivalently 



sin^ 29m 



-A + A cos 29 



(A sin 2^)^ 



(2.34) 



(2.35) 



{-A + A cos 29y + (A sin 29^ 

The mass and the mixing angle in matter are changed substantially compared to their 
vacuum values depending on the value of A relative to A. This change is most dramatic 
when the condition 



A = A cos 261 



(2.36) 



is satisfied. If this condition is attained in matter then sin^ 29m = 1 and the mixing angle 
in matter become maximal irrespective of the value of the mixing angle in vacuum. So even 
a small mixing in vacuum can be amplified to maximal mixing due to matter effects This 
is called the matter enhanced resonance effect or the the Mikhevey-Smirnov-Wolfenstein 
or the MSW effect TU, n| which is the favored solution to the Solar Neutrino Problem, 
a subject we shall address later in great details. Note that for A > one encounters this 
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resonance for neutrinos {A > 0) only ii 6 < 7i/4. Also note that the expression ( p.34| ) can 
be rewritten as 

tan 2^ 

tan2^^^ = - ^ 2.37 

1 — Tie /rfj^^ 

where is the electron density at the point of neutrino production and Tf^^'^ is the 
corresponding density at resonance. So that for vFJ ^ rfj'^ the mixing angle in matter 
Om — ^ 7r/2. 

We next look for the electron neutrino survival probability. For a constant density 
medium the survival probability is given by 

P,^,^ = l-sm^2e^sm^{^] (2.38) 



m ' 

2 



A. 



Hence for a constant density medium the expressions for transition and survival proba- 
bilities are exactly similar to ones for the vacuum case, with the mass squared difference 
and the mixing angle replaced by the corresponding terms in matter. But things get more 
complicated when the neutrinos move through a medium of varying density. We address 
this complex issue in the next section. 

2.2.1 Neutrino Survival Probability in Medium with Varying Density 

In most situations encountered in nature the neutrinos propagate through medium 
with varying density. Since the eigenstates and the eigenvalues of the mass matrix both 
depend on the density of the medium, the mass eigenstates defined in eq. ( |2.28| ) are 



no longer the stationary eigenstates. For varying density one can define the stationary 
eigenstates only for the Hamiltonian at a given point. Starting from the equation of 
motion in the flavor basis and dropping the irrelevant terms proportional to the identity 
matrix we have 



^ / A - A cos 2^ A sin 2Q 
"EV Asin26' 

l^^A n .^2 Li (2.40) 



d / Uf. 

d^Kiy^J ^ AE\ Asin2e -A + Acos2e J Ku,, 
J_ (Ml 
2E "I 
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where we have substituted x for t. Keeping in mind that the mixing matrix is now x 
dependent we obtain 



d 









'Ml 








2E^ 


V 





uU—u 



MI/2E -ide^/dx\( 



.m 
,ra 



(2.41) 



idQJdx Mil IE 

The off diagonal terms in eq. ( p.41| ) mix the states z/|" and z/™ and the mass eigenstates in 
matter keep changing with x. Hence one has to solve the eq. ( |2.41| ) to get the survival and 
transition probabilities. However the mass eigenstates travel approximately unchanged 
and unmixed as long as the off diagonal terms in ( |2.41| ) are small compared to the diagonal 
terms. Since the terms proportional to the unit matrix do not affect oscillation probabil- 
ities, the only physical parameter in the diagonal elements is their difference. Thus the 
condition for the off diagonal term to be small can be written as 



dQ^ 



dx 



M| - Ml 



4E 



(2.42) 



This is called the adiabatic condition. Since \d6m/dx\ has a maximum at resonance 
while |M| - Mil has a minimum, the adiabaticity condition becomes most stringent at 
the resonance. One can define an adiabaticity parameter as 



7 



Asin2 20 
2Ecos 26 



d 

dx 



Inrif 



(2.43) 



The adiabatic condition then reduces to 7 ^ 1. Depending on whether the adiabatic 
condition is satisfied or not neutrino propagation in matter may be of two types: adiabatic 
and non- adiabatic. 

•The Adiabatic Case: If one can neglect the off diagonal terms in eq. ( p.41| ) in 

comparison to the diagonal elements, i/™ and z/™ approximately become the eigenstate of 
the mass matrix. This is the adiabatic approximation. As long as the adiabatic condition 
is satisfied, the two mass states evolve independently and there is no mixing between 
them. The state created in matter then remains a z/[" always. This is the adiabatic 
propagation of neutrinos in matter. Adiabatic propagation of neutrinos in matter has very 
interesting consequences. For example, consider the neutrinos produced at the center of 
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the Sun or a supernova. Since the matter density at the point of production of these 
neutrinos is much higher that the resonance density, the mixing angle 9m ~ ti'/S and since 

Ve = cos 6mi^'^ + siu (2.44) 

the i^e are created almost entirely in the heavier z/^ state. Thereafter the neutrino moves 
adiabatically towards lower densities, crosses the resonance when A = A cos 26 and finally 
comes out into the vacuum. Since the motion was adiabatic the heavier state remains the 
heavier state even after it comes out. Since in vacuum the heavier state is 

1^2 = — sin + cos OUfj, (2-45) 

the survival probability of the electron neutrino is just sin^ 9. Hence for very small values 
of the mixing angle 6 one may have almost complete conversion of the z^e produced inside 
the sun. 

In the discussion above the survival probability is sin^ 9 for 9^ = vr/2. For any general 
angle 9m at the production point of the neutrino, the survival probability is 

P.... = ^ + ^ cos 2^ cos 29m (2.46) 

•The Non-Adiabatic case: One may encounter cases where 7 ~ 1 and the adiabatic 
condition ( |2.42D is not satisfied. This breakdown of adiabaticity is most pronounced at the 



position of resonance as discussed in the previous section. This violation of adiabaticity 
at the resonance signals that the off diagonal terms in eq. ( p. 41 ) become comparable to 



the diagonal terms and there is a finite probability of transition from one mass eigenstate 
to another. This transition probability between the mass eigenstates is called the level 
crossing or the jump probability. It is defined as 

Pj = \{u^ix^)\uT{x.))f (2.47) 

where x± refer to two faraway points on either side of the resonance. Pj can be determined 
by solving the equation of motion ( p.41| ) for a given matter density profile. For a linearly 



varying density, the jump probability is given by the Landau-Zener expression [12, 13 



Pj = exp(-^7) (2.48) 
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where 7 is the adiabaticity parameter given by eq. ( 2.43 ). For the Sun the density profile 
is roughly exponential and we will address that issue in the next chapter. 

The electron neutrino survival probability, taking into account the finite level crossing 
between the mass eigenstates due to breakdown of adiabaticity is given by 

P..u^ = \ + {\-Pj) COS 26 cos 26^ (2.49) 

This can be easily computed if one knows the form of Pj. 

2.2.2 Matter Effects with Sterile Neutrinos 

In the entire discussion above on the effect of matter on neutrino mass and mixing 
we had tacitly assumed that both the fiavors involved in oscillations were active fiavors, 
that is, can take part in weak interactions. But one may consider situations where one 
of the states involved is sterile. Sterile states have no interactions with the surrounding 
medium, neither charged nor neutral. Thus for the z/g — VsterWe mixing, relevant for the 
solar neutrino problem, the mixing angle and the mass squared difference in matter are 
given by, 

= \{-2V2GFneE + V2GFnnE + ^cos2ef + {^siYi2efY'^ (2.50) 
A sin 2^ 

tan2^„ = j= j= 2.51 

-2y2G^neE + y2GFn„E + Acos2^ ^ ' 

The survival probability Pu^u^ is still given by eqs. (|2.46|) and ( p.49|) . For the z/^ — Vstenie 
mixing, relevant for the atmospheric neutrino anomaly the corresponding expressions are 

A„ = [{V2G FUnE + ^cos2ef + {^siYi2efY^^ (2.52) 

tan 20^ = ^2.53) 

^GfUuE + a cos 2^ 

Note however that whether the neutrinos are oscillating into active or sterile species 
matter only when a material medium is involved. The oscillations in vacuum is the same 
for both the cases. 



17 



BIBLIOGRAPHY 

[1] B. Pontecorvo, Sov. Phys. JETP 6, 429 (1958). 

[2] Z. Maki et ai, Prog. Theor. Phys. 28, 870 (1962). 

[3] S. M. Bilenky, B. Pontecorvo, Phys. Rep. 41, 225 (1978). 

[4] S. M. Bilenky, S. T. Petcov, Rev. Mod. Phys. 59, 671 (1987). 

[5] T. K. Kuo, J. Pantaleone, Rev. Mod. Phys. 61, 937 (1989). 

[6] Neutrino Astrophysics, by J.N. Bahcall, Cambridge University Press, 1989. 

[7] Massive Neutrinos in Physics and Astrophysics, by R.N. Mohapatra and P.B. Pal, 
World Scientific, 1998. 

[8] Neutrinos in Physics and Astrophysics, by C.W. Kim and A. Pevsner, Harwood 
Academic Pubhshers, 1993. 

[9] L. Wolfenstein, Phys. Rev. D17, 2369 (1978); Phys. Rev. D20, 2634 (1979). 

[10] S. P. Mikheyev and A. Yu. Smirnov, Sov. J. Nucl. Phys. 42, 913 (1985); Nuovo 
Cimento, 9C, 17 (1986). 

[11] H. A. Bethe, Phys. Rev. Lett. 56, 1305 (1986). 

[12] W. C. Haxton, Phys. Rev. Lett. 57, 1271 (1986). 

[13] S. Parke, Phys. Rev. Lett. 57, 1275 (1986). 



18 



CHAPTER 3 



Bounds on Two Flavor Neutrino Mixing Parameters 



The question whether neutrinos are massive or not has been answered. The Super- 
Kamiokande (SK) atmospheric neutrino data has quelled all apprehensions about the 
existence of neutrino flavor mixing and has propelled neutrino physics to the forefront of 
particle phenomenology. The other puzzle which has warranted neutrino oscillations as a 
plausible solution is the long standing solar neutrino problem. For the last four decades 
solar neutrino detectors have recorded a flux far less than that predicted by solar models. 
Though the theory of neutrino flavor oscillations can offer the best possible solution to 
this anomaly, there are still a lot of issues that have to be settled. Finally there have been 
earnest searches for neutrino flavor mixing in the laboratory, using both reactors as well 
as accelerators as sources for neutrino beams. But all such quests have lead to negative 
results, with the exception of the LSND experiment at Los Alamos which has reiterated 
since 1996 to have observed positive neutrino oscillation signals. 

In this chapter we survey the current status of the solar neutrino problem, the atmo- 
spheric neutrino anomaly and the laboratory experiments. We briefly discuss the incident 
neutrino fluxes, the detection techniques, present the main experimental results and per- 
form detailed statistical analysis of the solar and the atmospheric neutrino data in terms 
of two-generation neutrino oscillations. We describe our solar and atmospheric neutrino 
code and discuss the method of analysis. Wc identify the bcst-flt solutions and display 
the allowed regions in the neutrino parameter space. We begin with the solar neutrino 
problem in section 3.1, next take up the case of the atmospheric neutrino anomaly in sec- 
tion 3.2 and finally move over to the description of the most stringent accelerator /reactor 
experiments in section 3.3. 
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3.1 The Solar Neutrino Problem 



Neutrinos are an essential byproduct of the thermonuclear energy generation process 
inside the Sun whereby four proton nuclei are fused into an alpha particle. 

4p — ^ ^i/e + 2e+ + 2z/e + 28 MeV (3.1) 

This process is called Hydrogen Burning and is responsible for the hydrostatic equilibrium 
of the Sun. About 2-3% of the total energy released in the process is carried away by 
the neutrinos. The rest is in the form of electromagnetic radiation, which diffuses out from 
the core to the surface of the Sun, getting degraded in frequency to appear as sunlight. 
It takes millions of years for the photons to emerge from the Sun due to interactions with 
the solar matter. As a result of these repeated scatterings, the photons cannot give us 
any direct information about the core of the Sun. The neutrinos on the other hand have 
typical scattering cross sections of about 10~^^ cm^, which for a solar density of about 
100 g/cc gives a mean free path of the order of 10^'' cm. This is much larger than the 
radius of the Sun. The neutrinos escape from the Sun unadulterated and bring along, 
all the information about the solar core imprinted on them. Thus a careful study of 
these neutrinos promises to provide detailed information about the solar interiors and the 
thermonuclear energy generation process inside the core and hence holds the potential to 
verify or refute the existing solar models. 

Keeping this in mind, Ray Davis began his pioneering experiment in which neutrinos 
are captured by the ^'^Cl atoms The corresponding work on theoretical predictions 
for the solar neutrino fluxes were done by Bahcall and his collaborators in the frame- 
work of what is called the "standard solar model". The first results of this experiment 
were declared in 1968 [|ll. This experiment has been running for almost four decades 
now Pi and the results have been remarkably robust; the experiment sees a depletion 
of the solar neutrino flux over the standard solar model (SSM) predictions. The solar 
neutrinos seem to vanish on their way from the Sun to the Earth and this mystery of 
missing solar neutrinos constitutes the Solar Neutrino Problem (SNP). Two decades later 
the Kamiokande water Cerenkov experiment corroborated this observed deficit of solar 
neutrinos 0]. This conflict between the SSM prediction and observation has been further 
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substantiated by the results from the '^^Ga experiments - the SAGE 0, the GALLEX 
and now the GNO [0, which is the upgraded version of the GALLEX. The advent 
of the Super-Kamiokande |T0|, the upgraded version the the original Kamiokande 

experiment, brought with it rich statistics, which provided further insight into the SNP in 
terms of the overall depletion of the solar flux |^ , the energy dependence of the suppres- 
sion rate [^] and the presence of any difference in the observed rate at day and at night, a 
phenomenon called the day-night effect |TD|. The recently declared results on the electron 
scattering events (ES) from the Sudbury Neutrino Observatory (SNO) |Tl| are consistent 
with the observations of SK while their charged current (CC) events when compared with 
the SK events signify the presence of a u^/ Uj. component in the resultant neutrino beam 
at the 3.3(7 level while the total solar flux is calculated to be 5.44 ± 0.99 x lO^cm^^s^^ 
which is in agreement with the SSM predictions 



Thus if the neutrinos are assumed to be the standard particle predicted by the Glashow- 
Weinberg-Salam standard model of particle physics then there seems to be an apparent 
discrepancy between the experimental observations and the SSM predictions. There have 
been several attempts in the literature to explain the experimental results by modifying 
the solar models, assuming neutrinos to be standard. These endevours are collectively 
called the astrophysical solutions. In the most general class of these solutions the solar 
neutrino fluxes are considered as free parameters, with the only requirement being the 



reproduction of the observed solar luminosity |]12| . But all such analyses fail to explain the 
observations from all the three experiments, CI, Ga and water Cerenkov simultaneously. 
In fact the best-fit for these solutions predict a negative flux for the ^Be neutrino which 
is an unphysical situation. The astrophysical solutions fail to explain even two of the 
solar neutrino results simultaneously. After the declaration of the SNO results the astro- 



physical solution has fallen into further disfavor [|T3|, |l4[. In ||T3[ Bahcall has shown that 
while prior to SNO the astrophysical solution failed to fit the data at the 2.5a level, after 
including the SNO data it is ruled out at 4a. Thus it is not the uncertainties in the solar 
models that are responsible for this anomaly. In fact the solar models have been remark- 
ably refined in the last four decades and the flux uncertainties have been reduced to a 
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large extent. Hence one needs to consider some non-standard property for the neutrino 
in order to be able to reconcile the data with the standard solar model predictions. 

Among the various particle physics solutions proposed till date, neutrino oscillations 
15] and/or in matter [16, l7|, have been the first and the most appealing 



m vacuum 



solution |T^, [T9|, ^ Neutrino flavor mixing has the potential to explain not only 
the total suppression of the solar flux, but also the energy dependence of the suppression 
factor and the day-night asymmetry. The analysis of the solar neutrino problem in terms 
of neutrino mass and mixing gives four pockets of allowed area in the neutrino oscillation 
parameter space. One of them has Am^ ~ 10~^ eV^ and mixing angles very small and 
is called the small mixing angle (SMA) region. Another has Am^ ~ 10^^ eV^ with large 
mixing angles and is referred to as the large mixing angle (LMA) solution. A third allowed 
zone has Am^ in the range 10^'' — 10^^ eV^ and mixing angle close to maximal. This is 
the LOW-QVO region, LOW stands for low Am^ and QVO for quasi vacuum oscillations. 
The last region is the one associated with vacuum neutrino oscillations with Am^ ~ 10"^'' 
eV^ and mixing close to maximal. While doing an analysis of the solar data one has to 
consider oscillations to either an active flavor {u^, u^) or some sterile species (z/gteriie)- 
The latter case will be different, firstly because the sterile neutrinos will not show up in the 
detectors in the neutral current interactions even though z/^, Vj. can. This feature affects 
both the MSW as well as the vacuum solutions. For the MSW solutions there is an extra 
difference as the sterile neutrinos do not have any interaction with the ambient matter, 
both in the Sun and in the Earth. Thus the effect of matter on the mixing parameters 
will be different for the fstenie as compared to inactive, as discussed in the previous chapter. 

We perform detailed -fits for both the t'e — ^'active and t'e — t'steriie transformations. We 
perform a global analysis of the solar neutrino data from all experiments in the frame- 
work of neutrino mass and mixing. We adopt an unified approach for the presentation 
of the MSW, the vacuum and quasi-vacuum oscillations solutions. We find the best-fit 
solution to the global data and show the allowed regions in the Am? -ia.r? 9 parameter 
space. From the global analysis of all available solar neutrino data, one finds that the 
large mixing angle Pe — inactive MSW solution gives the best fit. 
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Section 3.1.1 is a brief discussion on the solar neutrino flux predictions in the standard 
solar models and the various uncertainties involved. In section 3.1.2 we briefly present 
the essential features of the solar neutrino experiments and summarize their main results. 
In section 3.1.3 we discuss the solar neutrino code developed by us. In section 3.1.4 we 
introduce neutrino flavor mixing and develop the unified formalism for the analysis of the 
SNP in the context of neutrino mixing. Finally in section 3.1.5 we present the results of 
a comprehensive analysis of the solar neutrino data. We identify the best fit solutions 
and give CL. allowed areas in the neutrino parameter space. 

3.1.1 Neutrinos in Standard Solar Models 

"Standard Solar Model" is a solar model constructed with the best available physics 
and input data. Almost all solar models postulate the thermonuclear fusion of protons 
( ^.11 ) to be the main energy generation process in the Sun ^ ^ . The 

eq. (|3.1| ) is actually the compactified form for a chain of reactions in which four hydrogen 
nuclei are fused to form a helium nucleus. This chain of reactions, called the pp chain 



is shown in fig. p.l| . The other series of nuclear reactions in the core of the Sun that 

release neutrinos is the CNO cycle. However, since the CNO cycle becomes important 

only above core temperature ~ 10'' K, it produces only about 1.5% of the total solar 

neutrinos released. Nevertheless the CNO cycle is responsible for three sources of solar 

neutrinos and we call them the ^^A^, ^^O and neutrinos. So one has eight different 

types of solar neutrinos, five produced in the pp chain and three in the CNO cycle. 

From the observed solar luminosity and from the fact that 28 MeV energy are released 

per two electron neutrinos produced in eq. ( ^.1| ), one can make an order of magnitude 

estimate of the solar neutrino flux 

^ luminosity 
'^'^ AnD"^ X I (binding energy of ^He) 
4 X lO^^ergs/s 
47r X (1.5 X 10i3cm)2 x 14 MeV 
^ 10^°cm"2s~^ (3.2) 

The exact solar neutrino flux calculations depend on a number of factors such as the 
nuclear reactions rates, the metallicity (Z/X), the age of the Sun and the opacity. Detailed 

23 





e+ + l>e_ (pp) 


p + e + p ^ 


99.7% 






'^H + p - 


¥ ^He + 7 


84.9% 




15% 


^He -5. 


a + 2p 


^He + p — >■ a - 



0.23% 



^He + a -> ''Be + 7 



15% 

'Be + e- ^ 'Li + Ve_ C Be) 



0.02% 
'Be + p B + 7 



'Li + p — > 2q 
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The ^^-chain 



Figure 3.1: The reactions of the pp-chain. The probabihty of a particular reaction is shown 
as a percentage. The neutrinos are shown underhned. Those with double underlines are 
monoenergetic. We indicate in parenthesis the popular names by which these various 
neutrinos are addressed. 



solar neutrino flux calculations from the standard solar models are available p3, 24, 25, 



251 P^, E3, p3 and agree with the rough estimate made in eq. (^.21) in order of magnitude. 



The first two columns of Table ^]T| list the total solar neutrino fluxes along with their 
±1(7 uncertainties, from the eight different reactions, as given by the year 2000 model of 
Bahcall, Pinsonneault and Basu, which we shall henceforth refer to as BPBOO |2^. In flg. 



we show the energy spectrum of the neutrinos emitted in various reactions of the pp 
chain in BPBOO. Also shown are the ±lcr uncertainties in the various fluxesQ. 
^This figure has been taken from John Bahcall's homepage; www.sns.ias.edu/^jnb/ 

24 




Neutrino Energy (MeV) 

Figure 3.2: Solar neutrino spectrum in the standard solar model as a function of neutrino 
energy. The continuous spectra are in units of cm" -2MeV-^s-\ while the monoenergetic 
lines are in cm~^s~^. Figure shows the ±lo" uncertainties in the model predictions of the 
various fluxes. Also shown are the energy ranges over which the solar neutrino experiments 
are sensitive. 



Among the predictions for the neutrinos from different reactions in the Sun, the flux of 
the neutrino is the most uncertain. In fact the largest contribution to the differences 
in the predictions of the various standard solar models is their choice of different values 
for the Si7 which is the astrophysical S'-factor for the reaction (^Be{p, i')^B). While the 
SSM of Dar and Shaviv uses a value of Sn as low as 17 ± 2 eV barn, the earlier 1992 



and 1995 models of Bahcall and Pinsonneault p^ , |26[| had used 5*17 = 24±2 eV barn. The 
value of Sn used in BPBOO is 19lt eV barn which is the value accepted by the Institute 
of Nuclear Theory (INT) g. 

The solar fluxes are sensitive not just to the uncertainties in the nuclear reaction rates, 
but also to the value of the core temperature Tc. In fact the ^B flux $8^ oc T^^, the pp 
flux $pp oc T~^'^ while the ^Be flux ^r^e oc T^. Hence even a slight increase in the value 
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source 


Flux 
(IQio cm-2s-i) 


CI 

(SNU) 


Ga 
(SNU) 


pp 
pep 
hep 
''Be 

13iV 

15Q 


5.95(1.001°:°!) 

i.40xio-2(i.ooi[i:|ili) 

9.3x10"^ 
4.77xlO-ni-OOlo:lo) 
5 05xl0-^fl 00+S?2l 
5.48x10^2(1. oolg;?^) 
4.80X 10-2(1. 00l[J:?^) 
5.63xlO-4(1.00l[J:i) 


0.0 
0.22 
0.04 
1.15 
5 76 
0.09 
0.33 
0.0 


69.7 
2.8 
0.1 
34.2 
12.1 
3.4 
5.5 
0.1 


Total 






128t? 



Table 3.1: The BPBOO predictions for the solar neutrinos fluxes and neutrino capture 
rates in the CI and Ga detectors. The expected ^B flux in SK is 5.05xl0~^cm~2s~i. 

of Tc can seriously affect all the solar fluxes. In particular, it will sharply raise the ^B and 

^Be fluxes and lower the pp fluxQ. The value of Tc is sensitive to a number of factors 

including the value of the opacity of the solar core. If the value of the opacity is raised, 

it slows down the heat transport, leading to higher core temperatures. The value of the 

opacity in turn depends on the abundance of heavy elements in the Sun or the metallicity. 

Another very important ingredient in the solar models is the inclusion of element diffusion. 

Apart from convection, the two other mechanisms important for transporting solar matter 

are; (1) gravitational settling, which pulls heavier elements towards the center and (2) 

temperature gradient diffusion, which results in pushing lighter elements outward. Both 

of these cause the inward diffusion of ^ife and outward diffusion of ^H. Hence diffusion 

increases the opacity, which results in a higher Tc leading to an increase of ^B and ^Be 

fluxes and decrease of the pp flux. 

^This implies that just by adjusting the value of Tc alone one cannot solve the SNP since CI and SK 
experiment will require a lowering of Tc which will raise the pp flux making the Ga results look all the 
more puzzling. 
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In spite of the various uncertainties involved, only a fraction of which have been 



discussed above, it was shown in ||2^, that if the same input physics is used, then all 
the standard solar models agree with one another to an accuracy of better than 10%. 
Hence though we have presented the results on the total fluxes and the neutrino spectra 
from the standard solar model of Bahcall, Pinsonneault and Basu, predictions by almost 
all the standard models published so far are in reasonable agreement with each other. For 
our analysis of the SNP in terms of neutrino mass and mixing, we have used the latest 
SSM predictions by Bahcall, Pinsonneault and Basu (BPBOO) p9| . 

3.1.2 The Solar Neutrino Experiments 
The CI Experiment (Homestake) 

This is the first and the longest running experiment on solar neutrinos started in 
the sixties by Davis and his collaborators with 615 tons of C2CI4 (perchloroethylene) in 
the Homestake Gold mine in South Dakota |l|]. The neutrino detection process in this 
experiment is 

z/e + ^^C/ -^^'^ Ar + e- (3.3) 

The ^"^Ar atoms are extracted from the detectors at the end of a certain period of time 
and counted by detecting the Auger electron released when the ^^Ar decays by capturing 
a K-shell or an L-shell electron. The reaction ( |3.3| ) has a threshold of 0.814 MeV so that 
the CI experiment predominantly detects the ^Be and neutrinos. It misses out on the 
most abundant and least uncertain pp neutrinos which have a maximum energy of only 
0.42 MeV (cf. fig U). The third column of Table ^ shows the BPBOO predictions H] 
for the neutrino capture rates in the CI experiment for the different neutrino sources. Also 
shown for the CI detector are the total predicted rate and the ±lcr uncertainties in the 
model calculations. The numbers quoted are in a convenient unit called SNU, defined as, 
1 SNU = lO^^^events/target atom/second. The observed rate of solar neutrinos in the 
experiment is 

Observed Rateci = 2.56 ± 0.23 SNU (3.4) 
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Compared to the BPBOO prediction of 7.6tij as in Table [Ol , this gives a ratio of observed 
to expected SSM rate of 0.335 ± 0.029. 

The Ga Experiments (SAGE, GALLEX, GNO) 

These are also radiochemical experiments that use "^^Ga as their detector material. 
The ^^Ga captures a Ve to produce ^^Ge by the reaction 

Ue + '^^Ga Ge + e" (3.5) 

This reaction has a threshold of only 0.233 MeV. Hence the advantage that this detector 
has is that it is capable of seeing the pp neutrinos which are responsible for 98.5% of the 
energy generation of the Sun. Hence the fact that the Ga detector could detect these 
neutrinos confirms the basic postulate of all the solar models, that the Sun generates its 
energy through thermonuclear burning. This itself was a very significant achievement of 
the Ga detectors. 

The SAGE (Soviet American Gallium Experiment) in Russia and GALLEX (Gallium 
Experiment) in Italy are experiments that use this detection technique. The SAGE in 
Baksan Neutrino Observatory uses 60 tons of metallic Ga as the target. The ^^Ge pro- 
duced is separated and counted. The observed rate is 

Observed RaXesAGE = 75.4 ±5;° (stat.) (syst.) SNU (3.6) 

The GALLEX is located in the Gran Sasso laboratory in Italy and uses 30 tons of Ga in 
the form GaCl^ — HCl solution. The observed neutrino rate in GALLEX is 

Observed RaXec allex = 77.5 ±11 SNU (3.7) 

The GALLEX has now finished its run and has been upgraded to the GNO (Gallium 
Neutrino Observatory) which has already given results [0 

Observed Ratecjvo = 65.8 ±l^ol SNU (3.8) 

The combined SAGE and GALLEX+GNO results is 

Observed Rateca = 74.7 ± 5.0 SNU (3.9) 

which is more than 6a away from the SSM predicted rate of 128^7 SNU (cf. Table 
P). 
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The Water Cerenkov Experiments (Kamiokande and Super-Kamiokande) 

The water Cerenkov detectors detect solar neutrinos by tlie forward scattering of elec- 
trons 

i/g + — > z^e + e~ (3.10) 

As it moves, the scattered electron emits Cerenkov light, which is viewed by the huge 
number of photomultiplier tubes covering the entire detector volume. The water detector 
in general has a higher threshold so that it is sensitive to just the and the vanishingly 
small hep neutrinos. But it has many other advantages. It is a real time experiment 
which has directional information. The reaction (p.lO|) is forward peaked and the detector 



can reconstruct the direction of the incoming neutrino from the angle of the emitted 
Cerenkov cone. Earlier the Kamiokande |0] and now the Super-Kamiokande |^ have found 
an excess of events peaking broadly in the solar direction and have thus confirmed that 
the observed neutrinos are indeed coming from the Sun. The water detector can observe 
not just the Uf. as in the radiochemical experiments, but neutrinos and antineutrinos of all 
flavors. Thus it can detect z/^ and Ur through neutral current electron scattering though 
the neutral current scattering cross-section is about 1/6*^^ of the charged current scattering 
cross-section. This is important if one wants to distinguish between neutrinos oscillating 
out into either u^/ur or to some sterile species t'steriie, which does not have any standard 
model interactions. But the real strength of this experiment lies in its ability to provide 
information about the incident neutrino energy spectrum from the observed recoil electron 
energy spectrum. This piece of experimental observation tells us about the form of the 
energy dependence of the suppression rate which is extremely important in distinguishing 
between allowed solutions to the SNP. Also, since it is a real time experiment, it can 
divide its data set into day and night bins. Hence the detector can give information on 
the difference between the observed solar flux during day and night. 

The Kamiokande experiment, located in a deep mine at Mozumi, Japan, was a 4.5 
ktons detector with a threshold of 7.5 MeV. The observed solar neutrino flux reported by 
this experiment is 

Observed RatcKamiokande = (2.89 ± 0.42) x lO^cm-^g-^ (3.11) 
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Super-Kamiokande (SK) P is the upgraded version of Kamiokande. The first result of 
this experiment on solar neutrinos was released in 1998 [Q. The SK has now managed 
to reduce its threshold to 5.0 MeV [pT| and the observed solar fiux reported by SK after 



1258 day of data taking is [31] 

Observed RategK = (2.32 ± 0.08) x lO^cm-^g-i (3.12) 

The recoil electron energy spectrum released by the SK collaboration after 1258 day 
of data 1^ is consistent with no spectral distortion. This means that the suppression 
rate observed is essentially energy independent. The SK gives not just the total recoil 
energy spectrum but also the spectrum at day and the spectrum at night. In fig |3.3| we 
show the day and the night spectra separately for the 1258 day SK data. We can see that 
both the day and the night spectra are fiat upto Icr. The night bins have slightly more 
events than the day bins. The degree of difference between the day and night event rates 
is conveniently measured by the day-night asymmetry, defined as ^ = ($„ — ^d) /^average 
and ^average = ^{^n + ^d), whcrc $„($rf) is the obscrvcd fiux during night (day). The SK 
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A = 0.033 ± 0.022(stat.)lg:[!J^(sys.) (3.13) 

This is just a 1.3a effect which signifies almost no day-night asymmetry. 

In order to study the day-night effect in greater details, the SK divide their data 
on total rates into a day and five night bins according to the zenith angle at which the 
neutrinos arrive |Ty, ^ . They have now also divided their observed spectrum into zenith 
angle bins [^]. This helps to study the energy dependence of the suppression rate as well 
as the predicted day night asymmetry together in the most efficient manner. 

The Sun-Earth distance changes with the time of the year due to the eccentricity 
of the Earth's orbit. If the neutrinos oscillate in vacuum on their way from the Sun 
to Earth, then one should expect an extra modulation of the solar neutrino fiux due to 
oscillations with the time of the year, as the survival probability depends very crucially 
on the distance that the neutrinos travel. The SK have reported the seasonal variation of 
the solar fiux. The data is consistent with the expected annual variation due to orbital 
eccentricity of the Earth assuming no neutrino oscillations. 
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Figure 3.3: SK recoil energy spectra along with the ilcr errorbars for 1258 day data. 
The dotted errorbars are for the spectra at day while the solid errorbars are for the night 
bins. 



The Heavy Water Detector (SNO) 

The Sudbury Neutrino Observatory at Sudbury, Canada is the worlds first heavy water 
detector containing 1 kton of pure D2O surrounded by 7 kton of pure H2O. The main 
detection process is the charged current (CC) breakup of the deuteron 

Ue + d-^p + p + e' (3.14) 

It has now declared its first results on the observed flux [|lT| 

$^^0 _ ^ 75 ^ 0. 07 {stat)toil{sys) x lO^cm-^s"^ (3.15) 

The detector threshold for the kinetic energy of the observed electron released in (|3.14|) is 



6.75 MeV. SNO has also released its observed flux measured by the electron scattering 



(ES) reaction (cf. eq. ( 3.1C )) and they report |11 



$1^^^ = 2.39 ± OM{stat)ro [l{sys) x lO^cm-'s"' (3.16) 
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which agrees with the SK observation ( p.l2| ), though the errors are still large. Apart from 
the total rates, SNO also gives the recoil electron energy spectrum for the CC events 
and they do not report any significant distortion with energy. However the real strength 
of this detector is its ability to measure the flux of all the neutrino species with equal 
cross-section via the neutral current (NC) breakup of deuteron 

u^ + d—^n + p + jy^ (3-17) 

The observed NC rate from SNO is being eagerly awaited. Implications of the NC rate 
for the mass and mixing parameters is discussed in detail in . 

Summary of the Experimental Results 

In order to summarize the main results available from all the solar neutrino exper- 
iments, we present in Table |3.2| the ratio of the observed to the expected total rates^ 
in the Ga, CI, SK and SNO experiments^. The corresponding rough estimates for the 
compositions of the observed flux is also shown. Since the ES rate in SK and SNO is 
sensitive to both as well as u^/ u^-, we show in brackets separately the Uf, contribution to 
the observed rate assuming — tractive oscillations. We note that the observed rate have a 
strong nonmonotonic dependence on the neutrino energy since the Ga experiments which 
see the lowest energy z/g have the highest rate, the CI experiment observes intermediate 
energy neutrinos and reports the lowest rate, while the SK and SNO which are sensitive 
to the highest energy neutrinos, have a rate that is intermediate between the Ga and CI 
rates. 

In sharp contrast to the strong nonmonotonic energy dependence exhibited by the 
data on total rates, the recent SK data on the energy spectra at day and night show 
no evidence for any energy dependence. The SNO is also consistent with no spectral 
distortion, however the errorbars for the SNO spectrum is still high. Hence there is an 
apparent conflict between the total rates and the SK spectrum data since the former would 
prefer solutions with strong nonmonotonic energy dependence while the latter would favor 
•^From now onwards we shall call these ratios as the observed rates. 

^ We choose to neglect the Kamiokande data since the Kamiokande observations are consistent with 
the SK data which has much higher statistics. 
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pvnpriTTipnl" 


obsvd 
BPBOO 




CI 


0.335 ± 0.029 


(75%), '^Be (15%) 


Ga 


0.584 ± 0.039 


pp (55%), (25%), ^B (10%) 


SK 


0.459 ± 0.017 
(0.351 ±0.017) 


^B (100%) 


SNO(CC) 


0.347 ± 0.027 


^B (100%) 


SNO(ES) 


0.473 ±0.074 
(0.368 ± 0.074) 


^B (100%) 



Table 3.2: The ratio of the observed solar neutrino rates to the corresponding BPBOO 
SSM predictions. The rate due to Ve events in SK and SNO(ES) assuming z/e — t'active 
oscillations is shown within parentheses. The Ga rate corresponds to the combined SAGE 
and GALLEX±GNO data. Also shown is the composition of the observed fluxes. 

solutions with relatively weak dependence on energy. In addition the SK data is consistent 
with little or no day-night asymmetry which we shall see rules out large parts of the 
parameter space which predict strong Earth matter effects. 

3.1.3 The Solar Neutrino Code 

We perform a dedicated analysis of the global solar neutrino data on the total observed 
rate and the SK day-night recoil electron energy spectrum. This takes into account all 
available independent experimental features of the solar neutrino data. We take the rates 
from the CI, Ga (SAGE and GALLEX±GNO combined), SK and SNO CC experiments. 
The SNO ES data is not incorporated as it has large error. We also leave out the SNO CC 
spectrum for the same reason. We do not incorporate the Kamiokande rate as discussed 
before. We use the minimization technique to determine the best-fit parameters and 
draw the C.L. contours. For the statistical analysis for the total rates we define the 
function as 

xi=j: - (4)"' {K - R7) (3-18) 

where Rf' is the theoretical prediction of the event rate for the i^^ experiment and i?^^^ 
is the corresponding observed value shown in Table p.2| . The error matrix af^ contains 
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the experimental errors, the theoretical errors and their correlations. For the evaluation 
of the error matrix afj we have followed the procedure given in ||33|| . 

The expected event rate for the radiochemical experiments CI and Ga in presence of 
oscillations is 

8 

= E / MEMiE,){PUE,))dE, (3.19) 

where ai{Ey) is the capture cross section for the i^^ detector, E^^^ is the detector threshold, 
{P(.e{Ey)) is the neutrino survival probability averaged over the distribution of the neutrino 
production region inside the Sun, (j)k{E^) is the neutrino spectrum from the k*^^ source 
inside the Sun and the sum is over all the eight sources. For the SK experiment the 
corresponding event rate is given by 



R'sK=j^JEA I dETR{EA,ET) I dE,K{Eu 



{PeeiE,)) + —^{PeAE,)) 



dEr' ' " dEr 



(3.20) 



where At,^ is the normalized neutrino spectrum, Et is the true and Ea the appar- 
ent (measured) kinetic energy of the recoil electrons, i?^ is the detector threshold energy 
which is 5.0 MeV and Ii{EA,Erp) is the energy resolution function which is taken as 



In eq. ( p.20|) {Pee) is the time averaged Ue survival probability, (Pex) is the time averaged 
transition probability from Ue to z/^., where Ux is either or Ut, da^^/dEx is the differential 
cross section for (z/g — e) scattering while da^^/dEx is the corresponding cross section for 
{vx ~ ^) scattering. Note that if one has Ve ~ t'steriie transitions involved the second term 
will be absent and only the Ve contribution to the scattering rate will survive. 
For the Ve — d CC event rate in SNO we use 

r,th _ S dE,KXEu)(yCc{Eu){Pee{E,)) 

I dKX.SE.)acciE^) ^^-'^^ 

acc = f dEA rdETR{EA,ET)^^^^^^^^^^ (3.24) 
Je'^ Jo dEx 

For SNO E^J^ = (6.75 + rrie) MeV, where rUe is the mass of the electron and da,y^d/dET 

is the differential cross section of the Ve ~ d interaction. One of the major uncertainties 
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in the SNO CC measurement stems from the uncertainty in the u,, — d cross-section. We 



use the cross-sections from [p5| which are in agreement with p^. Both calculations give 



an uncertainty of 3% which is also the value quoted in [0]n- R{,Ea, Et) for SNO is given 



by the same functional form (|3.21| ) with the a = asNO given as |jTT| 

asNo = (-0.462 + 0. 547 V^ + 0.008722Et) (3.25) 



For the analysis of the day-night effect and the energy behavior of the suppression 
rate we define a function for the SK 1258 day day-night recoil electron energy spectra 
as 
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Xl = E (^n^f - {<ylr' {x^Sf - Sf') (3.26) 

where Sf^ are the theoretically calculated predictions for the i^^ energy bin, normalized 
to BPBOO, 5*^^^ are the corresponding observed values and the sum is over 19 day + 19 
night energy bins provided by SK. The error matrix for the spectrum analysis is defined 
as in |^8||. In eq. ( p.26| ) X„ is an overall normalization constant which is allowed to 
vary freely in the analysis. The SK provides information about three aspects of the solar 
neutrino flux suppression, (i) the overall suppression rate, (ii) the energy dependence of the 
suppression and (iii) the effect of Earth matter on the suppression rate. The information 
about the overall flux observed by SK is embodied both in the total rate and in the 
spectrum data. Since we have already accounted for this piece of information in the 
xl; we avoid the double counting of the total suppression rate in x| by introducing this 
floating normalization X„. Thus the SK day-night spectrum data provides information 
on only the presence of energy distortion, if any. It gives information on the the day-night 
asymmetry as well. 

For the global analysis we take into account the data on total rates as well the SK 

day-night spectrum data and define our total as = X/; + x|- If we assume no 

new property for the neutrino and use the flux predictions from BPBOO, then the value 

of = 89.27 which is definitely unacceptable. Even if the constraints on the solar 

^It was recently pointed out in |Q that the calculation of both and |Q underestimate the total 
Ve — d cross-section by 6%. We have not included this effect in our calculation. 
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models is relaxed, so that one allows the fluxes to take on any arbitrary value subject 
to the solar luminosity constraint, the fit is extremely poor if all the three experiments 
are considered together. The data cannot be explained by this approach, even if one 
takes only two experiments at a time. In fact as discussed in the introduction, all such 
fits predict "missing ^Be neutrinos" . This happens because the Ga observed flux can be 
almost accounted for by the pp and pep fluxes alone, given the luminosity constraint. If 
simultaneously the observations of the water Cerenkov experiments are to be accounted 
for, then there is an extra contribution from the flux in Ga leaving no room for the 
^Be flux. If on the other hand one considers CI and SK together, then the expected ^B flux 
in the former from the observation of the ^B flux in the latter, more than compensates 
the observed rate in CI, again demanding complete suppression of ^Be |^9|. With the 



advent of the SNO CC result the astrophysical solution gets comprehensively ruled out 



|T^. Thus one has to invoke some new property for the neutrinos beyond the standard 
model of particle physics in order to solve the solar neutrino problem. We probe the 
viability of neutrino mass and flavor mixing as a possible explanation of this discrepancy. 

We first find the best-fit solution to the data on only the total rates by minimizing Xr- 
Next we take into account the global data on rates as well as the SK day-night spectrum 
data so that our total is = + We minimize this for ^'e — '^active oscillations 
keeping the ^B flux normalization in the total rates fixed at the SSM prediction. We 
repeat the entire analysis for Ue — z^stcriic oscillations. For both these neutrino flavor 
mixing analyses we adopt a unified approach to which we turn our attention next. 

3.1.4 Unified Formalism for Analysis of Solar Data 

The general expression for the probability amplitude of survival for an electron neu- 
trino produced in the deep interior of the Sun, for two neutrino flavors, is given by 



^ee — ^el^ll ^le + ^e2^22 ^2e \O.ZI ) 

where Afi^{k = 1,2) gives the probability amplitude of z/g transition at the solar 

surface, A^l'^ is the survival amplitude from the solar surface to the surface of the Earth 
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and A®g denotes the — ^ z/g transition amplitudes inside the Earth. We can express 

A% = a%e-^^^^ (3.28) 



where 0f is the phase picked up by the neutrinos on their way from the production point 
in the central regions to the surface of the Sun and 

«®i' = ^ + (^-^j)cos2e„ (3.29) 

where 9m is the mixing angle at the production point of the neutrino and is given by eq. 
( p. 34 ) for transitions to active and by eq. ( |2.51| ) for transitions to sterile neutrinos, Pj is 



the non-adiabatic jump probability given by eq. ( |2.47| ) which for the exponential density 



profile of the Sun can be conveniently expressed as 



Pj = exp(-7csin^^)-exp(-7e) ^33^^ 
1 - exp(-7j 



Am 

The survival amplitude A'll'^ is given by 



d InUe 



dr 



(3.31) 



Alf = e-^^^(^-^'s) (3.32) 

where Ei^ is the energy of the state z/^, L is the distance between the center of the Sun 
and Earth and -R0 is the solar radius. For a two-slab model of the Earth — a mantle and 
core with constant densities of 4.5 and 11.5 gm cm~^ respectively, the expression for Af^ 
can be written as (assuming the flavor states to be continuous across the boundaries) 

= E f^eYe-^^^V.^^f/.V^'^^'f/gf/^^e-^f/^Jt/.^ (3.33) 

a,/3,(T 

where k) denotes mass eigenstates and (a, (3, a) denotes flavor eigenstates, U, f/*^ and 
U'" are the mixing matrices in vacuum, in the mantle and the core respectively and ip'^ 
and ip^ are the corresponding phases picked up by the neutrinos as they travel through 
the mantle and the core of the Earth. The z/g survival probability is given by 



I A 



ee 



„©2| /.e |2 I ^02| 4© 12 

"el l^lel + "e2 l^2el 

+2a®a®i?e[«*e^(^^-^^)(^-^0)e^(<^2''-<^?)] (3.34) 
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Identifying = a% and = l^feP eq. (lOP can be expressed as HOl H3l Hi] 



Pee = P0Pe + (l-P0)(l-Ps 



+2^P0(1 - P0)Pe(l - P®) cose (3.35) 

where we have combined all the phases involved in the Sun, vacuum and inside Earth 
in ^. This is the most general expression for survival probability for the unified analysis 
of solar neutrino data. Depending on the value of I\vp? jE one recovers the well known 
Mikheyev-Smirnov-Wolfenstein (MSW) ||16|, |l7l and vacuum oscillation (VO) ^ limits: 

• In the regime Am^/E ~ 5 X 10 eV^/MeV matter effects inside the Sun suppress 
flavor transitions and 6^ ~ 7r/2. Therefore, from ( p.29| ), we obtain Pq ^ Pj ^ cos^ 6 as 
the propagation of neutrinos is extremely non-adiabatic and likewise, P® = cos^ 6 to give 

p;;^ = 1 - sin^ 2^sin2(Am2 (L - Rq)/AE) (3.36) 

• For Aim? / E ~ 10^^ eV^/MeV, the total oscillation phase becomes very large and the 
cos^ term in eq. ( |3.35| ) averages out to zero. One then recovers the usual MSW survival 
probability 

pMsw ^ p^p^ ^ ^1 _ p^)(i _ (3 

• In between the pure vacuum oscillation regime where the matter effects can be 
safely neglected, and the pure MSW zone where the coherence effects due to the phase 
e can be conveniently disregarded, is a region where both effects can contribute. For 
5 X 10~^° eVVMeV ~ An? / E ~ lO"*^ eVVMeV, both matter effects inside the Sun and 
coherent oscillation effects in the vacuum become important. This is the quasi vacuum 
oscillation (QVO) regime [^, In this region, Pq Pj and P® = cos^ ^ and the 
survival probability is given by ^ 



Pee = Pj cos^ ^ + (1 - Pj) sin^ 9 + sin^ 2e^Pj{l-Pj) cos^ (3.38) 
We calculate Pj in this region using the prescription given in [HH] . 
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Day-Night Effect 



For the the range of /\m? / E for which matter effects inside the Earth are important 
(the pure MSW regime), one expects a significant day-night asymmetry. During day time 
the neutrinos do not cross the Earth and is simply the projection of the Vi state onto 
the Vf. state. Hence the Vf, survival probability during day is simply 



pg = \ + {\- p.j) cos2eco^2er, 



(3.39) 



D 

ee 



The probability during night is given by the full expression ( |3.37|) . If one factors out P^ 
in the complete expression (|3.37|) which includes the Earth matter effects then^ 

(3.40) 



P., 



j,D , (2Pi^-l)(sin^g-Al) 



COS 26 

where P2e = 1 — P®- From eq. ( |3.39| ) and ( |3.40| ) we see that the extra contribution coming 
due to the matter effects inside the Earth is 



P. e - P 



D 



(2Pi^-l)(sin^^-P®) 



(3.41) 



cos 2^ 

This is the total regeneration of Ve inside the Earth which we shall call Re- For convenience 
we shall define the regeneration factor 



fv, 



P2® - sin^ 9 



In the absence of any Earth matter effects, P2® = sin^ 6 and so f] 
Pg^ from eq. ( g^ . We also note that (with P® = |a®n 

(2Pe^ - 1) 



reg 



(3.42) 



and we get back 



cos 2d 



;i-2Pj) cos 20^ 
-{P%-Pfi) 



(3.43) 



In other words the above factor quantifies the amount of level crossing due to loss of 
adiabaticity at resonance. The Earth regeneration now can be conveniently expressed as 



(Pi - P^i)i 



reg 



1 -2Pj)cos2Mr, 



(3.44) 
(3.45) 



^Note that for the purpose of simphcity of presentation we show the expressions where the phases 
have been averaged out to zero. However for the actual calculation of the probabilities we use the unified 
expression ( |3.35 ). 
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Hence from eq. ( 3.44 ) and ( 3.45| ) we see that the z/g regeneration inside the Earth depends 
on 



1. The adiahatic factor Pj: Re is maximum for Pj = 0, decreases with increasing Pj, 
hits the minimum {Re = 0) for Pj = 1/2 and changes sign for Pj > 1/2. Which 
means that for Pj > 1/2 we have further depletion of Ve as they pass through the 
Earth. 



2. The value of cos 29 m-' We can see from eq. ( p.34|) that cos26'm > if the resonance 



density is more than the density at which the neutrinos are produced, cos 26^ = if 
the neutrinos are produced at the position of resonance and cos2^m < if the 
resonance density becomes less than the production density. As the resonance 
density decreases cos26'm decreases and reaches the value cos26'm — 1. As the 
value of Am^ decreases the position of resonance for the solar neutrinos shifts fur- 
ther outward and the value of cos2^m approaches —1. Thus for the LOW solution 
cos2^m = — 1 for almost all neutrinos energies and one gets maximum regeneration 
for the LOW solution. 

3. The regeneration factor: The net regeneration due to Earth matter effects depends 
crucially on the value of f^^g which determines quantitatively the actual effect of 
Earth matter. 



In fig 3.4 we show the Earth regeneration Pr as a function of . for three values 
I — ^ ° Am 

of the mixing angle in the range < 9 < ^. We note that the Earth matter effects are 
important only for 10^^ eV~^ ~ ^ ^^^"^ eV~^ which is the "pure" MSW regime 

and peaks at ~ 3 x 10^^ eV~^. Since around P = 5 — 15 MeV the SK day-night 

data allows for very small day-night asymmetry, most of regions around (2 — 8) x 10^^ 
eV^ for large mixing angles are disfavored. In fig |3.5| we show the regeneration factor 



freg and the total Earth regeneration Re vs energy at the SK latitude for typical values 
of the parameters in the SMA, LMA and LOW-QVO regimes. Since the latitude of the 
other detectors are not very different we do not expect freg and P^; to be very different 
for them. Noteworthy point is that while the regeneration factor freg is positive for all the 
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Figure 3.4: Earth reereneration of function of . ^ .) for three different values of 

the mixing angle. 



three cases considered, Re turns out to be negative for the SMA case. This is because 
for the SMA solution Pj > 1/2, or in other words — Pg® < 1, signifying large level 
crossing from the z/2 to the ui state in the solar matter at resonance. On the other hand 
for the LMA and LOW solutions the neutrino moves adiabatically inside the Sun, the 
z/g produced in the 1^2 state remains in a 1^2 state throughout and Pj — Pg® ~ 1, so that 
Re ^ frcg. Thus for both LMA and LOW solutions one has positive regeneration of z/g 
inside the Earth, the effect being more for the latter since for low Am^ all the neutrinos 
resonate far away from the production zone and cos2^m is closer to -1. Also note that 
for the LOW solution the regeneration is important at low energies while LMA has more 
regeneration for higher energy neutrinos. 



We finally present in fig p.6| the actual survival probability Pgg vs energy during day 
(shown by dotted lines), during night (shown by dashed lines) and the day-night average 
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Figure 3.5: The regeneration factor f^eg and the net regeneration Re as a function of 
neutrino energy for typical values of the parameters in the SMA, LMA and LOW regions. 



(shown by solid lines) for SMA, LMA and LOW case. In order to understand the nature 
of the probabilities we call cos 26' = e and note that: 



For the SMA region e ~ 1 and from fig. 3.5 we observe that freg is very small except- 



ing for two peaks at E ^ 6 MeV and E ^ 15 MeV corresponding to strong enhancement 
of the earth regeneration effect for the neutrinos passing through the core . Hence 



Pr " ^ Pe": (3.46) 

In this region for low energy (pp) neutrinos, resonance is not encountered (resonance 
density ^ maximum solar density) and hence Pj ~ and cos 29m ~ 1 giving P^J^^^ ~ 1 . 
For intermediate energy ijBe) neutrinos cos26'm —1 (resonance density ^ production 
density) and P^J^^ ^ Pj ^ for these energies. For high energy (*P) neutrinos also, 
cos2^m ~ — 1 and P^/^^ ~ Pj, with Pj rising with energy. 
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• For the LMA solution the motion of the neutrino in the solar matter is adiabatic 
for almost all neutrino energies and Pj ~ 0. For low energy neutrinos the matter effects 
are weak both inside the Sun and in Earth giving freg ~ and cos 26m ~ e so that for Ga 
energies 



Pr-'^^a + a (3.47) 
At SK and SNO energies matter effects result in cos2^^.m — 1 while freg is small but 



non-zero (~ 0.03 at 10 MeV as seen from fig. 3^) giving 



1 



2 



I reg 



= sin' 9 + freg (3.48) 

• In the LOW region cos2ftm ~ — 1 for all neutrino energies and Pj ~ (except for 
very high energy neutrinos) and thus for all neutrino energies 

Pef = - e) + freg (3.49) 

where freg is small for high energy neutrinos and large for low energy neutrinos (cf. fig. 



3.1.5 Results and Discussions 



We present in Table 3.3 the results of the x analysis for — inactive oscillations, using 



data from the CI, Ga, SK and SNOQ experiments. We use the total rates given in Table 



3^ and the 1258 day SK recoil electron energy spectrum at day and night. We show the 
best-fit values of the parameters Am^ and tan^ 9 , xLin ^^e goodness of fit (GOF) for 
the SMA, LMA, LOW-QVO, VO and Just So^ || solutions. 



The best-fit for the only rates analysis comes in the VO region which is favored at 
28.79%. Prior to SNO the SMA solution could explain the nonmonotonic energy depen- 
dence of the survival probability from the CI, Ga and SK experiments well and was the 
best-fit solution. But with the advent of SNO it falls into disfavor and is allowed at only 
6.59%. For the LMA solution on the other hand the survival probability is given by eqs. 

■"We incorporate only the SNO CC rate as the SNO ES rate and the SNO CC spectrum still have 
large errors. 
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a. 




Energy (MeV) 

Figure 3.6: The neutrino survival probability Pge as a function of neutrino energy for the 
LMA, SMA and LOW solutions. 



( p.47| ) and ( |3.48|) at Ga and SK/SNO energies respectively and for the values of e from 



Table 2]3| and /reg given in fig. ^]5|, it approximately reproduces the rates of Table . 
LMA is allowed at 18.27% while LOW-QVO is barely allowed at 1.55%. In fact the LOW 
solution gets allowed only due to the strong Earth regeneration effects at low energies 
which helps the LOW solution to explain the Ga data better. In addition to these four 
solutions we have a fifth solution called the Just So^ solution at Am^ ~ 5.38 x 10^^^ 



eV^. For these Am^ one gets a very small survival probability for the ^ Be neutrinos 
while for the neutrinos the survival probability is close to 1.0 |]SD|. Therefore it cannot 
explain the total rates data. 

We next perform a complete global analysis of the solar neutrino data taking the four 
total rates and the 1258 day SK day-night recoil electron spectrum data. We present 
in Table |3.3| the results obtained by minimizing = xl; + x|- We get five allowed 
solutions LMA, VO, LOW, SMA and Just So^ in order of decreasing GOF. The LMA 
solution can approximately reproduce the rates as discussed above and since the survival 
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Nature of 






tan2 6' 


• 

Arrwn 


Goodness 




Solution 


in eV^ 








of fit 




SMA 


7.71 X 10- 


-6 


1.44 X 10-3 


5.44 


6.59% 




LMA 


2.59 X 10' 


-5 


0.34 


3.40 


18.27% 


rates 


LOW-QVO 


1.46 X 10' 


-7 


0.67 


8.34 


1.55% 




VO 


7.73 X 10- 


11 


0.27 


2.49 


28.79% 




Just So^ 


5.38 X 10- 


12 


1.29 


19.26 


6.57 X 10-3% 




SMA 


5.28 X 10' 


-D 


3.75 X 10"^ 


51.14 


9.22% 


rates 


LMA 


4.70 X 10' 


-5 


0.38 


33.42 


72.18% 


+ 


LOW-QVO 


1.76 X 10' 


-7 


0.67 


39.00 


46.99% 


spectrum 


VO 


4.64 X 10- 


10 


0.57 


38.28 


50.25% 




Just So^ 


5.37 X 10" 


12 


0.77 


51.90 


8.10% 



Table 3.3: The best-fit values of the parameters, Xmin ' goodness of fit from 

the global analysis of rates and rates-|-spectrum data for MSW oscillations involving two 
active neutrino fiavors. 

probability ( |3.48| ) for SK is approximately energy independent it can account for the flat 
recoil electron energy spectrum. LMA thus gives the best-flt being allowed at 72.18%. 
The LOW solution with best-flt e = 0.2 and freg ~ 0.2 for Ga and ~ 0.025 for SK energies 
(cf. flg. |3.5| ) can just about reconcile the Ga and SK rates. However it provides a very 
good description of the flat SK spectrum and is allowed at 46.99%. The VO solution at 
Am^ ~ 4.64 X 10~^° eV^ gives a very low for the spectrum data and hence the overall 
flt for VO is very good. However for the SMA solution there is a mismatch between the 
parameters that give the minimum for the rates data and the spectrum data. The 
spectrum data prefers value of tan^ 6 which are one order of magnitude lower than those 
preferred by the rates data. Thus the overall flt in the SMA region suffers and it is allowed 
at only 9.22%. The Just So^ solution is very bad for the rates but since it gives a flat 
probability for the neutrinos the spectrum shape can be accounted for and the global 
analysis gives a GOF of 8.1%. 

In flg. ^]7|we show the 90%, 95%, 99% and 99.73% C.L. allowed areas from the analysis 
of the data on total rates (shown in the left hand panel) and the combined data on rates 
and the SK spectrum (shown in the right hand panel). For the only rates case we have 
allowed areas in the SMA, LMA, LOW-QVO and VO regions. For the global analysis we 
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Figure 3.7: The 90%, 95%, 99% and 99.73% C.L. allowed areas from the analysis of 
the total rates (left panel) and global analysis of the rates and the 1258 day SK recoil 
electron spectrum at day and night (right panel) , assuming MS W conversions to sequential 
neutrinos. 



get allowed zones in the LMA and the LOW-QVO zones. In the VO region we get just 
two small areas which are allowed. But the most significant feature is the disappearance 
of the SMA solution from the global fit even at 99.73% C.L. (3a). 

We have repeated the entire analysis for the Ug — z/gteriie case and have reported the 
results of the fit in Table |3.4| . We find that after the inclusion of the SNO data, all the 
solutions are disfavored with a probability of more than 99% from the total rates analysis 
while for the global analysis the GOF of these become much worse. 

3.2 The Atmospheric Neutrino Anomaly 

The atmospheric neutrinos are produced due to the collision of cosmic rays (Ncr) with 
the nuclei in the atmosphere [Nair) resulting in a chain of reactions which culminates 
in the production of neutrinos and antineutrinos with expected flavor ratio of roughly 
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Nature of 




tan2 6' 


y2 . 


Goodness 




Solution 


in eV^ 






of fit 




SMA 


4.18 X 10-6 


5.72 X 10-^ 


17.24 


1.80 X 10-2 % 




LMA 


4.98 X 10-5 


0.54 


23.96 


6.27x10-^% 


rates 


LOW-QVO 


1.00 X 10-^ 


0.94 


24.26 


5.40x10-^% 




VO 


1.07 X 10-1° 


0.27 


15.71 


3.88 X 10-2% 




Just So^ 


5.37 X 10-12 


1.28 


19.40 


6.13 X 10-3% 




SMA 


5.59 xlO-*^ 


2.83 X 10-^ 


54.21 


5.35% 


rates 


LMA 


6.13 X 10-5 


0.50 


52.93 


6.75% 


+ 


LOW-QVO 


2.93 X 10-s 


1.00 


53.18 


6.45% 


spectrum 


VO 


4.67 X 10-1° 


0.37 


46.28 


19.70% 




Just So^ 


5.37 X 10-12 


0.77 


52.09 


7.83% 



Table 3.4: The best-fit values of the parameters, xLin ' ^^e goodness of fit from the 
global analysis of rates and rates-|-spectrum data for MSW oscillations involving t'e" '^sterile 
oscillations. 



z/e : z/^ : z/^ ^ 1 : 2 : 0. 



± 



n\ K', 



(3.50) 



(3.51) 



The first pioneering measurement of events induced by these neutrinos were made in the 
Kolar Gold Fields in India and almost simultaneously in South Africa [^. But the 
interpretation of the results from these experiments remained ambiguous due to both 
theoretical as well as experimental uncertainties. It was only since late eighties that 
people earnestly started pursuing this problem |Q 0, Q and around 1988 

Kamiokande, which was originally designed to detect proton decay, declared its results 
on atmospheric neutrino measurements which showed a marked deficit of observed to 



expected z/^/z/g ratio ||53|. This embodies the Atmospheric Neutrino Anomaly. This 
was the second evidence of missing neutrinos after the solar neutrino problem and again 
called for the existence of neutrino mass and mixing for the correct interpretation of the 
experimental results. 
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However due to the large theoretical and experimental uncertainties involved in the 
atmospheric neutrino measurement, it failed to prove itself as a compelling evidence of 
neutrino oscillations until the arrival of the Super-Kamiokande atmospheric results in 1998 



58]. The SK atmospheric neutrino data |58, M, KG] not only confirmed the suppression 



of the muon type neutrinos, it also firmly established the fact that the observed deficit 
has a zenith angle dependence. That is, the neutrinos coming at larger zenith angles and 
hence traveling distances of the order of the diameter of the Earth were suppressed more 
compared to neutrinos arriving directly from the top of the detector. This differential de- 
pletion of the atmospheric z/^ flux is referred to as the observed up-down asymmetry. The 
z/g events on the other hand are reported to be consistent with theoretical expectations. 

The most convincing particle physics scenario which can explain all aspects of the SK 
data is z/^ — Ur oscillations with Am^ ~ 10^^ eV^ and sin^ 20 ~ 1 J^l], ^ |6^. Thus 



the theory of neutrino mass and mixing, which was conjectured in the late sixties as a 

plausible solution to the solar neutrino problem, was established as an accepted reality 

by the SK atmospheric neutrino results in 1998. This observation of neutrino mass by 

the SK (although indirect) is the first and till date the only evidence of physics beyond 

the standard model of particle physics. 

The observed depletion of atmospheric neutrinos was also reported by the 1MB 

earlier and more recently by the Soudan2 and MACRO J^] collaborations!. In order 

to do justice to all these experiments one should perform a global analysis taking all 

the experimental data into account. However due to its overwhelming statistics, the SK 

data dictates the atmospheric neutrino analysis, with the other experiments having little 

impact on the fit and the allowed parameter regions. Hence considering the rich statistics 

and low systematics of the SK, we choose to work with just the SK data in our analysis of 

the atmospheric neutrino anomaly and perform dedicated -fits to the SK 1144 day data. 

We describe two different established methods of analysis of the atmospheric data and 

discuss their merits and demerits, present the best-fit Am^ and sin^ 29 and display the 

^The Frejus ||5^ and Nusex experiments though were consistent with no deficit of atmospheric 
neutrinos. 
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90% and 99% CL. allowed zones in the Am? -sir? 26 parameter space for two-generation 
^'/x ~ oscillations. 

In section 3.2.1 we give an essence of the atmospheric neutrino flux calculations and 
the theoretical uncertainties associated with it. In section 3.2.2 we discuss the various 
features of the SK atmospheric experiment and state the main experimental results. In 
section 3.2.3 we describe our atmospheric neutrino code and discuss the advantages and 
disadvantages of the two methods of analysis that we have used. Finally in section 
3.2.4 we present the results of our fits. 

3.2.1 Atmospheric Neutrino Flux Predictions 

Atmospheric neutrinos are the result of a cascade of reactions following the interaction 
of the cosmic rays with the air nuclei (cf. eq. (|3.5CI|) , ( p.51|) ). In the decay chain ( |3.51|) , 



for every electron neutrino (or antineutrino) created one has two muon type neutrinos (or 
antineutrinos)^ produced. Thus one may naively guess that the ratio of fluxes 



n 



Ve + l^e 

2 (3.52) 



But for a realistic prediction for the atmospheric neutrino fluxes one needs to fold in 
the cosmic ray fluxes with the hadronic interaction model, which takes into account the 
complete reaction chain and perform a comprehensive calculation. Two such widely used 
atmospheric flux calculations are by Honda et al.^\ and Agarwal et a/.[|68|. 



The primary cosmic rays which are the main ingredient of the atmospheric neutrino 

flux calculations have large uncertainties in their over all normalization, composition, as 

well as in their spectrum, all of which are extremely crucial in the atmospheric neutrino 

flux calculations. Though the cosmic ray fluxes are relatively well known at ~ 30 GeV 

(these give rise to the sub-GeV atmospheric neutrino fluxes), there are few experimental 

measurements for the higher energy cosmic rays {E ~ 1000 GeV which result in multi-GeV 

neutrino fluxes) and hence this regime is plagued with large uncertainties. 

^For the purpose of discussion we shall henceforth use the term neutrino to mean both neutrino as 
well as antineutrino though for the actual calculation they are treated differently. 
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The resultant atmospheric neutrino fluxes also depend crucially on the hadronic in- 
teraction model. In their calculation, Honda et a/.use a full Monte Carlo method for the 



sub-GeV neutrinos while for the multi-GeV fluxes they employ a "hybrid model" ||67 
They make their atmospheric neutrino flux calculations using a one-dimensional approx- 
imation where they assume that all the secondary particles including the mesons and the 
leptons are collinear with the incident primary cosmic ray. Though this approximation 
works well for the higher energy calculations, it fails for the low energy sub-GeV fluxes, 
particularly for the horizontal zenith angles and therefore should be replaced by a full 



three-dimensional calculation ||69| 



Let us next discuss some of the most important factors which affect the atmospheric 
neutrino flux predictions. 

1. The solar activity: The effect of the solar wind on the cosmic rays entering the solar 
sphere of influence is to deviate their trajectories away from the Earth. This affects 
the low energy cosmic rays more than the higher energy ones and is known as the 
solar modulation of the cosmic ray fluxes. Since the solar wind depends on the solar 
activity, the fluxes are maximum at the solar minimum and minimum at the solar 
maximum. 

2. The geomagnetic field: Since the cosmic rays consist of charged particles they ex- 
perience a repulsion due to the Earth's magnetic field and may be deflected away. 
Hence only particles with momentum above a certain threshold can break this bar- 
rier and enter the Earth's atmosphere. This threshold is called the rigidity cut off' 
and is defined as, rigidity=momentum/charge. This effect results in cutting off the 
lower energy cosmic rays which affects the low energy neutrino fluxes. The geo- 
magnetic field bends the cosmic rays and thus introduces a directionality into the 
atmospheric fluxes, the effect obviously being more for lower energy neutrinos. This 
effect of the geomagnetic field results in predicting more neutrino flux from the west 
than from the east and is known as the east-west effect. SK has made observations 
of this predicted east-west anisotropy and has confirmed this estimated directional 
behavior |[70[] . 
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3. The density structure of the atmosphere: Once the pions and kaons are produced in 
the air, they may either decay creating muons and neutrinos or may interact with 
other particles. Which of these processes dominate is dictated by the energy of the 
meson and decay becomes comparable to interaction when 

E 1 mc^ 

c — -r ~ — , or _E ~ (3.53) 

mc"' an cran 

where r is the rest frame lifetime of the meson of energy E and mass m, a is the 
interaction cross-section with air and n is the number density of the air nuclei. 
Since the interaction rate depends inversely on the density of the atmosphere, and 
since the cosmic rays moving along the horizontal zenith angle interact with the air 
nuclei at a higher altitude (and hence lower density), the decay probability for these 
is more compared to interaction. Hence the resultant neutrino fluxes are maximum 
along the horizontal and minimum along the vertical directions. This gives rise to 
the zenith angle variation of the neutrino flux. Again, because the interaction rate of 
the meson is determined by its energy, the zenith angle dependence of the neutrino 
flux increases with energy. 

4. Decay lifetime of the mesons: The ratio 7^ ~ 2 only for vr decays at rest for which 
the decay (|3.51| ) is complete. However above ~ 5 GeV the decay lifetime of 
the produced muon becomes so large that they fail to decay in air. This results in 
raising the value of TZ from 2. Again, as the effective atmospheric depth along the 
horizontal direction is more than that along the vertical, the decay chain ( p.51| ) has 
comparatively a better chance of completion along the horizontal direction and this 
results in endowing TZ with a zenith angle dependence. 

Though the uncertainty due to primary cosmic rays for the sub-GeV neutrinos is small, 
they are significantly modulated by the solar activity and the geomagnetic field. Being 
low in energy, they have less zenith angle dependence. The multi-GeV fluxes on the other 
hand have less dependence on solar activity and geomagnetic field, but have large zenith 
angle variation and have to contend with huge cosmic ray uncertainties. We show in fig 

the atmospheric neutrino flux predictions by Honda et a/.[0 as a function of the 



51 




Figure 3.8: The atmospheric neutrino fluxes as a function of the neutrino zenith angle 
for three different fixed energies. Shown are the absolute fluxes as well as the ratio TZ. 



zenith angle for three typical values of energies. Also shown in the lower panels are the 
corresponding values of TZ. The fluxes shown are without the geomagnetic effects so that 
they are symmetric about the horizontal and hence the sign of cos^ is not important. We 
note that both Ve as well as decrease along the vertical, displaying the zenith angle 
behavior discussed above, the effect being more for the z/g. Thus the ratio TZ gets raised 
from the often quoted value of 2 along the vertical. The above effect is seen to intensify 
with energy. 

Among the major uncertainties involved in the atmospheric flux calculations, the most 
important contribution comes from the lack of correct estimates for the cosmic ray flux 
normalization, composition and spectrum. This results in about 10% error in the sub- 
GeV fluxes and 20% error in the multi-GeV flux range. The hadronic interaction model 
brings about 10% uncertainty above 300 MeV. The other sources of uncertainties include 
the method of calculation employed and the use of one-dimensional approximation. All 
these factors sum up to result in total uncertainty of about 20-30% in the absolute values 



of the fluxes [3/]. For the ratio 7Z on the other hand the uncertainties due to the primary 
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cosmic rays and the error in the calculation scheme cancel out and one has only about 
5% error for the sub-GeV and about 10% error for the multi-GeV neutrinos . 

3.2.2 Results from the Super-Kamiokande 

The Super-Kamiokande as described before is a large water Cerenkov detector. It 
contains 50 kton of ultra pure water, divided into two concentric cylinders. The inner 
detector with a fiducial volume of 22.5 ktons, has its surface lined by 11,146 inward 
facing 50cm photomultiplier tubes and 1885 outward facing 20cm photomultiplier tubes. 
The outward facing phototubes view the outer detector volume which serves both as a 
particle shield as well as a detector. The main process by which it detects the atmospheric 
neutrinos is the charged current quasi-elastic interaction 

ui + N ^l + N' (3.54) 

where the flavor of the incident neutrino is tagged onto the flavor of the released lepton. 
As this lepton moves in water it emits a Cerenkov cone which can be viewed by the 
phototubes lining the detector walls. The SK classifies its data into 

• e type or fi type events: The single-ring events can be classified into e or fi type 
depending on whether the observed Cerenkov ring is fuzzy or sharp respectively. 

• contained events: Events for which the products of neutrino interaction are con- 
tained within the detector volume are called contained events. If the products are 
fully contained within the inner detector volume itself, then these events are termed 
fully- contained events (FC). On the other hand if the produced leptons exit the 
inner detector volume and stop in the outer detector, then these events are called 
partially- contained events (PC). Only muons are penetrating enough to be able enter 
the outer detector and so only z/^ have PC events. 

• suh-GeV and multi-GeV events: Depending on the value of the visible energy of 
the released lepton {Ems) the FC events are categorized as sub-GeV (if E^s < 1-33 
GeV) or multi-GeV (if E^is > 1.33 GeV). 
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upward- going muons: If the muons are produced in the rock below the detector and 
penetrate into the detector volume, then such events are grouped together to be 
called upward-going muons If these muons traverse the entire detector volume 



and escape out then they are termed as upward through- going muons. If on the 
other hand they stop within the detector then one has upward stopping muons. 
These upward-going muons are recorded for zenith angle > 90" only. The neutrino 
induced muons with zenith angle less than this are not taken into account due the 
huge contamination resulting from the muons in the cosmic ray fluxes themselves. 

zenith angle bins: At the detector, the neutrino flux come from all directions. Thus, 
the total path length between the production point in the atmosphere and the 
detector varies from about 10 km to 13,000 km depending on the zenith angle. 
Neutrinos with zenith angle less than 90° {downward neutrinos) travel a distance of 
~ 10 - 100 km from their production point in the atmosphere to the detector while 
the neutrinos with larger zenith angles [upward neutrinos) cross a distance of up to 
~ 13,000 km to reach the detector. The SK can see the zenith angle of the lepton 
that emits the Cerenkov cone. The lepton zenith angle B is related to the incoming 
neutrino zenith angle ^ through the relation 

cos 9 = cos ^ cos + sin ^ cos sin ip (3.55) 

where i/) is the angle between the incoming neutrino ui and the scattered lepton / 
and is the azimuthal angle corresponding to the incident neutrino direction. The 
SK divides its events into bins corresponding to the zenith angle of the observed 
lepton cosO. 



After 1144 day of data taking the SK divide their observed number of sub-GeV and 
multi-GeV e type and /i type events into 10 zenith angle bins each. Instead of showing 
the actual number of events reported by SK, for the sake of clarity we present in fig 



3T9| the ratio of the number of observed (A^) to the number of expected (A"o) events. 



for the sub-GeV e-type (SGc ), multi-GeV e-type (MGc ), sub-GeV /i-type (SG^ ) and 
multi-GeV /x-type (MG^ ) data samples. We have included the PC events along with 
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the FC MG^ events in the panel for MG^ . Henceforth we will follow this procedure for 
representation of the PC and MG^ data samples. The ratio N/Nq are shown for the ten 
zenith angle bins along with the ilcr errorbars. Had there been no anomaly this ratio 
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Figure 3.9: The sub-GeV and multi-GeV event distributions vs. zenith angle. N is the 

number of observed e type or /i type events and Nq is the Monte Carlo expectation. The 
errorbars show the ±1(7 uncertainties. Also shown by the dashed lines are the zenith 
angle distributions for the various data samples predicted by the two-generation i/^ — Uj. 
oscillation best-fit solution. 



would have been one for all the data bins. However we observe that 

1. The e-type events, both SGe and MGe , are consistent with theoretical expectations. 

2. The observed /x-type events on the other hand show depletion compared to expec- 
tations. 

3. The deficit of /^-type events has a zenith angle distribution. The 1/^ depletion being 
most in the first zenith angle bin (upward) and least in the last one (downward), 
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for both the SG^ and MG^ data samples. Thus the data has a marked up-down 
asymmetry. 

4. In addition the deficit has an energy dependence. The up-down asymmetry being 
more for the multi-GeV than for the sub-GeV neutrinos. 

Since the uncertainties in the absolute fluxes are immense (~ 30%), one may instead use 
the double ratio, 

R = (3.56) 

where (A^^/A^e)|obs are the ratio of the total number of observed fi-tjpe to e-type events 
in the detector and {N^/Ne)\MC is the corresponding Monte Carlo expectation. Differ- 
ent atmospheric flux calculations agree to within better than 5% on the magnitude of 
this quantity. The up-down asymmetry can be then expressed in terms of an up-down 
asymmetry parameter Yi defined as 



Y ^ l^'l )\obs ryN 

-'^ — / Ar-0.2 /Ar+0.2v 1,0.0/ J 



Here Nj^^'"^ denotes the number of /-type events produced in the detector with zenith 
angle cos© < —0.2, i.e. the upward neutrino events while Nj^^''^ denotes the number of 
Z-type events for cosO > 0.2 i.e. events coming from downward neutrinos. The central 
bin has contributions from both upward and downward neutrinos and is not useful for 
studying the up-down asymmetry. For the 848 day and 535 day data the R and Yj are 
given in Table p.5| . The Yg is close to 1 but the Y^ and R are less than unity for both the 
sub-GeV as well multi-GeV cases. 



Quantity 


848 day data 


535 day data 


Sub-GeV 


Multi-GeV 


Sub-GeV 


Multi-GeV 


R 


0.69 ±0.05 


0.68 ±0.09 


0.63 ±0.06 


0.65 ±0.09 




0.74 ±0.04 


0.53 ±0.05 


0.76 ±0.05 


0.55 ±0.06 


Ye 


1.03 ±0.06 


0.95 ±0.11 


1.14 ±0.08 


0.91 ±0.13 



Table 3.5: The double ratio R and up-down asymmetry parameter YJ. 
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3.2.3 The Atmospheric Neutrino Code 

The expected number of / (e or /i) Hke 1 ring events recorded in the detector in presence 
of oscillations is given by 

Ni = TiT dE / dEi / dcostp / dcos^ — / d(j) 

Jo "'{-BOmin J-i J-i Svr Jo 

^ giEJ^,f^pE^ (3.58) 

dE d cos ^ dEi dcosijj 

where ut denotes the number of target nucleons, E is the neutrino energy, Ei is the 
energy of the final charged lepton, ip is the angle between the incoming neutrino z// and 
the scattered lepton /, ^ is the zenith angle of the neutrino and is the azimuthal angle 
corresponding to the incident neutrino direction (the azimuthal angle relative to the ip 
has been integrated out). The zenith angle of the charged lepton in terms of ip, ^ and (p 
is given by eq. ( p. 551 ). d"^ Fi /dEd cos is the differential flux of atmospheric neutrinos of 
type ui, d^cTi/ dEid cos ip is the differential cross section for viN IX scattering and e{Ei) 
is the detection efficiency for the 1 ring events in the detector. The efficiencies that were 
available to us are not the detection efficiencies of the charged leptons but some function 



which we call e'{E) defined as |^ 



e'iE) = ''f^'l (3.59) 
Pyi^i is the survival probability of a neutrino flavor / after traveling a distance L given by, 



{Re + hf - Re^sm^i- Re cos i (3.60) 

Re being the radius of the Earth and h is the height of the atmosphere where the neutrinos 
are produced. We use the atmospheric neutrino fluxes from For the sub-GeV events 
the dominant process is the charged current quasi-elastic scattering from free or bound 



nucleons. We use the cross-sections given in |74]. The events in multi-GeV range have 
contributions coming from quasi-elastic scattering, single pion production and multi pion 
production and we have used the cross-sections given in |]75[. For the multi-GeV events 



we assume that the lepton direction O is the same as the incoming neutrino direction ^. 
But actually they are slightly different. We simulate this difference in the zenith angles 

57 



by smearing the angular distribution of the number of events with a Gaussian distribution 



having a one sigma width of 15° for fi type events and 25° for the e type events |^3[. For 
the sub-GeV events, difference in direction between the charged lepton and the neutrinos 
are exactly taken care of according to eq. (|3.58|) and ( |3.55| ). 

In order to do a statistical analysis of the data we define a function as 



M Ml I I e ^ e 



(3.61) 



where the sum is over the sub-GeV and multi-GeV cases. The experimentally observed 
rates are denoted by the superscript "exp" and the theoretical predictions for the quanti- 
ties are labeled by "th" . AR^^^ is the error in R obtained by combining the statistical and 
systematic errors in quadrature. AV^^^ corresponds to the error in Y. For this we take 
only the statistical errors since these are much larger compared to the systematic errors. 
We include both the e-like and the /i-like up-down asymmetries in the fit so that we have 
4 degrees of freedom (6 experimental data - 2 parameters) for the oscillation analysis in 
the two parameters Am^ and sin^ 29. 



The use of these type of ratios for the x analysis test has been questioned in [76 



because the error distribution of these ratios is non-Gaussian in nature. However as it has 



been shown in ||6J] the use of the i?'s and y s as defined above is justified within the 3a 
region around the best-fit point for a high statistics experiment like SK and provides an 
alternative way of doing the x^-analysis. A comparison of the results of ||6^ with those 
obtained in |HT| shows that the best-fit points and the allowed regions obtained do 
not differ significantly in the two approaches of data fitting. The advantage of using the 
ratios is that they are relatively insensitive to the uncertainties in the neutrino fluxes and 
cross-sections as the overall normalization factor gets canceled out in the ratio. We have 
included the Ye in our analysis because to justify the z/^ — Ur oscillation scenario, it is 
necessary to check that including the data on electron events gives a low value and 
hence it is the standard practice to include these in the x^-analysis ^ |61[] . 

Even though the method of analysis defined above is acceptable, a better method 
is to use the absolute number of e or /i type events taking into account the errors and 
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their correlations properly BT], p3[ . Following this method we define our O ^3 



E (iVf -iVr)(4)''(^f -^r) (3-62) 

ij=l,40 

where the sum is over the sub-GeV and multi-GeV electron and muon bins. The ex- 
perimentally observed number of events are denoted by the superscript "exp" and the 
theoretical predictions for the quantities are labeled by "th". The element of the error 
matrix 0"^ is calculated as in including the correlations between the different bins. 



For two-generation analysis of the SK atmospheric neutrino data we employ this second 
more widely used method of analysis and present our results in the next section. 

3.2.4 Results and Discussions 

Since the e type events are consistent with the Monte Carlo expectations, the only 
possible two-generation oscillation schemes that one should consider to be responsible for 
the deficit are the — and z/^ — z/gteriie oscillation modesQ. The two cases differ 
from each other in the Earth matter effects on the mass and mixing parameters for the 
upward neutrinos. While for — oscillations there are no effects of the Earth matter, 
the — z/gteriie has its Am^ and sin^ 1Q modulated in presence of matter since the sterile 
neutrinos do not have any interaction with the ambient matter. For the z/^ — z/gteriie case the 
mass squared difference and mixing angle in matter are given by eq. (|2.52|) and eq. (|2.53|) 
of chapter 2. Since for the atmospheric neutrino anomaly the mixing angle in vacuum 
is close to maximal, the effect of matter in this case is mostly to reduce the mixing. 
Thus the — z/gteriie oscillations predict a lower suppression of the upward neutrinos 
than demanded by the data and fail to reproduce the correct zenith angle distribution, 
the effect being more for the higher energy neutrinos. Thus the z/^ — t'steriie solution 
gets disfavored by the zenith angle data. In addition, since the sterile neutrinos do not 
induce any neutral current event in the detector while both the /i and r neutrinos do, for 
^ii — t'steriie uiixiug ouc cxpccts a rcductiou of the neutral current events while for z/^ — z/^- 
oscillations they remain unchanged. In their paper , the SK collaboration have shown 
that their neutral current sample is consistent with expectations, there is no depletion of 
^''in more general four generation schemes the Vfj, can oscihate into a mixture of Vr and I'stcriic [l62| . 
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the events and hence transition of z/^ to z^steriie are disfavored. The SK collaboration have 
given 99% CL. exclusion plots for a combined analysis of the neutral current enriched 
multi-ring events, the high energy partially contained events and the upward-going muon 
data sets, in terms of z/^ — z/gteriie oscillations. They superimpose the 99% CL. allowed 
region that they obtain from the analysis of the FC data over these exclusion plots 
and show that there is no region of overlap between the two. They thus conclude that 
pure two-generation — z/gtcriio oscillations are disfavored by the data and hence the pure 
i^fi — mixing, which can provide excellent fit to all the atmospheric data samples without 



any inconsistency is the favored alternative fT^]. The latest 1289 day data sample on the 



upward-going events are consistent with Pj. appearance at the 2a level thereby further 
disfavoring the sterile option. We therefore concentrate on the — Vr oscillation scenario. 

In the two-flavor — picture the survival probability of an initial z/^ of energy E 
after traveling a distance L in vacuum is 

Pv,v, = 1 - sm^2esm^ (nL/Xosc) (3.63) 

where 6 is the mixing angle between the two neutrino states in vacuum and Xosc is the 
oscillation wavelength defined as, 

A„„ = (2.47 km) (3,64) 

where Am^ denotes the mass squared difference between the two mass eigenstates. If we 
insert eq. (|3.63| ) in eq. ( p.58|) and minimize the function defined by (|3.62|) , then for 
the 1144 day of SK data we get the following best-fits and xLin 



• x'ii^ = 36.23, Am2 = 0.0027 eV^, sin^ 2^ = 1.0 
For 38 degrees of freedom (40 — 2) this corresponds to a goodness of fit of 55.14%. We 



present in fig. |3.10| the 90% and 99% CL. allowed area in the neutrino parameter space 
for two-generation z/^ — Ur oscillations. In fig p.9| we show with dashed lines the histograms 
for the best-fit value for two-generation z/^ — oscillations. Both the absolute number of 
events as well as the zenith angle distributions of the SK data are seen to be reproduced 
well. 
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Figure 3.10: Region of the parameter space allowed by the z/^ — oscillation solution to 
the 1144 day SK atmospheric neutrino data. Shown are the 90% and 99% CL. allowed 
zones. 



3.3 Bounds from Accelerator and Reactor Experiments 

Neutrino oscillations experiments can be broadly classified into two types: 

• Disappearance experiments: In these one starts with a neutrino beam of a definite fiavor 
Ua and looks for a reduction in the original fiux, after the beam has traveled a certain 
distance. 

• Appearance experiments: These experiments search for a new fiavor in the original 
Ua beam, after it has been made to travel a distance L. 

Both the solar neutrino problem and the atmospheric neutrino anomaly that we have 
discussed above are disappearance experiments. But both these experiments have the 
intrinsic defect of large uncertainties in model predictions of neutrino fiuxes. The neutrino 
beams in the laboratory experiments on the other hand have well controlled spectrum as 
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well as path length and hence are better suited for precise study of neutrino mass and 
mixing parameters. The laboratory accelerator and reactor experiments may operate in 
either the disappearance or appearance modes. While the appearance experiments have 
the advantage of being apparently easier, since one is looking for a new flavor absent in the 
initial beam, the disappearance experiments have the strength to probe the oscillations of 
the initial z/q, to all channels including z/^ — z/gteriie which can never be seen in an appearance 
experiments. 

In two-generations the conversion probability of z/^ to another flavor z/^ after traveling 
distance L is given by eq. (|2.11| ) 



P,^,, = sin2 2^sin2(^) 

= sin2 2esinM 1.27 — - — j (3.65) 

A being the oscillation wavelength. Sensitivity of a given experiment can be ascertained 
in terms of the minimum Am^ and sin^ 29 that the experiment can probe. For Am^ very 
large A <^ L, and sin^(^) — > |. For the limiting value of Pj^^^^ seen by the experiment 

sm2 20„„„ = 2P;_„_ (3,66) 

On the other haud ,f Arn^ ,s such that .s very small their we can replace the 

sine function by its argument. 

Pl^^ = sin' 29 1^1.27^^ j (3.67) 

Hence the minimum value of Am^ that the experiment can constraint corresponds to the 
case where sin^ 29 — » 1, that is 

Using the above argument the terrestrial experiments give exclusion plots in the Am^ - 
sin^ 29 plane. All accelerator and reactor experiment till date barring the LSND in Los 
Alamos have failed to observe any positive signal of neutrino oscillations^ and give ex- 
clusion plots in the Am^ -sin^ 29 plane. We give below a list of the most prominent 

experiments and their main features. 

^^The preliminary results from the K2K experiment in Japan also gives a positive indication of neutrino 
oscillations. 

62 



• Reactor Experiments 

The neutrinos in tliese experiments are produced in (3 decay of fission products in the 
reactors. Thus typically one has an initial z/g beam with energies of the order a few MeV. 
Since the energies are so small these experiments operate in the disappearance mode. 
These Pe are detected by hquid scintillation detectors through the reaction 

Pe+P^e+ + n (3.69) 

Since these are necessarily disappearance experiments, the incident neutrino flux should 
be very accurately known to reduce the systematic uncertainties. Due to these high sys- 
tematic uncertainties the reactor experiments fail to probe very low values of sin^ 26 . To 
overcome this problem to some extent and to reduce the uncertainties, most experiments 
use two or more detectors placed at different distances and what is used in the analysis 
is the ratio of the neutrino flux observed at these different detectors. The characteristics 
of the most important reactor experiments and their limits on the oscillation parameters 
are listed in Table |3]^. Among the reactor experiments the most stringent bounds in the 
low Am? range are obtained from the CHOOZ long-baseline experiment in France P^ . 
This experiment has a baseline of 1 km so that it has an average L/£'~3 (i5~3 MeV) 
and hence is sensitive to Am^ as low as 10~^ eV^. It observes no neutrino oscillations in 
the z/g disappearance mode and rules out sin 26 > 0.1 for Am2 > 7 x 10"^ eV^ |2|. 



The Bugey experiment also in France is most constraining in the higher Am? region. 
It rules out large parts of the LSND allowed zone so that only Am^ > 0.2 eV^ of the 
allowed parameter space remains valid. 

• Accelerator Experiments 
The neutrino beam in these experiments consist of mainly u^j,, u^, Ur, i^t, produced from 
decay of pions and kaons in accelerators. These experiments have much higher energy 
neutrinos and may operate either in the disappearance or the appearance mode, depending 
on the energy of the initial neutrino beam. Table ^]6| reviews the characteristics of the 
main accelerator experiments. 

The Liquid Scintillator Neutrino Detector at the Los Alamos Meson Physics Facility 
has claimed to have seen neutrino mass and mixing P5| , pp] ]. In this experiment a high 
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experiment 



probes 
oscillation 
channel 



signal 
for 
oscillation 



Am^ in eV^ 



sin^ 26 



GOSGEN p 
Krasnoyarsk 
BUGEYiT 



CHOOZ 182 



Palo Verde 18 



Vrr. 



Vx 



negative 
negative 
negative 
negative 
negative 



< 0.02 

< 0.014 

< 0.01 

< 0.002 

< 0.002 



< 0.02 

< 0.14 

< 0.2 

< 0.1 

< 0.6 



U 



CDHSW 1841 



CHARM |85 
CCFR 
E776I 



E531 1^ 
E531 18 



LSND H 



LSND g 

KARMEN2 [| 

K2K ra 



yx 

Vx 
Vx 
Vx 
Vr 
Vt 

Ve 

V.,. 



negative 
negative 
negative 
negative 
negative 
negative 
positive 
positive 
negative 
positive 



< 0.23 

< 0.29 

< 15 

< 0.075 

< 0.9 

< 9 

1.2 

< 0.007 

awaited 



< 0.02 

< 0.2 

< 0.02 

< 0.003 

< 0.004 

< 0.12 

0.003 

< 0.0021 

awaited 



Table 3.6: The reactor (I) and accelerator (II) experiments that are running/completed. 
We display the minimum value of Am^ (cf. eq. ( |3.68| )) and minimum value of sin^ 26 (cf. 
eq. (p.66|) ) that are excluded by the experiments that observe null oscillations. For LSND 
we give the best-fit Am^ and sin^ 26 . For K2K an oscillation analysis is still awaited. 



energy proton beam is made to impinge on a target thereby predominately producing 
TT"'". The small number of 7r~ created are absorbed so that the dominant decay chain is 
7r+ f^^i^fi, and /i"*" e'^Ugi?^. Thus this beam has a paucity of z/g and in their first 
experiment, the LSND made a search for i/^ — z/g oscillations by looking for i?f, appearance 
by recording the UeP —>■ e^n events. They reported a detected transition probability of 
Pu^u, = (0.3ll[5:" ± 0.05)% m. In this case the 7r+ and /i+ decay at rest (DAR). On the 
other hand if they are made to decay in flight (DIE) (the energy for this has to be much 
higher) then the /i"*" do not decay much and the flux of Vf. in the beam is small making 
it possible to probe the — oscillation channel as well. The LSND again claims to 
have seen a positive signal with a transition probability of Py^^v^ = (0.26 ± 1.0 ± 0.5)% 
pO[] . They have recently done a reanalysis of their entire data sample collected between 
1993-1998 for both 1?^ — z/g and z/^ — z/g oscillations with a common selection criteria and 
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use a novel event reconstruction |9T| . The final transition probability that they report to 
have observed is (0.264±0.067±0.045)%. The best-fit value of the oscillation parameters 



Am and sin 29 that they give are shown in the Table 3.6 



On the other hand KARMEN, the Karlsruhe-Rutherford Medium Energy Neutrino 
experiment, at the ISIS spallation neutron facility, is an experiment almost similar to 
LSND which uses i/^, and beam from 7r+ and /i"*" decays at rest and looks for — v^, 
and — Vf. oscillations. This experiment has been consistently inconsistent with the 
positive neutrino oscillation observed at LSND |9^. It gives an exclusion plot in the 
Am^ -sin^ 26' plane which excludes the entire LSND allowed region above Am^ > 2 eV^ 
and part of the zones allowed for lower Am^ . 

The other experiment which has reported to have seen a depletion of expected flux due 
to oscillations is the still running K2K long-baseline experiment in Japan. It sends a pure 
and intense z/^ beam from KEK to SK with a baseline of 250 km and is the first long- 
baseline accelerator experiment. It has an average E ^ 1.3 GeV which gives L/E ^ 200 
and has two detectors. The near one at a distance of about 300m from the source is a 
one kton water Cerenkov detector, similar in technology to the SK. The far detector, 250 
km away, is the SK itself. The experiment can work both in the z/^ disappearance as 
well as in the appearance modes. This experiment promises to probe the region of the 
parameter space allowed by the SK atmospheric data and thus holds the potential to pin 
down the exact values of the Am^ and sin^ 29 responsible for the atmospheric neutrino 
anomaly. The K2K collaboration has already released their first results on muon events. 
They detect 28 events in SK against an expectation of 37.8j;3;g which is based on the 
interactions of the incident beam at the near detector [^. Thus there is a suppression 
of the flux, though the discrepancy is still within statistical errors. Even the neutrino 
energy spectrum that they observe is as of yet very poor in statistics. Robust results and 
oscillation analyses, on both the total fluxes as well as the observed energy spectrum from 
this experiment is eagerly awaited. 



65 



BIBLIOGRAPHY 



[1] R. Davis, Jr., D.S. Harmer, K.C. Hoffman, Phys. Rev. Lett. 20, 1205 (1968). 
[2] J.N. Bahcall, N.A. Bahcall, G. Shaviv, Phys. Rev. Lett. 20, 1209 (1968). 
[3] B.T. Cleveland et ai, Astrophys. J. 496, 505 (1998). 

[4] Y. Fukuda et al, (The Kamiokande collaboration), Phys. Rev. Lett. 77, 1683 (1996). 

[5] J.N. Abdurashitov et a/., (The SAGE collaboration), Phys. Rev. Lett. 77, 4708 
(1996); Phys. Rev. C60, 055801 (1999). 

[6] W. Hampel et ai, (The Gallex collaboration), Phys. Lett. B388, 384 (1996); Phys. 
Lett. B447, 127 (1999); Talk presented in Neutrino 2000 held at Sudbury, Canada 
(T.A. Kirsten for The Gallex collaboration), Nucl. Phys. B Proc. Suppl. 77, 26 
(2001). 

[7] M. Altmann et ai, (The GNO collaboration), Phys. Lett. B492,16 (2000); Talk pre- 
sented in Neutrino 2000 held at Sudbury, Canada ( E. Belloti for the GNO Collabo- 
ration) Nucl. Phys. B Proc. Suppl. 91 44 (2001). 

[8] Y. Fukuda et ai, (The Super-Kamiokande collaboration), Phys. Rev. Lett. 81, 1158 

(1998) ; erratum 81, 4279 (1998); 

[9] Y. Fukuda et ai, (The Super-Kamiokande collaboration), Phys. Rev. Lett. 82, 2430 

(1999) . 

[10] Y. Fukuda et ai, (The Super-Kamiokande collaboration), Phys. Rev. Lett. 82, 1810 
(1999). 

[11] Q.R. Ahmad et al, (The SNO Collaboration), Phys. Rev. Lett. 87, 071301 (2001). 

[12] N. Hata, S.A. Bludman and P. Langacker, Phys. Rev. D49, 3622 (1994); S. Parke, 
Phys. Rev. D74, 839 (1995). 

[13] J.N. BahcaU, hep-ph/0108147; hep-ph/0108148. 

[14] V. Berezinsky, hep-ph/0108166. 

[15] B. Pontecorvo, Sov. Phys. JETP 6, 429 (1958); V.N. Gribov, B. Pontecorvo, Phys. 
Lett, bf B28, 493 (1969). 

[16] L. Wolfenstein, Phys. Rev. D17, 2369 (1978); Phys. Rev. D20, 2634 (1979). 

66 



[17] S. P. Mikheyev and A. Yu. Smirnov, Sov. J. Nucl. Phys. 42, 913 (1985); Nuovo 
Cimento, 9C, 17 (1986). 

[18] J.N. Bahcall, P.I. Krastev, A. Yu. Smirnov, Phys. Rev. D60, 093001 (1999); Phys. 
Rev. D58, 096016 (1998). 

[19] G.L. Fogh, E. Lisi, D. Montanino, A. Palazzo, Phys. Rev. D62, 013002 (2000); G.L. 
Fogh, E. Lisi, D. Montanino, Astropart. Phys. 9, 119 (1998); Phys. Rev. D54, 2048 
(1996); Phys. Rev. D49, 3626 (1994). 

[20] M.C. Gonzalez-Garcia, P.C. Holanda, C. Pena-Garay J.W.F. Valle, Nucl. Phys. 
B573, 3 (2000); Talk presented at Neutrino 2000, Sudbury, Canada by M.C. 
Gonzalez-Garcia, hep-ph/0009041; C. Guinti, M.C. Gonzalez-Garcia, C. Pena-Garay, 
Phys. Rev. D62, 013005 (2000). 

[21] S. Choubey, S. Goswami, K. Kar, H.M. Antia and S.M. Chitre, Phys. Rev. D, (in 
press); S. Goswami, D. Majumdar, A. Raychaudhuri, Phys. Rev. D63, 013003 (2001); 
hep-ph/9909453. 

[22] G.L. Fogh, E. Lisi, D. Montanino, A. Palazzo, Phys. Rev. D, (in press); J.N. BahcaU, 
M.C. Gonzalez-Garcia and C. Pena-Garay, JHEP, 08, 014 (2001); A. Bandyapadhyay, 
S. Choubey, S. Goswami and K. Kar, Phys. Lett. B, (in press). 

[23] J.N. BahcaU, M.P. Pinsonneault, Rev. Mod. Phys. 64, 885 (1992). 

[24] J.N. BahcaU and M.H. Pinsonneault, Rev. Mod. Phys. 64, 885 (1992). 

[25] S. Turck-Chieze and I. Lopez, Astrophys. J. 408, 347 (1993); S. Turck-Chieze, et al, 
Astrophys. J. 335, 415 (1988). 

[26] J.N. BahcaU and M.H. Pinsonneault, Rev. Mod. Phys. 67, 781 (1995). 

[27] A. Dar and G. Shaviv, TECHNION preprint, Haifa, September 1995; Nucl. Phys. 
(Proc. Suppl.) B38, 81 (1995). 

[28] J.N. BahcaU, S. Basu, M.P. Pinsonneault, Phys. Lett. B433, 1 (1998). 

[29] J.N. BahcaU, S. Basu, M.P. Pinsonneault, Astrophys. J. 555, 990 (2001). 

[30] E. Adelberger et al, Rev. Mod. Phys. 70, 1265 (1998). 

[31] Y. Fukuda, et al, (The Super-Kamiokande Collaboration), Phys. Rev. Lett. 86, 5651 
(2001); Phys. Rev. Lett. 86, 5656 (2001). 

[32] J.N. BahcaU, P.I. Krastev and A.Yu. Smirnov, JHEP, 0105, 015 (2001); Phys. Rev. 
D63, 053012 (2001); Phys. Rev. D62, 093004 (2000). 

[33] G.L. Fogh, E. Lisi, Astropart. Phys. 3, 185 (1995). 

[34] J.N. BahcaU and E. Lisi, Phys. Rev. D54, 5417 (1996). 

67 



[35] S. Nakamura, T. Sato, V. Gudkov and K. Kubodera, Phys. Rev. C63, 034617 (2001). 

[36] M. Butler, J. Chen and X. Kong, Phys. Rev. C63, 034617, (2001). 

[37] J.F. Beacom and S.J. Parke, hep-ph/0106128. 

[38] The first reference in 



[39] Last reference in |T8 



[40] G.L. Fogh, E. Lisi, D. Montanino, A. Palazzo, Phys. Rev. D62, 113004 (2000). 
[41] S.T. Petcov, Phys. Lett. B200, 373 (1988). 

[42] S.T. Petcov, Phys. Lett. B434, 321 (1998); M. Narayan, G. Rajasekharan and R. 
Sinha, Mod. Phys. Lett. A13, 1915 (1998). 

[43] S.T. Petcov, Phys. Lett. B214, 139 (1988); Phys. Lett. B406, 355 (1997); S.T. 
Petcov and J. Rich, Phys. Lett. B214, 137 (1989). 

[44] A.de. Gouvea, A. Friedland, H. Murayama, Phys. Lett. B490, 125 (2000). 

[45] A. Friedland, Phys. Rev. Lett. 85, 936 (2000); A.de Gouvea, A. Friedland, H. Mu- 
rayama, Phys. Lett. B490, 125 (2000). 

[46] G.L. Fogh, E. Lisi, D. Montanino, A. Palazzo, S.T. Petcov, Phys. Rev. D63, 093002 
(2001). 

[47] E.K. Akhmedov, Nucl. Phys. B538, 25, (1999). 

[48] M.C. Gonzalez-Garcia, C. Pena-Garay, Y. Nir and A. Yu. Smirnov, Phys. Rev. D63, 
013007 (2001). 

[49] R. S. Raghavan, Science 267, 45 (1995); P.I. Krastev and S.T. Petcov, Phys. Rev. 
D53, 1665 (1996). 



[50] See fig. 5 of the first reference of |32 



[51] C.V. Achar et al, Phys. Lett. 18, 196 (1965); Phys. Lett, bf 19, 78 (1965). 
[52] F. Reines et al, Phys. Rev. Lett. 15, 429 (1965). 

[53] K.S. Hirata et al, (The Kamiokande Collaboration), Phys. Lett. B205, 416 (1988); 
Y. Oyama et al, (The Kamiokande Collaboration), Phys. Rev. D39, 1481 (1989). 

[54] K. S. Hirata et al, Phys. Lett. B280, 146 (1992); Y. Fukuda et al, Phys. Lett. 
B335, 237 (1994). 

[55] D. Casper et al, Phys. Rev. Lett. 66, 2561 (1991); R. Becker-Szendy et al, Phys. 
Rev. D46, 3720 (1992). 

68 



[56] Ch. Berger et al, Phys. Lett. B227, 489 (1989), Phys. Lett. B245, 305 (1990); Ch. 
Bcrgcr et al, Nucl. Instrum. Methods A261, 540 (1987); K. Daum et ai, Z. Phys. 
C 66, 417 (1995). 

[57] M. Agiletta et al, Europhys. Lett. 8, 611 (1989). 

[58] Y. Fukuda et ai, (The Super-Kamiokande Collaboration), Phys. Lett. B433, 9 

(1998) ; Phys. Lett. B436, 33 (1998); Phys. Rev. Lett. 81, 1562 (1998). 

[59] Y. Fukuda et al, (The Super-Kamiokandc Collaboration), Phys. Rev. Lett. 82, 2644 

(1999) ; Y. Fukuda et al., (The Super-Kamiokande Collaboration), Phys. Lett. B467, 
185 (1999). 

[60] Y. Fukuda et al, (The Super-Kamiokande Collaboration), Phys. Rev. Lett. 82, 5194 
(1999). 

[61] G.L. Fogh, E. Lisi, A. Marrone and G. Scioscia, Phys. Rev. D59, 033001 (1999) 

[62] G.L. Fogh, E. Lisi, A. Marrone and G. Scioscia, Phys. Rev. D63, 053008 (2001). 

[63] M.C. Gonzalez-Garcia, H. Nunokawa, O.L.G. Peres, T. Stanev and J.W.F. Vallc, 
Phys. Rev. D58, 033004, (1998); M.C. Gonzalez-Garcia, H. Nunokawa, O.L.G. Peres 
and J.W.F. Valle, Nucl. Phys. B543, 3 (1999); N. Fornengo, M.C. Gonzalez-Garcia 
and J.W.F. Valle, Nucl. Phys. B580, 58 (2000). 

[64] R. Foot, R.R. Volkas and O. Yasuda, Phys. Rev. D58, 013006 (1998). 

[65] W.W.M. Allison et al, Phys. Lett. B391, 491 (1997). 

[66] Talk presented by B.C. Barish in Neutrino 2000 held at Sudbury, Canada, Nucl. 
Phys. Proc. Suppl. 91, 141 (2001). 

[67] M. Honda, T. Kajita, S. Midorikawa, and K. Kasahara, Phys. Rev. D52, 4985 (1995). 

[68] V. Agarwal, T.K. Gaisser, P. Lipari, T. Stanev, Phys. Rev. D53, 1314 (1996). 

[69] M. Honda, T. Kajita, K. Kasahara, S. Midorikawa, Phys. Rev. D64, 053011 (2001). 

[70] T. Futagami et al, (The Super-Kamiokande Collaboration), Phys. Rev. Lett. 82, 
1529 (1999). 

[71] Y. Fukuda et al, (The Super-Kamiokande Collaboration), Phys. Rev. Lett. 82, 1810 
(1999); Phys. Lett. B467, 185 (1999). 

[72] Talk presented by The Super-Kamiokande Collaboration (H. Sobel for the Collabo- 
ration) in Neutrino 2000 held at Sudbury, Canada, Nucl. Phys. Proc. Suppl. 91, 127 
(2001). 

[73] O. Yasuda, private communication. 



69 



[74] T. K. Gaisser and J. S. O'Connell, Phys. Rev. D34, 822 (1986). 

[75] P. Lipari, M. Lusignoli and F. Sartogo, Phys. Rev. Lett. 74, 4384 (1995). 
[76] G.L. Fogli and E. Lisi, Phys. Rev. D52, 2775 (1995). 

[77] S. Fukuda et ai, (The Super-Kamiokande Collaboration), Phys. Rev. Lett. 85, 3999 
(2000). 

[78] T. Toshito, ptalk presented in Moriond 2001, France, 2001, hep-ph/0105023. 
[79] V. Zacek et al, Phys. Lett. B164, 193 (1985). 
[80] G.S. Vidyakin et ai, JETP Lett. 59, 390 (1994). 
[81] B. Achkar et ai, Nucl. Phys. B434, 503 (1995). 

[82] M. Appolonio et ai, Phys. Lett. B466, 415 (1999); Phys. Lett. B420, 397 (1998). 

[83] F. Boehm et al, Phys. Rev. D62, 072002 (2001); Phys. Rev. Lett. 84, 3764 (2000). 

[84] F. Dydak et al, Phys. Lett. B314, 281 (1984). 

[85] F. Bergsma et al, Phys. Lett. B142, 103 (1984). 

[86] I.E. Stockdale et ai, Z Phys. C27, 53 (1985). 

[87] L. Borodovsky et al, Phys. Rev. Lett. 68, 274 (1992). 

[88] N. Ushida et ai, Phys. Rev. Lett. 57, 2898 (1986). 

[89] C. Athanassopoulos et al, (The LSND Collaboration) Phys. Rev. Lett. 77, 3082 
(1996). 

[90] C. Athanassopoulos et al, (The LSND Collaboration) Phys. Rev. Lett. 81, 1774 
(1998). 

[91] A. Aguilar et al, (The LSND Collaboration) hep-ex/0104049. 

[92] Talk presented by The Karmen Collaboration (K Eitel for the Collaboration) in 
Neutrino 2000 held at Sudbury Canada, Nucl. Phys. Proc. Suppl. 91, 191 (2001). 

[93] S.H. Ahn et al, (The K2K collaboration), Phys. Lett. B511, 178 (2001). 



70 



CHAPTER 4 



The Energy Independent Solution to the Solar 

Neutrino Problem 



In the previous chapter we have seen that the suppression of the solar neutrino flux 
has been confirmed now by a number of experiments [|I], ^, 0, ^, covering a neutrino 
energy range of 0.2 — 20 MeV. However there is still considerable uncertainty regarding 
its energy dependence. An energy independent suppression factor has been advocated 
by several groups over the years This was shown to be disfavored however by the 
combined data on the total suppression rates [^]. More recently the weight of experimental 
evidence changed in favour of an energy independent solution following the SK data on 
the day/night spectrum, showing practically no energy dependence nor any perceptible 
day/night effect 

In this chapter we evaluate the experimental status of the energy independent solution 
vis a vis the MSW solutions after the inclusion of the SNO CC data. For the Ue survival 
probability Pu^ue we use the full expression given by eq. (|3.35| ) in chapter 3, including the 
coherent term containing the phase ^ which allows us to probe down the entire parameter 
space down to Am^ ~ 10~^^ eV^. For the composition of the neutrino fluxes seen in 
various experiments and the observed rates we refer to the Table p.2| presented in chapter 
3 where all the survival rates are shown relative to the standard solar model (SSM) 
prediction of BPBOO M. 



The apparent energy dependence in the survival rates of Table p.2| is conventionally 
explained in terms of the vacuum oscillation (VO), small and large angle MSW (SMA 
and LMA) as well as the LOW solutions The VO and SMA solutions show strong 
nonmonotonic energy dependence and SMA is essentially ruled out now with the inclusion 
of the SNO CC data and the SK day/night spectrum data. On the other hand the 
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spectrum data is compatible with the LMA and LOW solutions, which predict modest 
and monotonic decrease of the survival rate with energy. We will study the potential 
of an energy independent solution as a possible explanation of the global solar neutrino 
data. We shall see that with reasonable allowance for the renormalisations of the CI 
rate and the neutrino flux the data are described well by the energy independent 
solution. These renormalisations shall also be seen to improve the quality of fits with 
the oscillation solutions and enlarge the region of their validity in the parameter space 
substantially. Moreover we shall see that most of this enlarged region of parameter space 
shows weak matter effect on Pu^u^y implying practically energy independent suppression 
of the solar neutrino flux. Thus the energy independent solution can be looked upon as 
an effective parameterization of the oscillation solutions over this region. 

We shall perform a fit of the combined data on the survival rates of CI, Ga, SK and 
SNO (CC) and the SK spectrum at day and night with the energy independent solution as 
well as the four traditional solutions mentioned above. We have discussed the function, 
the error matrix and the method of analysis in chapter 3. In order to reconcile the energy 
independence of the spectrum with the apparent energy dependence in the rates of Table 
|3.2| , we shall consider the following variations in the CI rate and neutrino flux, since 
the CI experiment |0] has not been calibrated while the ^B neutrino flux is very sensitive 
to the underlying solar model. 

i) Since the CI experiment |@] has not been calibrated, there are several fits in the 
literature disregarding this rate |10|, |ll[] . We shall consider an upward renormali- 
sations of the CI rate by 20% (i.e. 2cr), which will push it marginally above the SK 
and SNO rates. This is favored not only by the energy independent solution but 
also the LMA and LOW solutions, showing monotonic energy dependence. 

ii) The ^B neutrino flux is very sensitive to the solar core temperature and hence to the 
underlying solar model. We shall consider a downward variation of the ^B neutrino 



flux of BPBOO 0, 




(4.1) 
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by upto 2(T. A downward renormalisation of this flux is favored by the energy 
independent solution and to a lesser extent by LOW, though not by the LMA 
solution. It is also favored by some helioseismic models, e.g. the model of [|1^ giving 
Jb = (4.16 ± 0.76) X 10^/cm^/sec. It may also be noted here that using la lower 



limits of the appropriate nuclear reaction rates Brun, Turck-Chieze and Morel [0 
have obtained a relatively low value of neutrino flux, Jb = 3.21 x 10^/cm^/ sec. 

In section 4.1 we present the results of the fit with the energy independent solu- 
tion. In section 4.2 we identify the regions of the parameter space with ~ 10% energy 
dependence and call these the quasi-energy independent regions. In section 4.3 we present 
the MSW solutions with renormalized CI and Xb (the flux normalization factor) and 
show the CL. contours in the Am^ — tan^^ plane. Most of the these allowed regions 
are seen to overlap with the quasi-energy independent zones. We end in section 4.4 with 
some discussions on the potential of future experimental measurements in distinguishing 
the energy independent solution with the MSW solutions. 

4.1 Energy Independent Solution 



Table ^]1| summarizes the results of fltting the global rate + spectrum data with an 
energy independent survival probability 

P.^.^ = l-\sm^2e. (4.2) 

We flt the data for CI observed and CI renormalized and for both the cases present the 
results for Xb flxed at SSM and Xb varying freely. It shows that even without any 
readjustment to the CI rate or the ^B neutrino flux the energy independent solution gives 
an acceptable goodness of flt (g.o.f.) of 24%. An upward renormalisation of the CI 
rate by 20% improves the g.o.f. to 42%. And varying the ^B neutrino flux downwards 
improves it further to 49%, which corresponds to a renormalisation factor Xb = 0.7 for 
the ^B neutrino flux. Note however that the g.o.f. and the best-flt value of the mixing 
angle for Xb = 1 are very close to those for Xb = 0.7. This is because the large error in 
the ^B neutrino flux of eq. (4J.) is already incorporated into the error matrix involved in 



the analysis for Xb = 1 case. 
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/ tan^ e \ 








Xb 


sin^ 26 \ or 


Xmin 


g.o.f 






\ cot a J 






Chlorine 


1.0 


0.93(0.57) 


46.06 


23.58% 


Observed 


0.72 


0.94(0.60) 


44.86 


27.54% 


Chlorine 


1.0 


0.87(0.47) 


41.19 


41.83% 


Renormalized 


0.70 


0.88(0.48) 


38.63 


48.66% 



Table 4.1: The best-fit value of the parameter, the Xmm ^'^d the g.o.f from a combined 
analysis of rate and spectrum with the energy independent solution given by eq. 



Traditionally the energy independent solution ( ^.2| ) is associated with the averaged 
vacuum oscillation probability at distances much larger than the oscillation wave-length, 
as originally suggested by Gribov and Pontecorvo As we shall see below however an 
effectively energy independent solution holds around the maximal mixing region over a 
wide range of Am? even after including all the matter effects in Sun and Earth. 

4.2 Regions of Energy Independent Solution in the Am^ — tan^ 9 
Plane 

The energy dependence of the survival probability arises from different sources in 
different regions of the parameter space. 

i) For Am? / E < 10"^'* eV the Earth regeneration effect can be safely neglected. Then 
the survival probability 

/ / Am^ L\ 

Pv.v. = Pi cos' 9 + P2 sin' 9 + 2^PiP2 sin ^ cos ^ I — ^ J , (4.3) 

where L is the distance between the neutrino production point at the solar core and 
its detection point on Earth; and P2 (= 1 — Pi) is the probability of the produced 
emerging from the Sun as the heavier mass eigen-state z/2 The coherent interference 
term, represented by the last term of eq. ( [4.3|) , is responsible for the nonmonotonic 
energy dependence of the VO and QVO solutions. 
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ii) For /S.iv} / E ~ 10~^^ — 10^^^ eV, the coherent term is neghgible but the Earth 
regeneration contribution can be significant leading to energy dependence of the 
survival probability. Besides over a large part of this region, represented by the 
MSW triangle, the z/g is adiabatically converted into z/2 in the Sun, i.e. P2 = 1- The 
the day/night averaged probability around these adiabatic zones is 

TjE sin^ 26 
4(1 -2riE cos2e + ril)' 

where 



P.,., = sin^ 9 + r , ^2^ ^ (4-4) 



Here p is the matter density in the EarthQ and Ye the average number of electrons 
per nucleon. The regeneration contribution is always positive and peaks around 
riE ~ 1, i.e. Am?' / E ~ 3 x 10"^'^ eV (cf. fig. p.4| ). Regions of the parameter space 
which give non-adiabatic neutrino propagation in the Sun pick up a strong energy 
dependence. 

iii) Finally for Am? / E > 10^^^ eV the survival probability can be approximated by 
the average vacuum oscillation probability of eq. ( [4.2|) . The MSW solutions (LMA 
and SMA) lie on the boundary of the regions ii) and iii), i.e. Am^ ~ 10^^ eV^ for 
E 1 MeV. The survival probability Pu^ue goes down from 1 — ^ sin^ 29 (> 0.5) to 
sin^ (< 0.5) in going up from Ga to SK (SNO) energy. 

All the solar neutrino rates except that of SK have ~ 10% error, which is also true for 
the SK energy spectrum. The SK normalization has at least similar uncertainty from the 

neutrino fiux. Therefore we shall treat solutions, which predict survival probability 
PueUe within 10% of eq. ( [4.2| ) over the Ga to SK energy range, as effectively energy 



independent solutions. Moreover the predicted Ga, CI and SNO rates will be averaged 
over the respective energy spectra, while the predicted SK rates will be averaged over 
0.5 MeV bins, corresponding to the SK spectral data, since experimental information is 
available for these averaged quantities only. 



^^Though in eqs. (4.4) and (4.5) we have shown the regeneration for a constant density Earth, we have 



used the two slab model for actual calculations. 
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Figure 4.1: The quasi-energy independent regions in Am^ — tan^ 9 parameter space where 
the solar neutrino survival probability agrees with eq. ( [4 .21) to within 10% over the range 
of Ga to SK energies. 



The shaded parts of fig. O shows the regions of effective energy independent solution 
as per the above definition. The parameter space has been restricted to Am^ < 10^^ eV^ 
in view of the constraint from the CHOOZ experiment |T3. One sees that the energy 
independent solution ([4.2|) is effectively valid over the two quasi- vacuum oscillation regions 
lying above and below the MSW range. Moreover it is valid over a much larger range 
of Am^ around the maximal mixing region, since the solar matter effect does not affect 
Pu^u^ at tan^ 9 = 1. It is this near- maximal mixing strip that is relevant for the observed 
survival probability, Pv^v^ ~ 1/2. The upper strip (Am^ = 10^^ — 10^^ eV^) spans the 
regions iii) and part of ii) till it is cut off by the Earth regeneration effect. The lower 
strip (Am^ = 10^'' — 5 x 10^^° eV^) spans parts of region ii) and i) till it is cut off by 
the coherent term contribution. Note that this near-maximal mixing strip represents a 
very important region of the parameter space, which is favored by the socalled bimaximal 
mixing models of solar and atmospheric neutrino oscillations |]T6|, p!7| . 
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One can easily check that averaging over the SK energy bins of 0.5 MeV has the effect 
of washing out the coherent term contribution for Am^ > 2 x 10~^ eV^. But including 
this term enables us to trace the contour down to its lower limit. It was claimed in ref. 



T8|| that the lower strip disappears when one includes the coherent term contribution. 
This may be due to the fact that their predicted rate in the SK energy range was not 
integrated over the corresponding bin widths of 0.5 MeV. 

To get further insight into the oscillation phenomenon in the maximal mixing region 
we have plotted in fig. [4.2| the predicted survival rates at maximal mixing against Am^ 
for the Ga, CI, SK and SNO experiments. In each case the rate has been averaged over 
the corresponding energy spectrum. This is similar to the fig. 7 of Gonzalez- Garcia, 



Pena-Garay, Nir and Smirnov ||T0|. As in ||T0| the predictions have been shown relative 



to the central value of the various neutrino fluxes of BPBOO along with those differing 
by ±1(7 from the central value, e.g. for the flux Xb = 1 ± 0.2. We have found that 
these curves are in good agreement with the corresponding ones of [|10[- The two regions 



of < 10% energy dependence for maximal mixing are indicated by vertical lines. One 
can easily check that in these regions the central curves lie within 10% of the energy 
independent prediction R = 0.5 (note that the SK rate apparently looks higher due to 
the neutral current contribution). One can clearly see the energy dependence due to the 
violent oscillations in the VO region (on the left) and the Earth regeneration peak in the 
LOW region (in the middle), particularly for the Ga experiment. It should be noted that 
the gap between the two energy independent regions due to the Earth regeneration effect 



in fig. |4.1| is a little narrower than here. This is due to a cancellation between the positive 
contribution from the Earth regeneration effect and the negative contribution from the 
solar matter effect at tan^ ^ < 1. It ensures that the resulting survival rate agrees with 
the energy independent solution ( [4.2|) over a somewhat larger range of Am^ slightly below 
the maximal mixing angle. 



The observed rates from the Ga, CI, SK and SNO experiments are shown in fig. 



as horizontal lines along with their la errors. With 20% downward renormalisation of 
the neutrino flux {Xb = 0.8) the energy independent prediction is seen to agree with 
the SK rate and also approximately with SNO. It is higher than the CI rate by about 
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Figure 4.2: The predicted survival rates at maximal mixing against Am? for CI, Ga, SK 
and SNO experiments (See text for details.) 
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Nature of 
Solution 


in eV^ 




■ 


g.o.f 




SMA 


5.28 X 10" 


-6 


3.75 X 10"^ 


51.14 


9.22% 


CI 


LMA 


4.70 X 10- 


-5 


0.38 


33.42 


72.18% 


Obsvd. 


LOW 


1.76 X 10- 


-7 


0.67 


39.00 


46.99% 




VO 


4.64 X 10- 


-10 


0.57 


38.28 


50.25% 




SMA 


4.94 X 10- 


-6 


2.33 X 10-^ 


50.94 


9.54% 


CI 


LMA 


4.70 X 10- 


-5 


0.38 


30.59 


82.99% 


Renorm. 


LOW 


1.99 X 10- 


-7 


0.77 


34.26 


68.57% 




VO 


4.61 X 10- 


-10 


0.59 


32.14 


77.36% 



Table 4.2: The best-fit values of the parameters, the Xmin ^"^^ g-O-f from a combined 
analysis of the CI, Ga, SK and SNO CC rates and the SK day-night spectrum in terms 
of Ve oscillation into an active neutrino, including the matter effects. Xb is kept fixed at 
the SSM value (=1.0). 



2a. The agreement with the Ga rate can be improved by going to a little smaller 6 and 
compensating the resulting deviation from the other rates by a somewhat smaller Xb as 



in Table [4.1| . Nonetheless the maximal mixing solution for Xb = 0.8, shown here, is in 
reasonable agreement with the observed rates over the energy independent regions. The 
Earth regeneration effect can be seen to improve the agreement with the Ga experiment 
for the LOW solution. 



4.3 The SMA, LMA, LOW and VO Solutions 



Tables [4.2| and |4.3| summarise the results of fits to the global rates + spectrum data 



in terms of the conventional oscillation solutions. Table ^]2| shows solutions to the data 
with observed and renormalized CI rate with the neutrino flux of BPBOO {Xb = 1), while 
Table [4.3| shows the effects of renormalizing this neutrino flux downwards by 25% 
{Xb = 0.75). The corresponding 90%, 95%, 99% and 99.73% (3cr) contours are shown in 

fig. M 



As we see from these tables and fig. the SMA solution gives rather poor fit in 



each case, with no allowed contour at 3a level. This result agrees with the recent fits of 
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Nature of 
Solution 


in eV^ 


tan2 6' 


• 


g.o.f 




SMA 


5.28 X 10- 


-6 


3.75 X 10"^ 


48.39 


14.40% 


CI 


LMA 


4.65 X 10- 


-5 


0.49 


38.90 


47.44% 


Obsvd. 


LOW 


1.74 X 10- 


-7 


0.71 


39.91 


42.95% 




VO 


4.55 X 10- 


10 


0.44 


37.17 


55.36% 




SMA 


8.49 X 10" 


-6 


1.78 X 10-^ 


50.77 


9.82% 


CI 


LMA 


4.64 X 10- 


-5 


0.51 


34.48 


67.61% 


enorm. 


LOW 


2.09 X 10- 


-7 


0.81 


33.47 


71.97% 




VO 


4.59 X 10- 


10 


0.53 


30.63 


82.86% 



Table 4.3: Best fits to the combined rates and SK day-night spectrum data in terms of 
Ue oscillation into active neutrino with a fixed Xb = 0.75. 



1 19(1 to the global data including SNO. The fit of to these data shows a small allowed 
region for SMA solution at 3a level due to a slightly different method of treating the 
normalization in the SK spectrum data, as explained there. 

The LMA and LOW solutions give good fits to the original data set, which improve 
further with the upward renormalisation of the CI rate by 2a. This is because the mono- 
tonic decrease of rate with energy, implied by these solutions, favors the CI rate to be 
marginally higher than the SNO and SK rates as mentioned earlier. For the renormal- 
ized CI case, downward renormalisation of the neutrino fiux by 25% is seen to give 
a modest increase (decrease) of g.o.f. for the best LOW (LMA) solution. On the other 
hand the allowed ranges increase in both cases as we see from fig. ^4.3| b and c. Combined 
together they imply that the downward renormalisation of the neutrino fiux modestly 
favors the LOW solution but makes little difference to the LMA. Its main effect on these 
two solutions is increasing their allowed ranges in the parameter space. Comparing fig. 
|4.3| c with fig. shows that much of these enlarged ranges of validity correspond to 
effectively energy independent solution. It is interesting to note that the best-fit values of 
parameters in the LMA region are same for CI observed and CI renormalized cases while 
the xiiiri improves for the latter. This shows that the best-fit already chose a probability 
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Figure 4.3: The 90%, 95%, 99% and 99.73% C.L. allowed area from the global analysis of 
the total rates from CI (observed and 20% renormalized) , Ga, SK and SNO (CC) detectors 
and the 1258 days SK recoil electron spectrum at day and night, assuming conversions to 
active neutrinos. 



at CI energy, which is a little higher than that at SK/SNO energy. Renormalizing the CI 
rate brought that point up to the fitted curve without changing the best-fit parameters. 

While the best vacuum solution seems to show remarkably high g.o.f. particularly for 
renormalized CI rate and neutrino flux, its regions of validity are two miniscule islands 
just below the lower energy independent strip of fig. [4.3| b,c. This solution has also been 
obtained in the global fits of ref. [|r^, as well as the SK fit to their rate and spectrum 
data [^. The position and size of its range of validity suggest this to be a downward 
fluctuation of the effectively energy independent quasi vacuum oscillation rather than a 
genuine VO solution of Just-So type. To get further insight into this solution we have 
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Nature of 
Solution 


Xb 


in eV^ 


tan2 6' 


• 


g.o.f 




SMA 


0.51 


5.25 X 10" 




3.44 X 10-4 


46.83 


15.41% 


CI 


LMA 


1.18 


4.73 X 10- 


-5 


0.33 


32.32 


72.89% 


Obsvd. 


LOW 


0.88 


1.75 X 10- 


-7 


0.67 


38.75 


43.57% 




VO 


0.70 


4.55 X 10" 


10 


0.44 


37.24 


50.44% 




SMA 


0.48 


4.66 X 10- 


-6 


2.32 X 10-4 


46.18 


17.01% 


CI 


LMA 


1.15 


4.71 X 10- 


-5 


0.36 


30.32 


80.80% 


Renorm. 


LOW 


0.83 


2.03 X 10- 


-7 


0.79 


33.18 


69.17% 




VO 


0.75 


4.63 X 10- 


10 


0.55 


30.56 


79.92% 



Table 4.4: Best fits to the combined rates and SK day-night spectrum data in terms of 
Ve oscillation into active neutrino with Xb free. 



analysed the resulting energy dependence. It shows practically no energy dependence 
below 5 MeV, but a 15% fall over the 5 — 12 MeV range. The latter seems to follow the 
SK spectral points rather closely amidst large fiuctuation. To check the stability of this 
trend we have repeated the fit to the SK spectral data points, plotted over 8 broad energy 
bins shown in [Q, which show much less fiuctuation than the 2 x 19 points sample. The 
solution completely disappears from this fit. This confirms that the above VO solution is 
simply an artifact of the sampling of the SK spectral data. 

For completeness we summarise in Table [4.4| the best fits of the above solutions with 
free neutrino fiux normalization. The SMA solution favors a very low neutrino 
fiux {Xb — 0.5), which raises the SK and SNO rates more than the CI, thus accentuating 



the nonmonotonic energy dependence of Table |3.2| . Still the g.o.f. of the SMA solution is 
rather marginal. On the other hand the LMA solution favors Xb > 1, which suppresses 
the SK and SNO rates more than the CI, resulting in a monotonic decrease of rate with 



energy. But the corresponding g.o.f. are no better than those of Table 4.2. The results 



of the LOW and VO fits are similar to those of Table 4.3 
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4.4 Comparisons and Discussions 



Let us conclude by briefly discussing whether some of the forthcoming neutrino ex- 
periments will be able to discriminate between the energy independent and the MSW 
solutions. In particular the SNO experiment is expected to provide both the charged 



current and neutral current scattering rates over roughly the same energy range as SK. 
Thus the neutrino flux can be factored out from their ratio, CC/NC. For oscillations 
into active neutrino the corresponding double ratio Rsno/^sno predicted to be larger 
than 0.5 for the energy independent solution and smaller than 0.5 for the LMA and LOW 
solutions. In the absence of the neutral current data from SNO one can try to make a 
similar comparison with the ratio of SK elastic and SNO charged current scattering rates, 

R'Ij, = XbPu^u. + r(l - P.^.JXb, r = a^Ja,, 0.17, (4.6) 



R'sNO - XsPu^Ue, (4.7) 

where we have assumed a common survival rate neglecting the small difference between 
the SK and SNO energy spectra One can eliminate Py^v^ from the two rates; and the 
resulting neutrino flux can be seen to be in good agreement with the BPBOO estimate 
1^. Alternatively one can factor out the flux from the ratio 

Table shows the best fit values of the above ratio for the LMA, LOW and the energy 
independent solutions along with the corresponding predictions for R'sno/ ^sno — Pi^ei^e- 
The predictions of the maximal mixing solution is also shown for comparison. While the 
LMA and LOW predictions for the CC/NC ratio differ by ~ 50% they differ by only 
about ~ 15% in the case of the ES/CC ratio. The observed ratio R^k / P'sno seen to 
favour the LMA over the LOW and energy independent solutions; but even the largest 
discrepancy is only ~ 1.5cr. With the expected sample of several thousand CC and NC 
events from SNO one expects to reduce the la error for each of these ratios to about 
5%. Then one will be able to discriminate between the three solutions meaningfully, 
particularly with the help of the CC/NC ratio from SNO. On the other hand the LOW 
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tan2 6' 


pC'C / jdNC 
^SNOl ^SNO 


1 pC'C 
^SKI ^SNO 


Expt. value of 

T-> RS 1 j-)CC 
^SKI ^SNO 


LMA 


4.7 X 10-^ 


0.38 


0.30 


1.36 




LOW 


1.99 X 10-^ 


0.77 


0.45 


1.19 




energy-independent 




0.47 


0.56 


1.13 


1.33 ± 0.13 


maximal mixing 




1.0 


0.5 


1.15 





Table 4.5: The values of the ratios Rsk / Rsno ^sno/ ^sno best-fit values for 

the LMA, LOW and energy independent solutions for the renormalized CI and Xb = 1.0 
case. Also shown are the predictions for the maximal mixing (Pee = 0.5) solution and the 
experimental value of Rsk/^sno- 



solution predicts a large Day-Night asymmetry of > 10% for the ^Be neutrino [^0], |Tl|] at 
the Borexino and the KamLAND experiments. This will be able to distinguish 



the LOW from the LMA and the energy independent solutions. Lastly it should be noted 
that the reactor neutrino data at KamLAND is expected to show oscillatory behavior for 



the LMA solution ||2^ , which will help to distinguish it from the LOW or a generic energy 
independent solution. 

In summary the recent SK data on day/night spectrum is in potential conflict with 
the apparent energy dependence in the suppression rates observed in Ga, CI, SK and 
SNO experiments. Including matter effects one can get acceptable oscillation solutions to 
both rates and spectrum data only over limited regions of mass and mixing parameters. 
However an upward renormalisation of the CI rate by 20% {2a) results in substantial im- 
provement of the quality of fit. Moreover a downward renormalisation of the neutrino 
flux by 25% (1.5(7) as suggested by the helioseismic model enlarges the allowed region 
of the parameter space substantially. Over most of this enlarged region the energy de- 
pendence resulting from the matter effects is too weak to be discernible at the present 
level of experimental accuracy. Hence with these renormalisations of the CI rate and the 
^B neutrino flux the data can be described very well by an energy independent solution. 
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CHAPTER 5 



A Three Generation Oscillation Analysis of the 
Super-Kamiokande Atmospheric Neutrino Data 
Beyond One Mass Scale Dominance Approximation 



In chapter 3 we have made detailed analysis of the solar and atmospheric neutrino 
data and presented our results in the two generation framework. In real world there are 
three active flavors of neutrinos [jl| . If one wants to do an analysis of the SK atmospheric 
neutrino data [0 in a three-generation framework, then there are two possibilities. The 
most popular three-generation picture in the context of the SK data is the one where one 
of the mass squared differences is in the solar neutrino range and the other is suitable 
for atmospheric neutrino oscillations @, ^ ||, In such a scheme one mass scale 

dominance (OMSD) applies for atmospheric neutrinos and the relevant probabilities are 
functions of two of the mixing angles and one mass squared difference. This picture 
however cannot explain the LSND results 0. In this chapter we perform a three flavor 
X^-analysis of the SK atmospheric neutrino data assuming a mass pattern with A12 — A13 
fixed in the eV^ range (henceforth Amfj = Aij) and allowing the other mass scale to vary 



arbitrarily. We also incorporate the CHOOZ reactor results [jTO[ into the analysis. Apart 
from being suitable to explain the SK atmospheric neutrino data this spectrum is also 
interesting for the laboratory based neutrino oscillation experiments as the higher mass 
scale is explorable in the short base line experiments, whereas the lower mass scale can 
be probed in the long base line experiments. In this scheme to a good approximation, 
neutrino oscillation in the short-base line accelerators and reactors will be governed by 



one (the higher) mass scale [TTI, |T2| - and only two of the mixing angles appear in the 
expressions for the oscillation probabilities. For the atmospheric and the long baseline 
experiments the characteristic energy and length scales are such that in general both 
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mass differences are of relevance and the probabilities involve all the three mixing angles. 
However the higher mass scale gives rise to Am^ independent average oscillations and it 
does not enter the fit directly. We determine the best-fit values of A23 and the three 
mixing angles by performing a analysis of the 

• SK atmospheric neutrino data 

• SK atmospheric and CHOOZ data 

We do a bin by bin analysis of the data taking into account the errors and their correla- 
tions and use the function described in eq. ( p.62| ) of chapter 3. Finally we compare 
the allowed values of the mixing angles as obtained from the above analysis with those 
allowed by the other accelerator and reactor neutrino oscillation data including LSND 
and KARMEN2. 

The mass scheme of this paper was first considered in |T3|, |14| after the declaration 



of the LSND result. These papers performed a combined three generation analysis of 
accelerator and reactor results as well as the Kamiokande atmospheric neutrino data. 
Three-generation picture with the higher mass difference in the eV^ range and the lower 
mass difference in the atmospheric range has also been considered in H^, ^ (pre-SK) and 



T7| , 0, |T9|, (post-SK). These papers attempted to explain both solar and atmospheric 
neutrino anomalies mainly by maximal z/^ <-> Ve oscillations driven by /S^atm ~ 10^'^ eV^. 



Although it was claimed in [|T7|, [T^ that this scenario can provide a good fit to all the 



available data on neutrino oscillations, it was shown in ||T9[ and also later in ||20| that this 
scenario cannot reproduce the zenith angle dependence of the SK atmospheric neutrino 
data. 

In this chapter our aim is to determine the allowed oscillation parameter ranges consis- 
tent with SK atmospheric, CHOOZ, LSND and other accelerator and reactor experiments. 
The solar neutrino problem can be explained by invoking a sterile neutrino. We discuss in 
the conclusions how the solar neutrino flux suppression can be explained in our scenario. 

The plan of the chapter is as follows. In section 5.1.1 we present the formalism for 
three-generation oscillation analysis in vacuum. We calculate the required probabilities 
including the earth matter effects in section 5.1.2. We use this expression of the probability 
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for the actual calculation of the number of events. In section 5.2 we present the three- 
generation analysis of only SK atmospheric neutrino data using the computer code 
described in detail in chapter 3. In section 5.3 we present the combined analysis of SK 
and CHOOZ data. In section 5.4 we compare the allowed values of mixing angles from 
the above analyses with those allowed by the other accelerator and reactor data including 
the latest results from LSND and KARMEN2. In section 5.5 we discuss the implications 
of our results for the future long baseline experiments and end in section 5.6 with some 
discussions and conclusions. 

5.1 Three-Flavor Analysis 

5.1.1 The Vacuum Oscillation Probabilities 

The general expression for the probability that an initial Va of energy E gets converted 
to a after traveling a distance L in vacuum is given by eq. ( |2.6| ) in chapter 2. The actual 
forms of the various survival and transition probabilities depend on the neutrino mass 
spectrum assumed and the choice of the mixing matrix U relating the flavor eigenstates 
to the mass eigenstates. We choose the flavor states a = 1,2, and 3 to correspond to e, 
yU and T respectively. The most suitable parameterization of U for the mass spectrum 
chosen by us is ?7 = R13R12R23 where Rij denotes the rotation matrix in the ij-plane. 
This yields: 



where Cij = cos 6ij and Sij = sin^ here and everywhere else in this chapter. We have 
assumed CP-invariance so that U is real. The above choice of U has the advantage that 
$23 does not appear in the expressions for the probabilities for the laboratory experiments 



Pf.i', = 1 - 2Ci3C^2 + 2C13C12 - 4(ci3Si2C23 - Si3S23)^(ci3Si2S23 + ^13023)^ S23 (5.2) 



u 




(5.1) 




12C23S23(C13S12C23 - S13S23) (C13S12S23 + S13C23) S23 (5.3) 
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1 _ e2 _ 4 A „2 2 c 

L ^'-12''12 ^'^12'^23''23 ^23 



(5.4) 



where S23 = sin^(7rL/A23), where Ajj is given by eq. (pTTp. Apart from the most general 
three generation regime, the following limits are of interest, as we will see later in the 
context of the SK data: 



1. The two-generation limits 

Because of the presence of more parameters as compared to the one mass scale 



dominance picture there are twelve possible two-generation limits with the os- 
cillations driven by either A^snd or Aatm- Below we list these hmits specifying 
the mass scales that drive the oscillations: 



S12 

Sl3 
Sl3 

S12 
S12 



2. = 0.0 



0, Si3 

1, -Sl3 

0, S23 

0, S23 
0, S23 

0,S23 








1 


1 



A, 

^ei Alsnd), 
^e, Aisnd), 

^T, Alsnd), 
lsnd), 



S12 - 
S12 - 

■S13 
■S13 

S12 - 
S12 - 



In this limit the relevant probabilities become 



0, Si3 ^1 (z/^ - z/e, Aatm 

1, si3 -^1 (z/e - z/^, Aatm] 

^ 1,S23 (z/^ 



'^T, A lsnd) 



^ 1, S23 ^1 (i^f,- i^T, Alsnd) 
1, -523^0 (z/e-z/^, Aatm) 



1,S23^1 {l^e-l^T, A 



ATM 



1 - 2c^3Si3 + 4:sl^cl^sl.^S23 



(5.5) 



4'5i3S23C235'23 



(5.6) 



1 ~ 4C23S235'23 



(5.7) 



Thus Puf^Uf, is the same as the two generation limit, Pui^ue is governed by two of the 
mixing angles and one mass scale and P^^^e is governed by two mixing angles and 
both mass scales. 
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For this case the probabihties take the form 



Pv^ve = 1 ~ 20^2^12 — 45^20235235*23 (5.8) 

PueU^ = 2C]^2'^12 ~ 4C]^2'^12'^23'^23'^23 (5-9) 
Puf^Uf, = 1 - 2c^2Si2 - 4C^2C23"523'S'23 (5.10) 

In this case the probabihties are governed by two mass scales and two mixing angles. 

We note that for cases (2) and (3) the probabilities are symmetric under the transforma- 
tion ^23 7r/2 — ^23- The probabilities for these cases are functions of at most two mixing 
angles as in the OMSD case but they are governed by both mass scales making these 
limits different from the OMSD limit. 

5.1.2 Earth Matter Effects 

Since on their way to the detector the upward going neutrinos pass through the earth, 
it is important in general to include the matter effect in the atmospheric neutrino analysis. 
The matter contribution to the effective squared mass of the electron neutrinos: 

A = 2V2GFEne (5.11) 

where E is the neutrino energy and rif, is the ambient electron density. Assuming a typical 
density of 5 gm/cc and E = 10 GeV, the matter potential A ~ 3.65 x 10^'^ eV^ and since 
this is of the same order as A23 in our case, matter effects should be studied carefully. 
The mass matrix in the flavor basis in presence of matter is given by 

M| = f/ f/t + Ma (5.12) 

where is the mass matrix in the mass eigenbasis, U is the mixing matrix and 

(A 0\ 

Ma = \ (5.13) 
VO 0/ 
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Since A12 ~ A13 ^ A23 ~ A, one can solve the eigenvalue problem using the degenerate 
perturbation theory, where the A23 and A terms are treated as a perturbation to the 
dominant A12 and A13 dependent terms. The mixing angle in matter is then given by 



tan 2^2! = r Tn .7 2 2 2% (5-14) 

A23 cos 2^23 -^(sf3-cf3S?2) 

while the mass squared difference in matter turns out to be 



AM 

^23 



(A23 COS 2^23 - A{sl, - cl,sl^)Y + (A23 sin 2^23 - Asi2 sin 2^13)'! (5-15) 



The mixing angles 612 and ^13 as well as the larger mass squared difference A12 remain 
unaltered in matter. From eq. ( |5.14|) and (|5.15| ) we note the following 



In the limit of both S12 —>■ and S13 0, the matter effect vanishes and we recover 
the two- generation z/^ — z/^ limit. 

The resonance condition now becomes A23 cos 26^23 = A^sl^, — Ci3sf2)- So that for 
A23 > 0, one can have resonance for both neutrinos - if 5^3 > 0^35^2 ^ well as for 
antineutrinos - if < 0^33^2- This is different from the OMSD picture where for 
Am^ > only neutrinos can resonate. 



In the limit of S12 —>■ 



. ooM A23 sin 2^23 

tan 26^23 = — (5.16) 

2' A23COs2^^23-^sf3 ^ ^ 



Here one gets resonance for neutrinos only (if A23 > 0). 
In the limit S13 — » 



tan2e^^ = (5.17) 

A23 cos 2023 + Asj2 ^ ' 

For this case for A23 > 0, there is no resonance for neutrinos but antineutrinos can 

resonate. 

In the limit where A23 — 

tan2g,^3^= y^^'f? , A^i = A{4, + cl,sl2) (5-18) 

■5l3 C]_3S]_2 
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Thus even for small values of A23 < 10"'* the mass squared difference in matter 
is ~ y4 and one may still hope to see oscillations for the upward neutrinos due 
to matter effects. The other point to note is that the mixing angle in matter 6^ 
depends only on 6*12 and 613 and is independent of the vacuum mixing angle 623 
and A23. Contrast this with the OMSD and the two-generation z/^ — z/g oscillations. 
For both the two-generation — z/g as well as the three-generation OMSD case, for 
Am^ 0, the mixing angle in matter 0, but for the mass spectrum considered 
in this paper the tan 2^23 maybe large depending on the values of s^2 ^13- Hence 
we see that the demixing effect which gives the lower bound on allowed values of 
Am^ in the two generation — z/g or the three-generation OMSD case, does not 
arise here and we hope to get allowed regions even for very low values of A23. On 
the other hand even small values of 623 in vacuum can get enhanced in matter. This 



special case where A23 ~ was considered in an earlier paper |22 



In the limit of S23 



.M -Asi2 sin 2013 



tan 2^2! = Tr-2 5^ (5-19) 



While for S23 — > 1 



, oaM -v4si2 sin 2^13 

tan 2^23 = 77^2 ^T-TT (5-20) 

-^23 — A[Si3 ~ C13S12J 



For the last two cases, corresponding to sin^ 26*23 ~^ 0, again the mixing angle 
^23 in matter is independent of its corresponding value in vacuum and hence for 
appropriate choices of the other three parameters, A23, 5^2 and 5^3, one can get 
large values for sin^ 2^^| even though the vacuum mixing angle is zero. 

The amplitude that an initial Ua of energy E is detected as z/^j after traveling through the 
earth is 

<j,X,p i,j,k,l 
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where we have considered the earth to be made of two slabs, a mantle and a core with 
constant densities of 4.5 gm/cc and 11.5 gm/cc respectively and include the non-adiabatic 
effects at the boundaries. The mixing matrix in the mantle and the core are given by 
f/*^™ and U'^^" respectively. Ef^ ^ mfx/'^E, X = core(mantle) and rriix is the mass of 
the i^^ neutrino state in the core(mantle). The neutrino is produced at time to, hits the 
earth mantle at ti, hits the core at t2, leaves the core at and finally hits the detector 
at time t. The Greek indices {a, A, p) denote the flavor eigenstates while the Latin indices 
{i,j,k,l) give the mass eigenstates. The corresponding expression for the probability is 
given by 

P{u^,to,up,t) = \A{u^,to,up,t)\^ (5.22) 

For our calculations of the number of events we have used the full expression given by 
eq. dOTD and ( ^1221) . 



5.2 x^"^^^lysis of the SK Data 



We do a bin by bin x analysis of the 1144 days SK atmospheric neutrino data [E^ 



following the procedure of statistical analysis given in We minimize the function 

defined in chapter 3, eq. (|3.62|) . For contained events there are forty experimental data 
points. The probabilities for the atmospheric neutrinos are explicit functions of one mass- 
squared difference and three mixing angles making the number of degrees of freedom 
(d.o.f) 36. The other mass squared difference gives rise to Am^ independent average 
oscillations and hence does not enter the fit as an independent parameter. 
For the three-generation scheme considered here the Xmin ^^e best-fit values of pa- 
rameters that we get are 

• xLn/d-o.f. = 34.65/36, A23 = 0.0027 eV^, s^g = 0.51, s?2 = 0.04 and sfg = 0.06 

This solution is allowed at 53.28% probability. 

The solid(dashed) lines in fig. ^T!] present the variation of the Ax^ = ~ Xmin the 
SK data, with respect to one of the parameters keeping the other three unconstrained, 
when we include (exclude) the matter effect. In fig. |5.1| (a) as we go towards smaller values 
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Figure 5.1: The variation of Ax^ = X^~Xmm with one of the parameters keeping the other 
three unconstrained. The sohd (dashed) hne corresponds to only SK data when matter 
effects are included (excluded) while the dotted curve gives the same for SK+CHOOZ. 
The dashed-dotted line shows the 99% C.L. limit for 4 parameters. 



of A23 around 10~^ eV^ the effect of matter starts becoming important as the matter term 
is now comparable to the mass term. If matter effects are not there then for values of A23 
~ 10~^ eV^ the 5*23 term in eq. (|5.4|) is very small and there is no up-down asymmetry 
resulting in very high values of is evident from the dashed curve. If the matter 

effects are included, then in the limit of very low A23 the matter term dominates and A2I 
is given by eq. ( ^.18| ). Since this term ~ 10~^ eV^ there can be depletion of the neutrinos 
passing through the earth causing an updown asymmetry. For A23 around 10~^ eV^, there 
is cancellation between the two comparable terms in the numerator of eq. ( ^.14| ) and the 
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mixing angle becomes very small and hence the around these values of A23 comes out 
to be very high. 

Fig |5.1j (b) illustrates the corresponding variation of Ax^ with S23 while the other three 
parameters are allowed to vary arbitrarily. For small and large values of S23 the inclusion 
of matter effect makes a difference. For S23 either very small or large (sin^ 26'23 0) 
the overall suppression of the z/^ flux is less than that required by the data if vacuum 
oscillation is operative and so it is ruled out. If we include matter effects then in the limit 
of S23 = and S23 = 1 the matter mixing angle is given by eqs. (|5.19| ) and (|5.2CI|) , which 
can be large for suitable values of sfg and s^2 ^^"^ hence one gets lower even for these 
values of S23. 

In figs |5]l|(c) and pJ](d) we show the effect of s^2 s^g respectively on A^^- From 
the solid and the dashed lines it is clear that matter effects do not vary much the allowed 
ranges of s^2 s^g. 

The dashed-dotted line in the figure shows the 99% CL. (= 13.28 for 4 parameters) 



limit. In Table |5]^ we give the allowed ranges of the mixing parameters, inferred from fig. 
O at 99% CL. for the 1144 day SK atmospheric data, with and without matter effects. 





A23 in eV2 


s2 
23 


■'12 


■^13 


with 
matter effects 


1.6xlO~^<A23< 7.0x10^3 
A23 < 6.5 X 10-5 


0.26 < < 0.77 
> 0.85 


s?2 < 0.21 


s?3 < 0.55 


without 
matter effects 


5xlO-^<A23< 7.0x10-3 


0.27 < < 0.74 


s?2 < 0.21 


s?3 < 0.6 



Table 5.1: The allowed ranges of parameters for the 1144 day SK data. 



5.2.1 Zenith-Angle Distribution 

Since the probabilities in our case are in general governed by two mass scales and 
all three mixing angles it is difficult to understand the allowed regions. To facilitate the 



qualitative understanding we present in figs. |5.2| and |5.3| the histograms which describe the 
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Figure 5.2: (a)The zenith angle distribution of the lepton events with A23 = 0.002 eV^ 
and S23 = 0.5 for various combinations of and 5^3. is the number of events as 
given by eq. ( p.58| ) and No is the corresponding number with survival probability 1. The 
panels labelled SGa and MGa (a can be e or /x) give the histograms for the sub-GeV and 
multi-GeV a-events respectively. Also shown are the SK experimental data points with 
± la error bars. (b)Same as in (a) for fixed s^2 = 0-1 s^g = 0.0 varying A23 and S23. 
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(d) 



0.8 



0.6 



0.4 



SGe 



MGe 



MG^ 



-1 





-1 


-1 





COS0 


COS0 


COS0 




COS0 


line type 




=23 


=2 


=2 

^13 


long dashed 


0.0027 eV^ 


0.5 


0.0 


0.0 


short dashed 


0.0027 eV^ 


0.51 


0.04 


0.06 


dotted 


0.00001 eV^ 


1.0 


0.2 


0.4 



'13 



0.1 for different A23 and s 



2 

23 



Figure 5.3: (c)Same as in (a) fixing — 0-0 ^ 
values. (d)Tfie long-dashed (short-dashed) line gives the zenith angle distribution of the 
lepton events for the best-fit cases of the two-generation (three-generation) oscillation 
solutions for SK. The dotted line gives the corresponding distribution for A23 = 10~^ 
eV^ = 0.2, s?, = 0.4 and s^, = 1.0. 



'13 



'23 
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zenith angle distribution. The event distributions in these histograms are approximately 
given by, 



J" ^Pu^u^+ j."Pu,u^ (5.23) 
p.... + ^P.... (5.24) 



where the quantities with suffix indicates the no-oscillation values. For the sub-GeV 
data N/j^^/NeQ ~ 2 to a good approximation however for the multi-GeV data this varies in 
the range 2 (for cosO =0) to 3 (for cos 9 = ±1) 0. 

In fig. |5]^(a) we study the effect of varying ^13 fixed values of A23 = 

0.002 eV^ and S23 = 0.5. From eq. ( ^.14| ), ( p.l5| ) and from fig. we see that for the 
values of the A23 and S23 considered in this figure the matter effects are small and we 
can understand the histograms from the vacuum oscillation probabilities. The thick solid 
line shows the event distribution for 5^2 — ■^13 — O-^- 5^3 increases from 0, 
keeping s^2 0, from eqs. ( |5.5| ) and ( |5.6| ) Pu^u^ decreases from 1 and Pu^u^ increases from 
zero resulting in a net electron depletion according to eq. ( p.24| ). The long dashed line 
corresponds to 5^3 = 0.3 for which the electron depletion is too high as compared to data. 
The muon events are also affected as Pui^ue increases with increasing 5^3 even though Pu^^u^ 
is independent of sf^. On the other hand for 5^3 = 0.0, the effect of increasing is to 
increase the number of electron events and decrease the number of muon events according 
to eqs. (|578|) , ( |5.9|) , (|5.10| ), (|5.24|) and ( |5.23|) . This is shown by the short-dashed and 
dotted lines in fig. pl2|(a). For 5^2 = 0-2 the electron excess and muon depletion both 
becomes too high as compared to the data. For the case when both 5^2 ^13 
0.1 the electron depletion caused by increasing 5^2 ^^"^ the excess caused by increasing 
s^3 gets balanced and the event distributions are reproduced quite well, shown by the 
dashed-dotted line. 

In fig. p.2| (b) we study the effect of varying S23 and A23 in the limit of — 
with s^3 fixed at 0.1. Although we use the full probabilities including the matter effect, 
for 0.004 eV^ this is not so important and one can understand the histograms from the 
vacuum oscillation probabilities. For fixed A23 as S23 increases, Pui^u^ decreases, making 
the muon depletion higher. This is shown in the figure for two representative values of 
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A23- The electron events are not affected much by change of S23. The shght increase 
with is due to increase of both P^^y^ and Pu^^ue- To understand the dependence on 
A23 we note that for S23 = 0.2, if one looks at the vacuum oscillation probabilities, 
N^^/N^a ^ 1 - 0.65523. For 0.004 eV^ the contribution of ^23 is more resulting in a 
lower number of muon events. For the electron events however the behavior with A23 is 
opposite, with Ne/N^o = 0.82 + 0.125*23. Thus with increasing A23 the number of electron 
events increase. Also note that since the contribution of ^23 comes with opposite sign the 
zenith-angle distribution for a fixed A23 is opposite for the muon and the electron events. 

In fig. ^]^(c) we show the histograms in the limit of sf^ = 0.0, keeping as 0.1 and 
varying A23 and S23. As S23 increases all the relevant probabilities decrease and therefore 
both N^/N^Q and Nf./Neo decrease giving less number of events for both. For this case 
the ^23 term comes with the same sign (negative) in both N^/N^q and Ne/Neo- Therefore 
the depletion is more for higher A23 for both muon and electron events. 

Finally, the long dashed line in fig. |5]^(d) represent the histograms for the best-fit 
value for two-generation z/^ — Ur oscillations, for which Pu^u^ = 1- The short dashed 
line gives the histograms for the three-generation best-fit values. Both give comparable 
explanation for the zenith angle distribution of the data. The dotted line gives the event 
distribution for A23 = 10~^ eV^. As discussed earlier even for such low value of A23, we 
find that due to the unique feature of the beyond OMSD neutrino mass spectrum, earth 
matter effects ensure that both the sub-GeV as well as the multi-GeV upward muon events 
are very well reproduced, as are the electron events. But since is high, the downward 

are depleted more than the data requires. 

5.2.2 Allowed Parameter Region 

In fig. the solid lines give the 99% CL. allowed area from SK data in the A23-SI3 
plane keeping the values of sf^ and fixed in the allowed range from fig. ^]l] and Table 



The first panel represents the two-generation z/^ — oscillation limit modulo the 
difference in the definition of the CL. limit as the number of parameters are different. 
We have seen from the histograms in fig. |5.2| (a) that raising 5^2 results in electron excess 
and muon depletion. On the other hand increase in s^g causes electron depletion. The 
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Figure 5.4: The allowed parameter regions in the A23 — S23 plane for various fixed values 
of "5i3; shown at the top of each panel. The solid lines corresponds to the 99% 

CL. contours from the SK data alone, while the dotted line gives the 99% contour from 
the combined analysis of the SK+CHOOZ data. 



above features are refiected in the shrinking and disappearance of the allowed regions in 
the first row and column. In the panels where both ^13 nonzero one may get 

allowed regions only when the electron depletion due to increasing s^g is replenished by 
the increase in 5^2- 

In fig. ^.5| we present the 99% CL. allowed areas in the bilogarithmic tan^ ^12 — tan^ 613 
plane for various fixed values of the parameters A23 and S23. We use the log(tan) repre- 
sentation which enlarges the allowed regions at the corners and the clarity is enhanced. 
The four corners in this plot refer to the two-generation limits discussed before. The 
extreme left corner {612 0, ^^13 0) correspond to the two generation z/^ — Ur oscillation 
limit. As we move up increasing ^13, one has i/g ~ J^^i and z/g — mixing in addition and 
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Figure 5.5: Same as fig. O but in the bilogarithmic tan^ 612 
values of A 23 and S23. 
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for 1 one goes to the two generation z/^ — oscillation region. For the best-fit 

values of A23 and S23 if we take s^2 sf^ to be and 1 respectively, then the Xmin 
66.92 which is therefore ruled out. Both the right hand corners in all the panels refer to 
pure Ue — Vt oscillations and therefore there are no allowed regions in these zones. For 
the panels in the first row, A23 = 0.006 eV^ and the 2-generation — Vj. oscillation limit 
is just disallowed. The small area allowed for the middle panel of first row (between the 
solid lines) is due to the fact that for non-zero s^2 ^'^d ■^^3 the electron events are better 
reproduced, while S23 = 0.5 takes care of the muon events. Hence for this case slight 
mixture of — Vf. and z/g — oscillations is favoured. This feature was also reflected in 
the fact that in the fig. ^.4| , the panel for 5^2 = 0.1 and = 0.3 has more allowed range 
for A23 than the panel for the 2-generation — limit. For the panels with A23 = 0.002 
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Figure 5.6: Same as fig. p]^but in the s 
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plane for fixed values of sf^ and A23 



eV^, both the pure — Ur limit as well as full three-generation oscillations, give good 
fit. For the last two rows with A23 = 0.0007 eV^ and 0.0004 eV^ the matter effects are 
important in controlling the shape of the allowed regions. Infact the allowed region that 
one gets for 0.0004 eV^ and S23 = 0.5 is the hallmark of the matter effect in this particular 



three-generation scheme. As can be seen from fig. ^.l| (a) and Table |5]l|, if one does not 
include the matter effect, then there are no allowed regions below A23 = 0.0005 eV^ for 
any arbitrary combination of the other three parameters. Even for the first and the last 
panels with A23 = 0.0007 eV^, one gets allowed areas solely due to matter effects. 

In fig. ^.6| the solid lines show the 99% CL. allowed regions from SK data in the 
S23 — 5^2 plane for fixed values of A23 and 5^3. In contrast to the previous figure, here 
(and in the next figure) we use the sin — sin representation because the allowed regions are 
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around ^^23 = 7r/4 and this region gets compressed in the log(tan)— log(tan) representation. 
For explaining the various allowed regions we separate the figures in two sets 

• For = 0.0, the four corners of the panels represent the no-oscillation limits 
inconsistent with the data. Also as discussed in the earlier sections for S23 = 0.0 or 1.0 
one goes to the limit of pure z/^ — z/g conversions driven by Alsnd, which is not consistent 
with data. One obtains allowed regions only when S23 is close to 0.5 with 5^2 small, so that 
^'^ ~ conversions are dominant. The allowed range of is controlled mainly by the 
electron excess as has been discussed before while the allowed range of S23 is determined 
mostly by the muon depletion. 

• For 7^ 0, the four corners represent the two-generation Ue — i't oscillation limit 
and hence these corners are not allowed. For S23 = 0.0 or 1.0 and 7^ or 1 one has 
'Alsnd driven z/^ — z/g and — Vr conversion and Aatm driven Ve — conversions. This 
scenario is not allowed as it gives excess of electron events and also fails to reproduce the 
correct zenith angle dependence. For a fixed A23 as increases the electron depletion 
increases which can be balanced by increasing 5^2 which increases the number of electron 
events. Hence for a fixed A23 the allowed regions shift towards higher 5^2 values. 

As in fig. |5.5| the allowed area in the middle panel of the last row is due to the inclusion 
of the matter effect. 



In fig. |5.7] the solid contours refer to the 99% CL. allowed areas from SK atmospheric 
neutrino data in the s^3 — S23 plane for various values of A23 and s^2- 

• For s^2 — 0-0 ^he corners represent no oscillation limits. In the limit S23 ^ or 1, 
one gets u^. — oscillation driven by Alsnd which is also not allowed. For 5^3 = 0.0 and 
S23 ~ 0.5 one has maximal two-fiavour z/^ — oscillation limit which is therefore allowed 
(not allowed for A23 = 0.006 eV^ as discussed before). As s^g increases the electron 
depletion becomes higher and that restricts higher 5^3 values. 

• For s^2 7^ O5 fo^^ corners represent two-generation limits driven by A^snd- This 
is the regime of average oscillations and cannot explain the zenith angle dependence of 
the data. For a fixed A23 the allowed region first expands and then shrinks in size and 
also shifts towards higher 5^3 values as 5^2 increases just as in fig. 
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Figure 5.7: Same as fig. |5]^ but in the s^g — s^^ plane for various fixed values of 

A23. 



Matter effect is important for the last two rows and the increase in the allowed areas 
for the last two panels of A23 = 0.0004 eV^ are typical signatures of matter effect. 

In fig. |5]^ we present the allowed range in the A23 — S23 plane with A23 in the 
10^^ — 10~^ eV^ range and s^2' ^13 fixed at 0.185 and 0.372 respectively. We get allowed 
regions in this range of small A23 and small mixing due to matter effects - a feature 
unique to the mass spectrum considered in this chapter. 

5.3 analysis of the SK + CHOOZ Data 

The CHOOZ experiment can probe upto 10~^ eV^ and hence it can be important 
to cross-check the atmospheric neutrino results. In particular a two-generation analysis 
shows that CHOOZ data disfavours the z/^ — z/g solution to the atmospheric neutrino 
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Figure 5.8: The allowed parameter space in the A23 — S23 plane with A23 in the range 
10-5 - 10"^ eV2 and with fixed values of sj^ = 0.185 and s^^ = 0.372. 



problem. The general expression for the survival probability of the electron neutrino in 
presence of three fiavours is 



P.,., = 1 - AU^,{1 - U^,)sm\7rL/Xu) - AU^,U^,sm'{7TL/ X 



23) 



(5.25) 



This is the most general expression without the one mass scale dominance approximation. 
We now minimize the defined as 



2 _ 2 I 2 
X — Xatm ' XcHOOZ 



(5.26) 



where Xatm calculated as before using eq. ( p.62| ) and we define x'^chooz as |2^ 



XcHooz — X! 



j=l,15 ^•'^3 



(5.27) 



where Xj are the experimental values, yj are the corresponding theoretical predictions and 
the sum is over 15 energy bins of data of the CHOOZ experiment [|l^. For the CHOOZ 
experiment the sin^(7rL/Ai2) term does not always average out to 0.5 (for SK this term 
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Figure 5.9: Same as fig. O but for smaller values of ^13^ chosen from the range 



determined by the SK+CHOOZ dashed line in fig. 5.1. 



always averages to 0.5) and one has to do the energy integration properly. For our analysis 
we keep the A12 fixed at 0.5 eV^ and do a four parameter fit as in SK. The xLin ^^e 
best-fit values of parameters that we get are 

• xLJd-o.f. = 42.22/51, A23 = 0.0023 eV^, 5^3 = 0.5, sj^ = 0.0022 and = 0.0. 

Thus the best-fit values shift towards the two-generation limit when we include the 
CHOOZ result. This provides a very good fit to the data being allowed at 80.45% prob- 
ability. 



The dotted lines in fig. 5.1 give the combined SK-I-CHOOZ Ax (= X Xmin) given by 



eq. (|5.26| ), as a function of one of the parameters, keeping the other three unconstrained. 



We find that the CHOOZ data severely restricts the allowed ranges for the parameters 
and 3^3 to values ~ 0.047, while A23 and S23 are left almost unaffected. Since CHOOZ 
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is consistent with no oscillation one requires Pu^ue close to 1. So the second and the third 
terms in eq. ( p. 251 ) should separately be very small. The second term implies U^i to be 



close to either or 1. U^^ close to zero implies either '^^ ^13 close to 1 which is not 
consistent with SK. Therefore U^i is close to 1. Then from unitarity both U^2 ^^cl ?7g3 
are close to and so the third term goes to zero irrespective of the value of A23 and S23. 
Hence contrary to expectations, CHOOZ puts almost no restriction on the allowed values 
of S23 and A23, although A23 ~ 10^^ eV^ - in the regime in which CHOOZ is sensitive. 
On the other hand it puts severe constraints on the allowed values of s^2 ■^13 order 
to suppress the average oscillations driven by A12. Because of such low values of 3^2 and 
s^3 the matter effects for the atmospheric neutrinos are not important and the additional 
allowed area with low A23 and high S23 obtained in the SK analysis due to matter effects 
are no longer allowed. The 99% CL. regions allowed by a combined analysis of SK and 
CHOOZ data is shown by the dotted lines in figs. |5.4| - |5.7| . It is seen that most of the 
regions allowed by the three-flavour analysis of the SK data is ruled out when we include 



the CHOOZ result. None of the allowed regions shown in fig. are allowed excepting 
the two-generation z/^ — Ur oscillation limit because CHOOZ does not allow such high 
values of either 5^3 or 5^2- Hence we present again in fig. |5.9| the allowed regions in the 
A23 — S23 plane for various fixed values of s^2 ^"^^ -^135 determined from the dotted lines 
in fig. |5]l|. The solid lines in fig. ^]9| give the 99% C.L. area allowed by the SK data 
while the dotted lines give the corresponding allowed region from the combined analysis 
of SK-I-CHOOZ. We find that for the combined analysis we get allowed regions in this 
plane only for much smaller values of s^2 ^13^ which ensures that the electron events 
are neither less nor more than expectations. 

5.4 Combined Allowed Area from Short Baseline Accelerator 
and Reactor Experiments 

As mentioned earlier the higher mass scale of this scenario can be explored in the 
accelerator based neutrino oscillation search experiments. For the mass-pattern considered 



the most constraining accelerator experiments are LSND 0, CDHSW |2^, E531 



and KARMEN ||28|. Among these only LSND reported positive evidence of oscillation. 
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Other experiments are consistent with no-oscillation hypothesis. Also important in this 
mass range are the constraints from the reactor experiment Bugey [B^|. The relevant 



probabilities are [|T^ 

• Bugey 

Pp,^, = 1 - 4c?34sin2(7rL/Ai2) + 444sin2(7rL/Ai2) (5.28) 

• CDHSW 

.2 „2 „• 2 



P^^^^ = 1 - 4cf2sf2sin^(7rL/Ai2) (5.29) 



LSND and KARMEN 



E531 



Pv,v^ = ^ci^si,ci^sm'{7TL/Xu) (5.30) 



Aci^si,si.,sm'{7rL/Xrs) (5.31) 



We note that the probabilities are functions of one of the mass scales and two mixing an- 
gles. Thus the one mass scale dominance approximation applies. There are many analyses 
in the literature of the accelerator and reactor data including LSND under this one mass 
scale dominance assumption [|^, ^ . These analyses showed that when one considers the 
results from the previous (prior to LSND) accelerator and reactor experiments there are 



three allowed regions in the 9i2 — 9i3 plane [|T4|, ^ 

• low 9i2 - low 9i3 

• low 9i2 - high 

• high 9i2 - ^13 unconstrained 

When the LSND result was combined with these results then only the first and the third 
zones remained allowed in the mass range 0.5 < A12 < 2 eV^. In these earlier analyses of 



the accelerator and reactor data p4| , pO| E776 |^ was more constraining than KARMEN. 
But with the new data KARMEN2 gives stronger constraint than E776. Also the results 
from the KARMEN2 experiment now rule out most of the region allowed by the LSND 
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Figure 5.10: The area between the dashed hnes is the 90% CL. region allowed by LSND 
while the light shaded zone gives the 90% CL. allowed region from the non-observance 
of neutrino oscillation in the other short baseline accelerator and reactor experiments 
except KARMEN2. The corresponding area which includes KARMEN2 as well is shown 
by the region shaded by asterix. The 90% CL. allowed region from SK+CHOOZ analysis 
is within the dotted line. The dark shaded area corresponds to the combined allowed 
region. 



experiment above 1 eV^ |2^. The LSND collaboration has also now done a reanalysis 
of their entire data sample and report their final transition probability in [^]. We have 
repeated the analysis with the latest LSND and KARMEN results for one representative 
value of Ai2 = 0.5 eV^ and present the allowed region in fig. |5.1U . 

The light-shaded area in fig. |5.10| shows the 90% CL. allowed area in the bilogarith- 
mic tan^ 612 — tan^ ^13 plane from the observance of no-oscillation in all the other above 
mentioned accelerator and reactor experiments except KARMEN2. The inclusion of the 
KARMEN2 results as well gives the 90% CL. region shown by the area shaded by asterix. 
The 90% allowed region by the LSND experiment is within the dashed lines. The KAR- 
MEN2 data severely restricts the LSND allowed regions. The dotted line shows the 90% 
CL. (x^ < Xmin + 7.78) region allowed by the combined analysis of the SK+CHOOZ 
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data keeping A23 and s^g at 0.002 eV^ and 0.5 respectively. The combined SK atmospheric 
and the CHOOZ reactor data rule out the third zone (high 612 with 6*13 unconstrained 
) allowed from LSND and other accelerator and reactor experiments. Thus if one takes 
into account constraints from all experiments only a small region in the first zone (small 
^12, ^13) remains allowed. This common allowed region is shown as a dark-shaded area in 
the fig. [5.10|. As evident from the expression of the probabilities for the accelerator and 



reactor experiments the combined allowed area of all the accelerator reactor experiments 
remains the same irrespective of the value of A23 and S23. Even though the combined 
area in fig. ^.10| shows that in the first zone (small 6*12,6*13), SK+CHOOZ data allows 
more area in the 6^12 — 6^13 plane for A23 = 0.002 eV^ and S23 = 0.5, from fig. ^]5| we see 



that for some other combinations of A23 and S23 one does not find any allowed zones from 
the SK+CHOOZ analysis, even at 99% C.L.. For those sets of values of A23 and S23 the 
SK+CHOOZ analysis is more restrictive than the LSND and other accelerator reactor 
data. 

5.5 Implications 

From our analysis of the SK atmospheric data the explicit form for the 3x3 mixing 
matrix U at the best-fit values of parameters is 



U 



( 0.95 -0.039 0.31 

-0.2 0.686 0.7 I (5.32) 



\-0.24 -0.727 0.644 
From the combined SK+CHOOZ analysis the mixing matrix at the best-fit values of 
the parameters is 

0.999 0.033 0.033 \ 
U = I -0.047 0.706 0.706 (5.33) 
-0.0 -0.707 0.707/ 



From the combined allowed area of fig. |5.10| the mixing matrix at A12 = 0.5 eV^, 



A23 = 0.0028 eV\ = 0.005, sj^ = 0.001 and 5^3 = 0.5, is 



0.997 0.028 0.072^ 
U = I -0.071 0.705 0.705 | (5.34) 
-0.032 -0.708 0.705 
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Figure 5.11: 90% C.L. regions in the A23 — S23 plane that can be explored by the i'^ — 
(solid line) and v^ — Ve (dashed line) oscillation channels in the K2K experiment. The area 
inside the dotted line shows the 90% C.L. region allowed by SK+CHOOZ. The curves are 
presented for fixed values of and s\-^ with A12 = 0.5 eV^. 




Thus the allowed scenario corresponds to the one where (z/i|z/e) is close to 1 while the 
states z/2 and are combinations of nearly maximally mixed and v^- 

Long baseline (LBL) experiments can be useful to confirm if the atmospheric neutrino 
anomaly is indeed due to neutrino oscillations, using well monitored accelerator neutrino 
beams. Some of the important LBL experiments are K2K (KEK to SK, L 250 km) 

MINOS (Fermilab to Soudan, L ^ 730 km ) [||] and the proposed CERN to Gran 

Sasso experiments (L 730 km) ||3^ . In this section we explore the sensitivity of the LBL 

experiment K2K in probing the parameter spaces allowed by the SK+CHOOZ and other 

accelerator and reactor experiments including LSND. K2K will look for disappearance 

Thus this scenario is the same as the one termed 3a in Table VI in the pre-SK analysis of In 
their notation the states 2 and 3 were 1 and 2. It was disfavoured from solar neutrino results. 

^'*K2K has already presented some preliminary results |Q. 
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as well as Ve appearance. In fig. p. 11 we show the regions in the A23 — S23 plane that can 



be probed by K2K using their projected sensitivity from ||3^. The top left panel is for the 

-.2 

'13 



two-generation — limit. The other panels are for different fixed values of sfg and s\ 



while A12 is fixed at 0.5 eV^. For LBL experiments the term containing A12 averages to 
0.5 as in the atmospheric case. The solid lines in the panels show the region that can be 
probed by K2K using the disappearance channel while the dotted lines give the 90% 
C.L. contours allowed by SK+CHOOZ. One finds that for for A23 > 2 x lO^^ eV^, the 
whole region allowed by SK+CHOOZ can be probed by the z/^ disappearance channel in 
K2K. The dashed lines show the 90% C.L. area that K2K can probe by the v^, appearance 
mode. As 5^2 increases the constraint from the Pv^^v^ channel becomes important as is seen 
in the top right panel of fig. |5.11| . However such high values of 5^2? although allowed by 



SK+CHOOZ, is not favoured when one combines LSND and other accelerator and reactor 
results. For lower 5^2 values allowed by all the accelerator, reactor and SK atmospheric 
neutrino experiment the projected sensitivity in the z/^ — z/g channel of K2K is not enough 
to probe the allowed regions in the A23 — S23 plane as is shown by the absence of the 
dashed curves in the lower panels. 

In fig. |5.12| we show the regions in the bilogarithmic tan^ Qyi — tan^ 6*13 plane which 
can be probed by K2K. For drawing these curves we fix A23 = 0.002 eV^, S23 = 0.5 and 
A12 = 0.5 eV^. Shown is the area that can be explored by the — z/^ (left of the solid 
line) and — v^, (hatched area) channels in K2K at 90% C.L.. The light-shaded area is 
allowed by SK+CHOOZ and the dark shaded area is allowed by the combination of all the 
accelerator, reactor and SK atmospheric neutrino data at 90% C.L.. It is clear from the 
figure that even though the sensitivity of the z/g appearance channel is not enough, K2K 
can still probe the combined allowed region in the Q\2 — ^13 plane from z/^ disappearance. 

The projected sensitivities of MINOS and the CERN to ICARUS proposals are lower 
than K2K and it will be interesting to check if one can probe the regions allowed in this 
picture better in these experiments. However since in our case the OMSD approximation 
is not applicable one has to do the energy averaging properly to get the corresponding 
contours in the three-generation parameters space, and one cannot merely scale the al- 
lowed regions from the two-generation plots. For K2K we could use the fig. 5 of 
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Figure 5.12: Sensitivity of the K2K experiment in the tan^ ^12 — tan^ ^13 plane for 
A23 = 0.002 eV^, S23 = 0.5 and A12 = 0.5 eV^. The area that can be explored by the 
i^ti ~ (left of solid line) and — Ve (hatched area) channels in K2K at 90% CL. is 
shown. The light-shaded area is allowed by SK+CHOOZ and the dark-shaded region is 
the combined area allowed by all accelerator and reactor data at 90% C.L.. 

to circumvent this problem. However since the analogous information for MINOS and 
CERN-Gran Sasso proposals is not available to us we cannot check this explicitly. 

An important question in this context is whether one can distinguish between the 
OMSD three generation and this mass scheme. In both pictures the SK atmospheric 
neutrino data can be explained by the dominant z/^ — oscillations mixed with little 
amount of — ^pii^r) transition. However the mixing matrix U is different. A distinction 
can be done if one can measure the mixing angles very accurately. 

What is the prospect in LBL experiments to distinguish between these pictures? We 
give below a very preliminary and qualitative discussion on this. If we take — 0.02, 
s?3 = 0.02 and s^g = 0.5, P^^^^ would be (0.038 + 0.0004 (^23))- As the second term 
is negligible one has average oscillations. This is different from the OMSD limit where 
Py^y^ = iU^^U^^^'is is energy dependent. If one combines the other accelerator and reactor 
experiments including LSND then the allowed values of of ^13 even less and 
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choosing sl^ = 0.005, sj^ = 0.001 and s^^ = 0.5 we get P^^^^ = 0.01 - 0.004(,S23). Here 
also the term involving (S'23) is one order of magnitude smaller and the oscillations will 
be averaged. Thus this channel has different predictions for the OMSD limit and beyond 
the OMSD hmit. 

5.6 Discussions and Conclusions 

In this chapter we have done a detailed analysis of the SK atmospheric neutrino 
data going beyond the OMSD approximation. The mass spectrum chosen is such that 
A12 = Ai3 ~ eV^ to explain the LSND data and A23 is in the range suitable for the 
atmospheric neutrino problem. We study in details the implications of the earth matter 
effects and bring out the essential differences of our mass pattern with the OMSD scenario 
and the two-generation limits. 

We first examine in detail what are the constraints obtained from only SK data con- 
sidering its overwhelming statistics. The allowed regions include 

• the two-generation i/^ — 1/^- limit (both *i3 zero) 

• regions where either s^2 or s^g is zero; in this limit the probabilities are functions in 
general of two mixing angles and two mass scales. 

• the three-generation regions with all three mixing angles non-zero and the proba- 
bilities governed by both mass scales. 

The last two cases correspond to dominant i/^ — Ur oscillation with small admixture 
of Vfj, — Vf. and — oscillation. 

• regions with very low A23 (< 10~^ eV^) and close to 1, for which the earth matter 
effects enhance the oscillations of the upward neutrinos and cause an up-down flux 
asymmetry. This region is peculiar to the mass spectrum considered by us and is 
absent in the two-generation and the OMSD pictures. 

We present the zenith angle distributions of the events in these cases. With the inclusion 
of the CHOOZ result the allowed ranges of the mixing angles ^^<i •^is constrained 
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more (~ 0.047), however the allowed ranges of A23 and do not change much (see fig. 

except that the low A23 region allowed by SK due to matter effects is now disallowed. 
The inclusion of the constraints from LSND and other accelerator and reactor experiments 
may restrict the allowed area in the ^12 — ^13 plane for certain values of A23 and S23, but for 
some other combinations of A23 and S23, SK+CHOOZ turns out to be more constraining. 
We have included the latest results from LSND and KARMEN2 in our analysis. 

In order to explain the solar neutrino problem in this picture one has to add an extra 
light sterile neutrino. With the new LSND results the allowed 4 neutrino scenarios are 

• the (2+2) picture where two degenerate mass states are separated by the LSND gap 
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the (3+1) scheme with three neutrino states closely degenerate in mass and the 



fourth one separated from these by the LSND gap [p9l ffO, 41 



Our scenario can be easily extended to the (2+2) scheme which is still allowed after the 
inclusion of the SNO results The (3+1) scheme is however shown to be ruled out 



from the atmospheric data in ET 



To conclude, one can get allowed regions from the SK atmospheric neutrino data where 
both the mass scales and all the three mixing angles are relevant. The beyond one mass 
scale dominance spectrum considered in this chapter allows new regions in the low mass 
- low mixing regime due to the earth matter effects. With the inclusion of the CHOOZ, 
LSND and other accelerator reactor results, the allowed regions are constrained severely. 
It is, in principle, possible to get some signatures in the LBL experiments to distinguish 
this picture from the OMSD limit. 
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CHAPTER 6 



Decay of Atmospheric Neutrinos 



In the previous chapters we have seen that neutrino flavor oscillations in vacuum, both 
in two as well as three flavors, gives an excellent flt to the SK atmospheric neutrino data. 
However there are other possibilities 0, ^ |, ||, ^, that have the potential to explain 
the data, though may be not as well as — oscllations do. Nevertheless one has to 
rule out all these possibilities before coming to any deflnite conclusion about the fate of 
the atmospheric neutrinos. One among these is neutrino decay 0. 

One aspect peculiar to the oscillation hypothesis is the periodicity of the resultant 
neutrino beam at the detector. But so far one has no direct experimental evidence for 
this periodicity. The other way to remove any ambiguity is to directly observe the number 
of Vr predicted by the z/^ — oscillation scenario. The latest 1289 day SK data for the 
upward-going sample is consistent with i/^ appearance at the 2a level. But one needs 



more statistics before coming to any flnal conclusion. 

If one assumes the existence of a characteristic wavelength A for the vanishing such 
that A has a power law dependence on energy (A oc E'"") then the — conversion 



probability can be assumed to be parameterized as O, 12 



P,^,^ = asm^PLE^) (6.1) 

where a, (3 and n are free parameters. If then one performs a most general analysis, 
keeping the energy exponent n a free parameter, the best-flt comes out for n ^ — 1 



TT|, Q2[. In ref. |11] Fogli et al. further corroborate this feature by showing that only an 



L/E distribution can flt the muon zenith angle data ||T3|. Hence it can be conclusively 



said that the data strongly prefers theoretical models which predict a, L/ E dependence for 
the z/^ survival probability. Oscillations in vacuum are hugely favored since they predict 
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this L/E behavior. But so does decay and it will be interesting to check if neutrino decay 
can really offer a decent fit to the observed atmospheric neutrino deficit. 

In ref. it was shown that neutrino decay gives a poor fit to the data. However they 
considered neutrinos with zero mixing. Barger et al. considered the situation of neutrino 
decay in the general case of neutrinos with non-zero mixing angle [^]. They showed that 
the neutrino decay fits the L/E distribution of the SK data well. The Am^ taken by them 
was > 0.1 eV"^ so that the Am? dependent term in the expression for the neutrino survival 
probability averages out to zero. As pointed out in |^ such a constraint on Am? is valid 
when the unstable state decays into some other state with which it mixes. If however 
the unstable state decays into a sterile state with which it does not mix then there is no 
reason to assume Am^ > 0.1 eV^. 

In this chapter we present our results of the neutrino decay solutions to the atmo- 
spheric neutrino problem by doing x^"fit to the 848 day of sub-GeV and multi-GeV 
Super-Kamiokande data fl^. We also present the results of x^-fit to the 535 day SK 



data and compare it with the results for the 848 day data. For the we use the defini- 
tion given by eq. ( |3.61| ) in chapter 3, where we use the double ratio R and the up-down 
asymmetry parameter {a = e,fi) [|1^, |T^. We first present the results of this fit for 



the two-generation z/^ — z/,- oscillations. For the neutrino decay analysis we take the most 
general case of neutrinos with non-zero mixing and consider two pictures 

• Am^ > 0.1 eV"^ (scenario (a)) 

• Am^ unconstrained (scenario (b)) 

We also explicitly demonstrate the behavior of the up-down asymmetry parameters in 
both scenarios. 

Our analysis shows that scenario (a) is ruled out at 100%(99.99%) C.L. by the 848(535) 
day of SK data. However if we remove the constraint on Am^ and consider the possibility 
of decay into a sterile state then one can get an acceptable fit for Am^ ~ O.OOle^^ and 
sin^ 29 large. 

In section 6.1 we present the most general expression for the i/^ survival and transition 
probabilities with unstable neutrinos. We first consider neutrinos to be stable and display 
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our results for two-generation — u^. oscillation analysis in section 6.2.1. In section 6.2.2 
we present our results for the neutrino decay solution constraining Am^ to be > Q.leV'^. 
In section 6.2.3 we do a three parameter analysis by removing the constraint on Am^. 
In section 6.3 we conclude by performing a comparative study of the three cases and 
indicate how one can distinguish experimentally between the scenario (b) and the z/^ — Vr 
oscillation case though both give almost identical zenith-angle distribution. 

6.1 Neutrino Survival Probabilities 

Neutrinos are assumed to be stable in the standard model of particle physics. But 
if one allows for neutrino decay then the the analysis of neutrino oscillation experiments 
become quite different. We will here assume that the only unstable component is V2 which 
decays into some other lighter state Vj which may be an active or sterile species. Since 



radiative decays of neutrinos are severely constrained [0 we consider the two possible 
non-radiative decay modes discussed in the literature. 

• Model 1: If neutrinos are Dirac particles one has the decay channel V2 ^jR + 
where VjR is a right handed singlet and is an iso-singlet scaler. Thus all the final 
state particles for this model are sterile and there is no distinct signature of this 
decay apart from in disappearance experiments. This model is discussed in |TB|. In 



this model a light scalar boson with lepton number -2 and a singlet right handed 
neutrino is added to the standard model. The neutrino coupling to this scalar boson 
is given by g2jv]i.C~^i'ji^, C being the charge conjugation operator. 

• Model 2: If neutrinos are Majorana particles, the decay mode is z/2 Vj + J, where 
J is a Majoron, produced as a result of spontaneous breaking of a global U{X)L^-L^, 
symmetry [^]. In this model the neutrino masses are generated by extending the 
higgs sector of the standard model. 

In both the decay scenarios the rest frame lifetime of is given by 

-0 - -T (6-2) 
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where g is the couphng constant, is the i>i mass and = rn^ — rn^, the mass squared 
difference between the states that are involved in decay. Assuming 1712 » rrij the equation 
r2D can be written as 



g'^I^d ~ 167ra (6.3) 
where a is the decay constant related to tq as a = m2/ro. We assume a scenario where 

Ve ~ (6.4) 

~ ^2 COS 6 + sin 6 (6.5) 

From eq. ( |6.5| ) the survival and transition probabilities of the of energy E, with an 
unstable component z/2, after traveling a distance L is given by, 

P,^,^ = sin^0 + cos'^exp(-47rL/Ad) 

+ 2 sin^ e cos^ exp(— 27rL/Ad) cos(27rL/Aosc) , (6.6) 



P,^,^ = ^ sin2 2^{1 + exp(-47rL/Ad) - 2 exp(-27rL/Arf) cos(27rL/Ao,,)} (6.7) 

where A^ is the decay length (analogous to the oscillation wavelength Xosc given by eq. 
and is defined as, 

A, = 2.47km-^— (6.8) 
GeV a 



We see that the neutrino survival probability (|6.6|) depends on the decay constant a 
and the mass squared difference Am^ between the states that mix, apart from the mixing 
angle 9. Hence one may consider either a = 0, which would give pure z/^ — Vr oscillations 
of stable neutrinos or a non-zero, which would corresponds to the case of decay along 
with oscillations. The latter may again be subdivided into two cases depending on the 
value of Am^ : (a) Am^ > 0.1 eV^ and (b) Am^ unconstrained. We present below the 
results for the fits to all the three above mentioned cases. The data used for R and Y 
is shown in Table |3.5| in chapter 3. 
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6.2 Results of the analysis 

6.2.1 Two-Generation — v^. Oscillations 

For the two flavor z/^ — Vj. oscillations with the defined as in eq. ( p.61| ) of chapter 
3, the Xmm that we get for the 848 day data is 1.21 with the best-fit values as Am^ = 
0.003 eV"^ and sin^ 2^ = 1.0. With six data points and two parameters this provides a 
good fit to the data being allowed at 87.64%. If we use the 535 day data then the Xmm 
that we get is 4.25 with the best-fit values as /S.rn? = 0.005 eV'^ and sin^ 20 = 1.0, the 
g.o.f being 37.32%. Thus the fit becomes much better with the 848 day data with no 
significant change in the best-fit values. Though we have used a different procedure of 
data fitting, our results agree well with that obtained by the SK collaborationj^. 

6.2.2 Neutrino Decay with Am^ > 0.1 eV^ 

If the unstable component in the state decays to some other state with which 
it mixes then bounds from K decays imply Am^ > 0.1 eV^ [^. In this case the 
cos(27rL/Aosc) term averages to zero and the probability becomes 

P^^^^ = sin^ e + cos^ 9 exp(-47rL/Ad) . (6.9) 



In figs. ^]1] and |6.2| we show the variation of R and Y with a for various values of sin^ 9 for 
the sub-GeV and multi-GeV cases. For higher values of a, the decay length given by 
eq. (^) is low and the exponential term in the survival probability is less implying that 



more number of neutrinos decay and hence R is low. As a decreases the decay length 

increases and the number of decaying neutrinos decreases, increasing R. For very low 

values of a the exponential term goes to 1, the neutrinos do not get the time to decay so 

that the probability becomes 1 — | sin^ 29 and remains constant thereafter for all lower 

values of a. This is to be contrasted with the — Vj. oscillation case where in the no 

oscillation limit the sin^(7rL/Aosc) term and the survival probability — * 1. For multi- 

GeV neutrinos since the energy is higher the A^ is higher and the no decay limit is reached 

for a larger value of a as compared to the sub-GeV case. This explains why the multi-GeV 

i^The best-fit values that the SK coUaboration had got for the 848 day data are Q Am^ = 0.003 eV^, 
sin^ 29 = 0.995 and = 55-4 for 67 d.o.f. This corresponds to a g.o.f of 84.33%. 
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Figure 6.1: The variation of R and Y with a for the sub-GeV neutrinos (denoted by 
the subscript sg) assuming neutrino decay with Am^ > O.leV^. The curves are drawn at 
fixed values of sin ^=0.03 (sohd hue), sin2^=0.04 (long dashed line), sin2^=0.06 (short 
dashed line), sin^ ^=0.08 (dotted line) and sin^ 9=0.1 (long dashed-dotted line). The short 
dashed-dotted lines give the SK 848 day results within a ±lcr band. Also shown by a 
cross are the R and Y at the best-fit point. 



curves become flatter at a higher a. The behavior of the up-down asymmetry parameter 
is also completely different from the only oscillation case |^T|. In particular the plateau 
obtained for a range of Am^ which was considered as a characteristic prediction for up- 
down asymmetries is missing here. For the decay case even for a as high as 0.001 eV"^, the 
decay length = 2500 (E/GeV) km so that the exponential term is 1, there is almost no 
decay for the downward neutrinos and the survival probability is P = 1 — | sin^ 29 while 
the upward going neutrinos have some decay and so Y is less than 1. As a decreases, 
the Xd increases, and the fraction of upward going neutrinos decaying decreases and this 
increases Y. For very small values of a even the upward neutrinos do not decay and Y —>■ 
1 being independent of 9. 

We also perform a analysis of the data calculating the "th" quantities in ( p.61| ) for 
this scenario. For the 848 day data the best-fit values that we get are a = 0.33 x 10~^ in 
eV"^ and sin^ 9 = 0.03 with a Xmm of 49.16. For 4 degrees of freedom this solution is ruled 
out at 100%. The best-fit values for the 535 day of data that we get are a = 0.28 x 10~^ in 
eV"^ and sin^ 9 = 0.08 with a Xmm 31.71. For 4 degrees of freedom this solution is ruled 



out at 99.99% ||22|. Thus the fit becomes worse with the 848 day data as compared to the 
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Figure 6.2: Same as in fig. BTTI but for multi-GeV neutrinos. 



535 day data. We liave marked tlie R and Y corresponding to tfie best-fit value of tlie 
parameters a and sin^ 9 in figs. |6]l| and It can be seen that the best-fit value of R for 



the sub-GeV neutrinos is just below and that for the multi-GeV neutrinos is just above 
the ±1(7 allowed band of the SK 848 day of data. The up-down asymmetry parameter Y 
is quite low for the sub-GeV neutrinos and extremely high for the multi-GeV neutrinos as 
compared to that allowed by the data. The fig. |6.1| shows that for the sub-GeV neutrinos 



the data demands a lower value of a while from fig. |6.2| we see that the multi-GeV 
neutrinos need a much higher a to explain the SK data. It is not possible to get an a that 
can satisfy both the sub-GeV and the multi-GeV SK data, particularly its zenith angle 
distribution. In this scenario, decay for the sub-GeV upward neutrinos is more than that 
for the multi-GeV upward neutrinos (downward neutrinos do not decay much) and as a 
result Y for sub-GeV is lower than the Y for multi-GeV, a fact not supported by the data. 
Since the 848 day data needs even lesser depletion of the sub-GeV flux as compared to 
the multi-GeV flux, the fit gets worse. 

6.2.3 Neutrino Decay with Am^ unconstrained 

In this section we present the results of our x^-analysis removing the constraint on 
Am^. This case corresponds to the unstable neutrino state decaying to some sterile state 
with which it does not mix H]. The probability will be still given by eq. 
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Figure 6.3: The variation of R and Y with Am^ for the sub-GeV neutrinos (denoted by 
the subscript sg) assuming neutrino decay with Am? unconstrained. In these curves the 
a is fixed at its best-fit value of 0.3 x 10~^el^^. The curves are drawn at fixed values of 
sin^^=0.7 (dotted line), sin^6'=0.6 (short dashed line) and sin^6'=0.5 (long dashed line). 
The solid lines give the curves for the best-fit value (sin^ 9 = 0.5) of the — Vj. oscillation 
case. The dotted-dashed lines give the SK 848 day results within a ±la band. Also shown 
are the R and Y at the best-fit point. 



In fig. |6.3| and |6.4| we plot the R vs. Am and Y vs. Am for the sub-GeV and multi- 
GeV data for a = 0.3 x lO^^eV^ (which is the best-fit value we get for the 848 day data) 
and compare with the curve obtained for the best-fit value of sin^ ^ (=0.5) for the only 
oscillation case (solid line). For the best-fit value of a that we get, the downward neutrinos 
do not have time to decay while the upward neutrinos undergo very little decay. Thus the 
curves are very similar in nature to the only oscillation curves. In the sub-GeV case (fig. 
|6.3| ), for high values of Am^ around 0.1 eV"^ both upward and downward neutrinos undergo 
Am^ independent average oscillations and R stays more or less constant with Am^ . For 
the upward going neutrinos in addition to average oscillation there is little amount of 
decay as well and hence Y ~ Nup/Ndown is ~ 1- As Am^ decreases to about 0.05 eV"^ the 
oscillation wavelength increases - for upward neutrinos it is still average oscillation but 
for the downward neutrinos, the cos(27rL/Aosc) term becomes negative which corresponds 
to maximum oscillation effect and the survival probability of these neutrinos decreases, 
and hence R decreases; while the upward neutrinos continue to decay and oscillate at 
the same rate and Y becomes greater than 1. As Am^ decreases further, the downward 
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Figure 6.4: Same as in fig. but for multi-GeV neutrinos. 



neutrino oscillation wavelength becomes greater than the distance traversed and they are 
converted less and less and thus R increases and Y decreases. Below Am^ = 0.001 eV"^ 
the downward neutrinos stop oscillating completely while for the upward neutrinos the 
cos(27rL/Aosc) term goes to 1, and R and Y no longer vary much with Am^ . 



For the multi-GeV case (fig. |6.4|) the oscillation wavelength is more than the sub- 
GeV case and for Am^ around 0.1 eV"^ the cos{27t L/Xosc) term stays close to 1 for the 
downward neutrinos; while the upward neutrinos undergo average oscillations and slight 
decay and Y is less than 1. As Am^ decreases the downward neutrinos oscillate even 
less and the upward neutrinos also start departing from average oscillations and hence R 
increases and Y decreases. Around 0.01 eV'^ the downward neutrinos stop oscillation while 
for upward neutrinos the oscillation effect is maximum (A ~ L/2) and the cos{2ttL/Xosc) 
term is ~ -1 and Y stays constant with Am^ . As Am^ decreases further the upward 
neutrino oscillation wavelength increases and they oscillate less in number making both R 
and Y approach 1 for Am^ around 0.0001 eV"^. For multi-GeV neutrinos the decay term 
contributes even less as compared to the sub-GeV case. 

We perform a minimization in the three parameters Am^ , sin^ 29 and a. The 
best-fit values that we get for the 848 day data are Am^ = 0.003e\^^, sin^ 26' = 1.0 
and a = 0.3 x 10~^eV^^. The minimum that we get is 1.11 which is an acceptable 
fit being allowed at 77.46%. For the 535 day data the best-fit values that we get are 
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Am2 = 0.002e1/2, sin^ 26 = 0.87 and a = 0.0023e\/2 with a xLn of 4.14 which is allowed 
at 24.67%. Thus compared to the 535 day data, the fit improves immensely and the best- 
fit shifts towards the oscillation limit, the best-fit value of the decay constant a being 
much lower now. It is to be noted however, that the best-fit in this model does not come 
out to be a = 0.0, viz the only oscillation limit. In Table ^?1] we give the contributions to 

from the i?'s and Y^s at the best-fit value of a and for the a = 0.0 case. 

Thus from the contributions to we see that for the best- fit case there is improvement 
for the multi-GeV R and Y as compared to the a = 0.0 case. The fo^' sub-GeV Y also 
improves. In fig. |6.5| we plot Ax^ = ^ Xmm ^s. a with Am^ and sin^ 29 unconstrained. 
There are two distinct minima in this curve - one for lower values and another at higher 
values of a. The best-fit Am^ in both cases is ~ 0.001 eV"^. In this model there are two 
competing processes - oscillation and decay. For lower values of a the decay length is 
greater than the the oscillation wavelength and oscillation dominates. The decay term 
exp{—aL/E) is close to 1 and does not vary much with the zenith distance L. As a 
increases the exponential term starts varying very sharply with L and the variation is 
much more sharp for the sub-GeV as compared to multi-GeV. This behavior is inconsistent 
with the data and that is why one gets a peak in Ax^ for higher a. As a increases further 
the exp{—aL/E) term goes to zero for the upward neutrinos and there is complete decay 
of these neutrinos while the downward neutrinos do not decay, the exponential term still 
being 1. Whenever the exponential term is or 1 for the upward neutrinos, the wrong 



Quantity 



a 



0.3 X 10-5 eV^ 



a 



0.0 eV2 



e 

ymg 



0.085 
0.011 
0.48 
0.014 
0.344 
0.176 



0.021 
0.033 
0.56 
0.073 
0.344 
0.176 



Table 6.1: The various contributions to the x^min the best-fit value of a and at a=0.0 
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Figure 6.5: The Ax^ 
day (solid line) and 535 day data (dotted line). 



— Xmin ^s. a with Am^ and sin^ 29 unconstrained for the 84^ 



energy dependence of this term does not spoil the fit and these scenarios can give good 
fit to the data. Even though fig. |6.5| shows that the data allows a wide range of a, we get 
the two distinct minima in the Ax^ vs. a curve for high and low a values, for both the 
535 day (dotted line) and 848 day (solid line) data. But while the 848 day data prefers 
the lower a limit, the 535 day data gives a better fit for the high a limit. The reason 
behind this is that for the 848 day data the R is much higher than for the 535 day data. 
Hence the 848 day data prefers lower a and hence lower suppression. 

In fig. |6.6| we show the 90% and 99% CL. allowed parameter region in the Am^ - 
sin^ 29 plane for a range of values of the parameter a. In fig. |6.7| we show the 90% and 
99% CL. contours in the a - sin^ 29 plane fixing Am^ at different values. These contours 
are obtained from the definition ^ X^mm + ^X^; with Ax^ = 6.25 and 15.5 for the 
three parameter case for 90% and 99% CL. respectively. The bottom left panel in fig. 



Kdj is for the best-fit value of a. For high a (the top left panel) no lower limit is obtained 
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Figure 6.6: The allowed parameter region for the 848 day data in the Am? -sin^ 29 plane 
for 4 different values of a shown at the top of each panel. The solid and the dashed lines 
correspond to the area allowed at 90% C.L. and 99% C.L. respectively. 



on Am^, because even if Am^ becomes so low so that there is no oscillation the complete 
decay of upward neutrinos can explain their depletion. As we decrease a the allowed 
parameter region shrinks and finally for a = we get the two parameter limit modulo the 
small difference in the C.L. definitions for the two and three parameter cases. The upper 
right panel of fig. |6.7| corresponds to the best-fit value of Am^ . For very low a, even 
though there is no decay, we still have oscillations and that ensures that when Am^ is 
large enough there is no lower bound on a as evident in the fig. |6.7| . For Am^ = lO'^^eV"^ 
the neutrinos stop oscillating and hence we get a lower bound on a beyond which the 
depletion in the neutrino flux is not enough to explain the data. 
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Figure 6.7: The allowed parameter region for the 848 day data in the a-sin^ 26 plane for 
4 different values of Am? shown at the top of each panel. The solid and the dashed lines 
correspond to the area allowed at 90% C.L. and 99% C.L. respectively. 



6.3 Comparisons and Discussions 



In fig. |6.8| we show the histogram of the muon event distributions for the sub-GeV 
and multi-GeV data under the assumptions of — oscillation, and the two scenarios 
of neutrino decay for the best-fit values of the parameters both for the 535 and the 848 
day of data. From the figures it is clearly seen that the scenario (a) (big dashed line) 
(Am^ > O.ley^) does not fit the data well there being too much suppression for the sub- 
GeV upward going neutrinos and too less suppression for the multi-GeV upward going 
neutrinos. The scenario (b) (Am^ unconstrained, small dashed line), however, reproduces 
the event distributions well. However with the 848 day data the sub-GeV events are 
reproduced better as compared to the 535 day data and the quality of the fit improves. 
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The neutrino decay is an interesting idea as it can preferentially suppress the upward 
flux and can cause some up-down asymmetry in the atmospheric neutrino data. However 
the intrinsic defect in the decay term exp(— aL/i?) is that one has more decay for lower 
energy neutrinos than for the higher energy ones. Thus neutrino decay by itself fails 
to reproduce the observed data 0]. If however one considers the most general case of 
neutrinos with non-zero mixing then there are three factors which control the situation 

• the decay constant a which determines the decay rate 

• the mixing angle 9 which determines the proportion of neutrinos decaying and mix- 
ing with the other flavour 

• the Am? which determines if there are oscillations as well 

If the heavier state decays to a state with which it mixes then Am? has to be > Q.leV"^ 
because of bounds coming from K decays pO[. The best-flt value of a that one gets is 
0.33 X lO-^eV^ with the 848 day SK data. At this value of a the e""^/^ term tends to 1 
for the downward going neutrinos signifying that they do not decay much. The survival 
probability goes to (1 — ^ sin^ 29) which is just the average oscillation probability. In order 
to suppress this average oscillation the best-flt value of sin^ 9 comes out to be small in 
this picture. For the upward going neutrinos, in scenario (a), there will be both decay and 
average oscillations. If one had only average oscillation then the probability would have 
stayed constant for a flxed value of the mixing angle 9. But because of the exponential 
decay term the survival probability drops very sharply as we go towards cosO = —1.0. 
The drop and hence the decay is more for lower energy neutrinos. As a result the sub-GeV 
flux gets more depleted than the multi-GeV flux, a fact not supported by the data. In fact 
the 848 day data requires the sub-GeV flux to be even less suppressed than the multi-GeV 
flux as compared to the 535 day data and the flt worsens with the 848 day data. The 
small mixing signifles that the z/^ has a large fraction of the unstable component V2 (see 
eq. ( |6.5| )). Hence the constant a comes out to be low so that the decay rate is less to 
compensate this. However even at the best-flt a of 0.33 x lO^'^eV'^ the survival probability 
in the bin with cosG between —1.0 to —0.6 comes out to be 0.15 for E=l GeV, much 
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Figure 6.8: The sub-GeV and multi-GeV jj event distributions vs. zenith angle for the 
various scenarios considered. A^^ is the number of fi events as given by eq. (|3.58| ) and A^^o 
is the corresponding number with survival probability 1. The panels labelled SG(535) and 
MG(535) give the histograms for the sub-GeV and multi-GeV 535 day data respectively, 
while the SG(848) and MG(848) give the corresponding histograms for the 848 day data. 
For the both the sets the solid line corresponds to the best-fit — Vr oscillation solution, 
the long dashed line is for the best-fit for scenario (a) and the short dashed line for the 
best-fit for scenario (b). Also shown are the SK /x event distributions with ilex error bars 
for both the sets. 



lower than the value of ~ 0.5 as required by the data. Thus scenario (a) fails to explain 
the upward going neutrino data properly because of two main reasons 

• ^ is low in order to suppress the average oscillations of the downward neutrinos 

• the energy dependence of the exponential decay term is in conflict with the data 

In the scenario (b), in addition to mixing with z/^, the unstable component in decays 
to some sterile state with which it does not mix. In this case there is no restriction on 
Am? and it enters the fit as an independent parameter. We find that: 

• The best-fit Am? does not come out naturally to be in the Am^ independent average 
oscillation regime of > 0.1 eV"^, rather it is 0.003e\^^. 
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• The best-fit value of tlie decay constant a = 0.3 x 10 ^eV"^ implying that the decay 
rate is small so that the mixing angle is maximal (sin^ 9 = 0.5). 

• Large values of a giving complete decay of upward neutrinos are also allowed with 
a high CL. In fact with 535 day data the best-fit was in this region. 

• The best-fit value of the decay constant a is non-zero signifying that a little amount 
of decay combined with Am^ dependent oscillations gives a better fit to the data. 

At the best-fit values of the parameters there is no oscillation of the downward neutrinos so 
that the cos(27rL/Aosc) term goes to 1. The decay term also goes to 1 signifying that there 
is not much decay either for the downward neutrinos and the survival probability is ~ 1 
without requiring the mixing angle to be low. On the other hand for the upward neutrinos 
there are oscillations as well as little amount of decay. The sub-GeV upward neutrinos 
have smaller oscillation wavelength and they are close to the average oscillation limit 
(survival probability ~ 0.5) while for the multi-GeV neutrinos the oscillation wavelength 
is such that one has maximum oscillations and the survival probability is less than 0.5. 
Thus this scenario reproduces the correct energy dependence of the suppression - namely 
sub-GeV is suppressed less as compared to multi-GeV neutrinos. The best-fit value of a 
being even smaller now than the scenario (a) the decay term e""-^/^ does not vary very 
sharply with the zenith distance L or the energy E so that its wrong energy dependence 
does not spoil the fit. 

The conversion probability of to Vj. is given by eq. ( p.?] ). The value of P^^v^ 
integrated over the energy and the zenith angle, for a = 0.3 x lO'^eV"^ (the best fit for 
scenario (b)) is 0.33 for sub-GeV and 0.26 for multi-GeV. For a = 0.44 x lO'^eV^ (the 
second minima in the Ax^ vs. a curve) the corresponding numbers are 0.21 and 0.15, 
while for the only z/^ — Ur oscillation case, the corresponding values are 0.37 and 0.26 
respectively. The values of Am^ and sin^ 29 for all the cases are taken as 0.003ey^ and 
1.0 respectively. 

The fig. |6.8| shows that the zenith angle dependence of the scenario (b) is almost 
similar to the case of z/^ — Ur oscillation. But the two cases are very different in principle. 
For the oscillation case a larger 9 implies a larger conversion whereas in scenario (b) a 
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larger 6 means the fraction of the unstable component is less in z/^ and the depletion is 
less. If one compares the conversion probability as given by eq.( |6.7| ) with the one for 
the Ufj^ — v-r oscillation case, then the scenario (b) considered in this chapter would have 
smaller number of i/^s in the resultant flux at the detector, especially for the larger values 
of a which are still allowed by the data and the two cases might be distinguished when 
one has enough statistics to detect r appearance in Super-Kamiokande |]TD|, ^ or from 
neutral current events 



In this chapter we have followed the procedure of data fitting as done in [jl6l . Thus we 
use the ratios for which the common systematic errors get canceled out. Strictly speaking 
one should use the absolute number of events and include all the correlations between 
bins and e-like and /x-like events. But the best-fit points and the allowed regions are not 
expected to change significantly. We have compared the scenarios of neutrino oscillation 
and decay with the same definition of and for this purpose of comparison neglecting the 
correlation matrix will not make much difference. Apart from the statistical analysis we 
have given plots of R and Y for various values of the parameters. The allowed parameter 
ranges from these plots are consistent with what we get from our statistical analysis. The 
histograms that we have plotted are also independent of our definition of x^- We have 
checked that if we estimate the allowed ranges from the histograms these are consistent 
with what we get from our definition of x^- Thus we agree with the observation in ref. 



Iql that although this method of data fitting is approximate it works well. 

In [25| the authors have considered another scenario where Am^ < 0.0001 eV^ and 
have obtained good fits to the SK data. They have also discussed in details the neutrino 
decay models and the consequences of such models for astrophysics and cosmology. 
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CHAPTER 7 



Massive Neutrinos in Supernova 



The core of a massive star (M > 8Mq) starts collapsing once it runs out of nuclear 
fuel. The collapse continues to densities beyond the nuclear matter density after which 
a bouncing of the infalling matter takes place leading to supernova explosion and the 
formation of a protoneutron star Only a small fraction of the huge gravitational 
energy released in the process goes into the explosion and all the rest of the energy is 
carried away by neutrinos and antineutrinos of all three flavors. About 10^^ neutrinos, in 
all three flavors carrying a few times 10^^ ergs of energy are released in a type II supernova. 
The luminosities of all the neutrino species are almost same while the average energies 
are approximately 11 MeV for z/g, 16 MeV for Uf. and the average energy of and z/,- and 
their antiparticles is 25 MeV p|. These neutrinos for galactic supernova events can be 
detected by the current water Cerenkov detectors, the Super-Kamiokande (SK) and the 
Sudbury Neutrino Observatory (SNO) 0. In contrast to the solar, the atmospheric as 
well as the accelerator/reactor neutrinos where one has neutrino flux of a single flavor at 
the source, postbounce supernova neutrinos (antineutrinos) start from the source in all 
three flavors but with u^/ut {u^/ut-) having average energies more than that for z/e(t'e) and 
it is an interesting problem to study whether their flux and their signal at the terrestrial 
z/ detectors get appreciably altered in reaching the earth if neutrinos do oscillate. 

The detection of the SN1987A neutrinos by the water Cerenkov detectors at Kamioka 
1^ and 1MB settled many important issues in the subject of type II supernova theory. 
The observation of neutrinos from any future galactic supernova event will answer the 
remaining questions regarding the understanding of the supernova mechanisms. A galactic 
supernova event will also bring in a lot of information on neutrino mass, which of late, 
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has been an issue of much discussion. The effect of neutrino mass can show up in the 
observed neutrino signal in these detectors in two ways, 

• by causing delay in the time of flight measurements 

• by modifying the neutrino spectra through neutrino flavor mixing 

Massive neutrinos travel with speed less than the speed of light and for typical galactic 
supernova distances ~ 10 kpc, even a small mass results in a measurable delay in the 
arrival time of the neutrino. Many different analyses have been performed before to give 
bounds on the neutrino mass by looking at this delay in the arrival time of the massive 
neutrinos (||^, |^ and references therein). Neutrino oscillations on the other hand convert 
the more energetic i^^/i^rii^ii/i^T) into z/e(z/e) thereby hardening the resultant z/e(z/e) energy 



spectra and hence enhancing their signal at the detector y, |T0|, |TT], . In this chapter 
we give quantitative predictions for the number of neutrino events coming from a typical 
type II supernova at a distance of lOkpc in both SNO and SK and show how the number 
of events for each detection process would change in case oscillations do take place. We 
study in detail the effect of neutrino mass and mixing on the total number of events 
recorded in the detector, the distortion of the neutrino spectra due to oscillations, the 
effect of delay due to neutrino mass on the time response of the signal and the effect of 
oscillations on the delay due to mass recorded at the detector. 

The water Cerenkov detectors detect neutrinos through various charged and neutral 
current processes. The differential number of neutrino events at the detector for a given 
reaction process is 

"^'^ j:j^2N^MUE)a{E)e{E) (7.1) 



dEdt ^ A-kD^ 

where i runs over the neutrino species concerned, N^-{t) = Li,-{t)/ {Ei,-{t)), are the number 
of neutrinos produced at the source where L^^{t) is the neutrino luminosity and {Ej,^{t)) 
is the average energy, cr(£') is the reaction cross-section for the neutrino with the target 
particle, D is the distance of the neutrino source from the detector (taken as lOkpc), n 
is the number of detector particles for the reaction considered and fvX^) is the energy 
spectrum for the neutrino species involved, while e{E) is the detector efficiency as a 
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function of the neutrino energy. By integrating out the energy from eq.{7A_) we get the 
time dependence of the various reactions at the detector. To get the total numbers both 
integrations over energy and time has to be done. 

In section 7.1 we use the eq. ( [7.1| ) to estimate the signal that a future galactic super- 
nova event would register in SK and SNO, using the luminosities and average energies 
from a realistic 20Mq supernova model. We consider a scheme of neutrino mass and 
mixing such that one has almost pure vacuum oscillations and study its effect on the 
neutrino spectrum and hence on the signal at the detector. We next consider in section 
7.2 a neutrino mass spectrum where one of the neutrino masses is in the eV range and 
we study the effect of delay in the arrival time on the neutrino signal for both with and 
without the presence of neutrino flavor mixing. We end this chapter by drawing the main 
conclusions in section 7.3. 

7.1 Effect of Neutrino Oscillations in Vacuum 

There have been various attempts before to estimate the effect of non-zero neutrino 
mass and mixing on the expected neutrino signal from a galactic supernova. Matter 
enhanced resonant flavor conversion has been observed to have a large effect on the Uf, 
signal 1^, |T^, |Tl], |I2|, The Ue events of course remain unchanged in this case. With 



vacuum oscillations we can expect an increase in both the z/g and Ue signal. Burrows et 
al. have considered for SNO, the effect of vacuum oscillations as well and have found 
that with two-flavors the effect of vacuum oscillations on the signal is small, using their 
model predictions for the different v luminosities. 

We have considered a three-generation mixing scheme and have calculated the effect of 
neutrino oscillations on the signal from a 20 supernova model developed recently [|1^ 
by Totani et al. based on the hydrodynamic code developed by Wilson and Mayle. For 



the neutrino luminosities and average energies we use the model predictions from |n 
Though in their paper Totani et al. observe that the neutrino spectrum is not a pure 
black body, but we as a first approximation use a Fermi-Dirac spectrum for the neutrinos, 
characterized by the u temperature alone for simplicity. The effect of a chemical potential 
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is to cut the high energy tail of the neutrino spectrum and we also study its effect on the 
the V signal and on the enhancement of the signal when oscillations are introduced. 
For the mass and mixing parameters we consider two scenarios 

• Scenario 1: First we do our calculations for the threefold maximal mixing model 
T5| , |16|, consistent with the solar (Am^2 ~ 10^^° eV^) and the atmospheric neutrino 



data (Am23 ~ Am^g ~ 10^^ eV^). For Am^ ~ 10~'^ eV^ normally we expect matter 
enhanced resonance in the supernova. But for the particular case of maximal mixing 
it has been shown before, both numerically |jT9| and analytically that there are 



no discernable matter effects in the resultant neutrino spectrum on earth. Though the 
arguments in both these previous papers are for solar neutrinos, extension to the case of 
supernova neutrinos is straightforward. Hence we are concerned with vacuum oscillations 



only. Since the oscillation wavelengths (cf. eq. ( p.7| )) corresponding to both the mass 
scales are much smaller compared to the Earth-supernova distance L, the oscillation 
probabilities reduce to [|1^, ^ 

P.... = I (7.2) 

Pu^Ue + PvT-Ue ~ ^ Pfel'e (7-3) 

We call this case 1. 

• Scenario 2: Here we set Amfg ~ 10~^^ eV^ for which A ~ L and hence the 
oscillation effects are observable in the neutrino spectrum. The other mass range in kept 
in the solar vacuum oscillation regime Am^g ^ Arrigg ~ lO^^*' eV^. For this case the 
oscillations due to Am^g and Am^g are averaged out as the neutrinos travel to Earth but 
those due to Am^2 survive. In this scenario again there is no matter effects and one has 
vacuum oscillations. The transition and survival probabilities in this case are 

irT 1 

P,^,^ = 1 - sin^ 2^12 cos^ sin^ — - - sin^ 2^^i3 (7.4) 

Ai2 2 

We use the for the mixing matrix U = R23R13R12, where Rij have been defined before in 
chapter 5. For ^^13 we consider two sets of values allowed by the solar u data. We have 
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done our calculations for sin^ 26'i3 = 1.0 and sin^2^i3 = 0.75. The first set is called Case 
2a while the second is called Case 2b. Since nothing constrains Am^2 this scenario 
we can vary 612 and study its effect on the u signal. We have tabulated our results for 
sin^ 29i2 = 1.0 since it gives the maximum increase in the signal from the no oscillation 
value. 

The corresponding expressions for the antineutrinos will be identical. We note that 
because the energy spectra of the and Uj. are identical, we do not need to distinguish 
them and keep the combination Puf,ue + PurUe- We have made here a three-generation 
analysis where all the three neutrino flavors are active. Hence if both the solar u problem 
and the atmospheric u anomaly require u oscillation solutions, then in the scenario 2, 
the atmospheric data has to be reproduced by z/^ — z/gteriie oscillations^. We are interested 
in this scenario as only with neutrinos from a supernova can one probe very small mass 
square differences ~ 10~^^ eV^. To find the number of events with oscillations we will have 
to fold the expression ( [7. 1| ) with the expressions for survival and transition probabilities 
for the neutrinos for all the cases considered. 

In Table [7.1| we report the calculated number of expected events for the main reactions 



in H2O and D2O. Column 2 of Table |7[T] gives the expected numbers for the supernova 
model under consideration when the neutrino masses are assumed to be zero. Column 
3,4,5 give the corresponding numbers for the two scenarios of neutrino mixing that we have 
considered (see Table ^]l]for details). All the numbers tabulated have been calculated for 
1 kton of detector mass. To get the actual numbers we have to multiply these numbers 
with the relevant fiducial mass of the detector. The efficiency of both the detectors (SNO 
and SK) is taken to be 1 p, 0, [l^. The energy threshold is taken to be 5 MeV for both 
SK H and SNO 0. The energy threshold of SNO in the recently declared solar neutrino 



results is 6.75 MeV [|23|. For the cross-section of the {ue — d), (z/g — d), {u^ — d) and (ue—p) 



reactions we refer to |2^. The cross-section of the (z/e(i/e) — e ) and (i^x — e ) scattering 
has been taken from while the neutral current {u^ — O) scattering cross-section 



is taken from For the ^^0(z/e,e Y^F and ^^0{i'f;,e~^y^N reactions we refer to 

^^The pure — I'stcriic solution to the atmospheric neutrino problem is now disfavored from the SK 
atmospheric data pl| . However schemes in which j/^ can oscillate into a combination of active and sterile 
species are still allowed |2^ . 
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Table 7.1: The expected number of neutrino events for a 1 kton water Cerenkov detector 
(H2O or D2O). 



where we have used the cross-sections for the detector with perfect efficiency. From a 



comparison of the predicted numbers in Table |7]T|, it is evident that neutrino oscillations 
play a significant role in supernova neutrino detection. For the neutral current sector the 
number of events remain unchanged as the interaction is fiavor blind. 

The 32 kton of pure water in SK detects neutrinos primarily through the capture 
of Ue on protons {ueP — ^ ne~^) and (t'e(^'e) — e~) scattering. The energy threshold for 
^^0{ue, e-y^F is 15.4 MeV and that for ^^0{ue, e+Y^N is 11.4 MeV, hence these reactions 
are important only for very high energy neutrinos. The typical average energies of z/g and 
z/g from a type II supernova is about 11 MeV and 16 MeV respectively, so we do not expect 



significant contribution from these two reactions. This is evident from Table |7.1| where 
the ^^O events are only 2.1% of the total charge current signal at SK. As a result of mixing 
the mu and tau neutrinos and antineutrinos oscillate (with average energy ~ 25 MeV) 
into z/g and z/g during their fiight from the galactic supernova to the detector resulting in 
higher energy z/g and z/g and the number of ^^O events are increased appreciably (for Case 
1 (z/g — O) events go up by 13 times) so that after oscillations they are 7% (Case 1) 
of the total charge current events at SK. The effect of oscillations on the (z/g-p) capture 
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is to enhance the expected signal by about 25% (Case 1). In all previous studies where 



the effect of MSW transition on the neutrino signal has been studied [|T3|, , there is no 
enhancement in the number of expected events for the (i^e-p) sector while we do get a 
significant change in the expected signal with vacuum oscillations. For the (z/e(t'e) — e~) 
scattering the effect of oscillation is very small. 

The SNO is the world's first heavy water detector made of 1 kton of pure D2O sur- 
rounded by ultra pure H2O. We find about 99% increase in (i/g — d) events and about 
46% increase in (z/g — d) events for the Case 1. From the column 2 of Table |7.1| we can see 
that there are more (z/g — d) than (z/g — d) events even though there are more u^. than 
coming from the supernova. This is because the reaction cross-section a ~ E"^'^ and the 
Ue spectrum is harder than the spectrum. This also results in a greater enhancement 
due to oscillations for the (z/g ~" d) events, as the difference between the energies of the 
z/e and z/^(z/t-) is greater than those between z/g and z/^(z/t-) and hence the effect on the Uf, 
events is more. As a result after oscillations are switched on the number of (z^e — d) events 
supersede the (z^e ~ d) events. We observe a similar effect for the ^^O events, where the 
z/g signal without oscillations is more than the z/e signal, while the effect of oscillations 
is more for the latter. The effect is more magnified in this case due to the very strong 
energy dependence of the reaction cross-section and also due to the fact that the energy 



threshold for (z/g — O) event is lower than for the — ^® O) event. In fig. [7.1| we plot 
the signal due to the (z/g — d) events as a function of neutrino energy, without oscillations 
and with oscillations for the Case 1 and Case 2b. All the features mentioned are clearly 
seen. The plot for the Case 2b clearly shows oscillations. 

In fig. |7]^ we plot the cumulative fluence of the coming from the supernova at 10 
kpc without oscillations and with oscillations for Case 1 and Case 2b. It is seen that the 



result of oscillation in fact is to reduce the total number of z/g. Yet as seen from Table 7.1 



we have obtained significant increase in the (z/g — d) events and the (z/g — O) events. The 
solution to this apparent anomaly lies in the fact that the cross-section of these reactions 
are strongly energy dependent. As a result of oscillations the z/g flux though depleted in 
number, gets enriched in high energy neutrinos. It is these higher energy neutrinos which 
enhance the z/ signal at the detector. This also explains the difference in the degree of 
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Figure 7.1: The {u^ — d) signal at SNO vs neutrino energy without (sohd hne) and with 
oscillations for the Case 1 (short dashed line) and Case 2b (long dashed line). 

enhancement for the different processes. For the (i/g — d) and (i/g — O) events, especially 
for the latter, the effect is huge while for the {ue — e") scattering it is negligible as its 
reaction cross-section is only linearly proportional to the neutrino energy. Due to their 
high energy dependent cross-sections the ^^0{i'e,e^y^F events turn out to be extremely 
sensitive to oscillations. A similar argument holds true for the case of the antineutrinos, 
only here the effect of oscillations is less than in the case for the neutrinos as the difference 
between the energies of the Pg a-nd P^/ P^ is comparatively less as discussed earlier. 

For the scenario 2 we have studied the effect of the mixing angles on the signal. For 
a fixed ^13 the effect of oscillations is enhanced if we raise 612- The effect of 613 is more 
subtle. The effect of oscillations increase with ^13 initially and then decrease. We have 
also checked the effect of a chemical potential /i on the neutrino signal. A non-zero /i 
cuts the high energy tail of the neutrino signal as a result of which the total signal goes 
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Figure 7.2: The cumulative z/g fluence as a function of the neutrino energy without (sohd 
hne) and with oscillations for the Case 1 (dotted line) and Case 2b (short dashed line). 
Also shown is the fluence (long dashed line) for comparison. 



down for both with and without oscillations, the effect being greater for the more energy 
sensitive reactions. 

With the supernova model of Totani et al. [H], we have obtained oscillation effects 
in the expected v signal which are significantly larger than those obtained by Burrows 
et al. 0, ^. In the model that Burrows et al. use in their study, the v luminosities 
Ly are more than those for Totani et al. model, but the average energy is much smaller, 
particularly for the z/g and ty^^rii^fj^T)- Hence their z/^ spectra lacks in high energy neutrinos 
which results in almost negligible effect of oscillations in their case. Again in the model 
of Burrows et al. the average energies decrease with time while in the model of Totani et 
al. not only the average energies but also the difference between the average energies of 
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Ve{^e) and v^^ri^^^r) incieases with time. The effect of all these is to magnify the effect 
of oscillations in our case. 

7.2 The Effect of Delay in the Time of Flight 

In the previous section we studied the effects of neutrino flavor oscillations on the 
supernova neutrino spectrum and the number of charged current events at the detector 
using a realistic supernova model. In this section we study the neutral current signal as 
a function of time in the water Cerenkov detectors, for a mass range of the neutrinos 
where both the phenomenon of delay and flavor conversion are operative. That the time 
response of the event rate in the detector is modified if the neutrinos have mass alone and 
hence delay is a well known feature p, 0]. Here we stress the point that since neutrino 
flavor conversions change the energy spectra of the neutrinos, and since the time delay of 
the massive neutrinos is energy dependent, the time dependence of the event rate at the 
detector is altered appreciably in the presence of mixing. We suggest various variables 
which act as tools for measuring this change in the time response curve of the neutral 
current events and in differentiating the cases of (a) massless neutrinos (b) neutrinos with 
mass but no mixing and (c) neutrinos with mass as well as mixing. In particular we study 
the ratio of the charged current to neutral current ratio R(t), as a function of time in 
the SNO detector and show that the change in the value and the shape of R(t) due to 
flavor mixing cannot be emulated by uncertainties. We also study other variables like the 
normalized n-th energy moments of the neutral current events and the ratio of charged to 
the neutral current n-th moments as important diagnostic tools in filtering out the effects 
of neutrino mass and mixing. 

For the neutrino luminosities and average energies, though it is best to use a numerical 
supernova model, but for simplicity, we will here use a profile of the neutrino luminosities 
and temperatures which have general agreement with most supernova models. We take 
the total supernova energy radiated in neutrinos to be 3 xlO^^ ergs. This luminosity, 
which is almost the same for all the neutrino species, has a fast rise over a period of 
0.1 sec followed by a slow fall over several seconds in most supernova models. We use a 
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luminosity that has a rise in 0.1 sec using one side of the Gaussian with a = 0.03 and 
then an exponential decay with time constant r = 3 sec for all the flavors |^. 

The average energies associated with the z/g, Ue and z/^ (the z/^, z/^, z/^ and z/^ have the 
same energy spectra) are 11 MeV, 16 MeV and 25 MeV respectively in most numerical 
models. We take these average energies and consider them to be constant in time. We 
have also checked our calculations with time dependent average energies and estimated its 
effect. The neutrino spectrum is taken to be a pure Fermi-Dirac distribution characterized 
by the neutrino temperature alone. 

If the neutrinos are massless then the time response of their signal at the detector 
reflect just the time dependence of their luminosity function at the source, which is the 
same for all the three flavors and hence the same for the charged current and neutral 
current reactions. If neutrinos have mass ~ eV then they pick up a measurable delay 
during their course of flight from the supernova to the earth. For a neutrino of mass m 
(in eV) and energy E (in MeV), the delay (in sec) in traveling a distance D (in 10 kpc) is 

At{E) = 0.515{m/EfD (7.6) 

where we have neglected all the small higher order terms. The time response curve then 
has contributions from both the luminosity and the mass. We will now consider a scheme 
of neutrino masses such that Amfj ~ 10^^ eV^ consistent with the LMA MSW solution 
of the solar neutrino problem and Am^g ^ Amgg ~ 1 — 10^ eV^. For Am^ ~ 10^^ 
eV^ the effect of Earth matter can show up in the neutrino signal but we neglect 
it here for simplicity. The neutrino mass model considered here is one of several, given 
for the purpose of illustration only. In this scheme the atmospheric neutrino anomaly 
will have to be explained by the z/^ — z/gteriie oscillation mode. The mass range for the 
neutrinos as the hot component of hot plus cold dark matter scenario in cosmology is a 
few eV only which will conflict with the higher values in the range of m^^ = 1 — 100 



eV that we consider here if z/3 is stable. Hence, we assume that the z/3 state is unstable 
but with a large enough life time so that it is does not conflict with the observations of 
SN 1987A ||29[ (even though SN1987A observations did not correspond to any Ur event. 



one can put limits on the Z/3/Z/3 lifetime as the z/g/z/g state is a mixture of all the three 
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mass eigenstates) and is also consistent with Big Bang Nucleosynthesis. In fact, from 
the ref. we know that using the time delay technique, the SK and SNO can be 
used to probe neutrino masses down to only 50 eV and 30 eV respectively. Hence we 
have presented all our results for a particular representative value of rriu-i — 40 eV. There 
have been proposals in the past for an unstable neutrino with mass ~ 30 eV and lifetime 
~ 10^'^ s |^D[. Since direct kinematical measurements give m^^ < 5 eV ||^, we have taken 



the z/g to be massless and the charged current events experience no change. But since 
the ^'T(z^r) pick up a detectable time delay (for the mass spectrum of the neutrinos that 
we consider here, the v^iv^) do not have measurable time delay), the expression for the 
neutral current events gets modified to, 
(jQd n r 

^ = ^ y dEa{E){N^St)USE) + N.MUAE) + K^U^E) + 

+ N,^{t)USE) + KAt-^m)hAE) + N,^{t-At{E))U^{E)} (7.7) 



where dS^^/dt denotes the neutral current (nc) event rate with delay (d). Delay therefore 
distorts the neutral current event rate vs. time curve. By doing a analysis of this 
shape distortion one can put limits on the Ur mass [§, 0]. 

We next consider the neutrinos to have flavor mixing as well. The mixing angle sin^ 9i2 
can be constrained from the solar neutrino data (cf. chapter 3) while for sin^ 6*13 there 
is no experimental data to fall back upon, but from r-process considerations in the 'hot 
bubble" of the supernova, one can restrict sin^ ^13 ~ 10^^ [10[. In this scenario there 
will be first a matter enhanced u^. — resonance in the mantle of the supernova followed 
by a z/g — resonance in the envelope. The MSW mechanism in the supernova for the 



neutrino mass scheme that we consider here is discussed in details in ref. |Tn|. As the 
average energy of the v^jv^ is greater than the average energy of the i^e? neutrino flavor 
mixing modifies their energy spectrum. Hence as pointed out in the previous section the 
Vf. flux though depleted in number, gets enriched in high energy neutrinos and since the 
detection cross-sections are strongly energy dependent, this results in the enhancement 
of the charged current signal. The total number of events in SNO, integrated over time 
in this scenario are given in Table IJl. In the third column of Table |7.2| (marked A) we 



show the number of events for galactic supernova neutrinos with luminosities and average 
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Table 7.2: The expected number of neutrino events in SNO. To get the number of events 
in SK, one has to scale the number of events in H2O given here to its fiducial mass of 
32 kton. The column A corresponds to massless neutrinos, column B to neutrinos with 
complete flavor conversion The Vi here refers to all the six neutrino species. 



energies considered in this section while the last column (marked B) gives the number of 
events with complete flavor oscillations [Pv^v^ = 0)- Of course since the z/g do not have any 
conversion here, the z/g signal remains unaltered. Also as the neutral current reactions are 
flavor blind, the total neutral current signal remains unchanged. But whether the time 
response curve of the neutral current signal remains unchanged in presence of mixing, in 
addition to delay, is an interesting question. 

If the neutrinos have mass as well as mixing, then the neutrinos are produced in their 
flavor eigenstate, but they travel in their mass eigenstate. The neutrino mass eigenstates 
will travel with different speeds depending on their mass and will arrive at the detector 
at different times. For the scenario that we are considering only z/3 and z/3 will be delayed. 
Hence to take this delay in arrival time into account, the eq.( |7.7| ) has to be rewritten in 
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terms of the mass eigenstates. It can be shown that expression for the neutral current 
event rate in terms of the mass eigenstates is, 



dEa{E){K,{t)U,{E) + N,,{t)U,{E) + K,{t)U,{E) 



dt AttD^ 

+ N,,{t)U,{E) + K.^{t - /\t{E))U,{E) + N,,{t - /\t{E))U,{E)} (7.8) 

where N^- is the Vi flux at the source. If the neutrinos are produced at densities much 
higher than their resonance densities, all the mixings in matter are highly suppressed, 
and the neutrinos are produced almost entirely in their mass eigenstates. For the three 
generation case that we are considering, ~ J^s, ^'/^ ~ and For the antineutri- 

nos on the other hand, at the point of production in the supernova z/g ~ ^i, ^'^t ~ ^2 and 
z/t- ^ z/3. Hence the above expression for the neutral current event rate in the presence of 
delay and mixing can be written as. 



do 

' dEa{E){N,^{t)U^{E) + N,Xt)UAE) + iV,^(t)/..(i?) + 



dSt 



dt AttD'^ 

+ N,^{t)U^{E)+N,Xt-^m)fuAE) + N,Xt-^t{E))U^{E)] (7.9) 

Note that the above expression does not depend on the neutrino conversion probability 
as the neutral current interaction is flavor blind. 



In fig. [7.3| we have plotted the neutral current event rate for the reaction (z/j + d ^ 
n + p + Vi, where Vi stands for all the 6 neutrino species) as a function of time for massless 
neutrinos along with the cases for mass but no mixing (eq.( |7.7| )) and mass along with 
mixing (eq.( |7.9| )). The figure looks similar for the other neutral current reactions as well, 
apart from a constant normalization factor depending on the total number of events for 
the process concerned. The curves corresponding to the massive neutrinos have been 
given for mj^g = 40 eV. As expected, the shape of the neutral current event rate changes 
due to the delay of massive Vr- Since the delay given by eq.( [7!6D depends quadratically 
on the neutrino mass, the distortion is more for larger masses |^2|. But the noteworthy 
point is that the presence of mixing further distorts the rate vs. time curve. The reason 
for this distortion can be traced to the fact that the time delay oc l/E"^. As the energy 
spectrum of the neutrinos change due to flavor mixing, the resultant delay is also modified 
and this in turn alters the neutral current event rate as a function of time. In fact the 
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Figure 7.3: The neutral current event rate as a function of time in D2O in SNO. The sohd 
hne corresponds to the case of massless neutrinos, the long dashed line to neutrinos with 
only mass but no mixing, while the short dashed line gives the event rate for neutrinos 
with mass as well as flavor mixing. 



flavor conversion in the supernova results in de-energising the u^/ Vj. spectrum and hence 
the delay given by eq.( [7.6|) should increase. As larger delay caused by larger mass results 
in further lowering of the neutral current event rate vs. time curve for early times, one 
would normally expect that the enhanced delay as a result of neutrino flavor conversion 
would have a similar effect. But the fig. [7.3| shows that during the first second, the curve 
corresponding to delay with mixing is higher than the one with only time delay. This 
at first sight seems unexpected. But then one realizes that while the flavor conversion 
reduces the average energy of the massive Vr increasing its delay and hence depleting its 
signal at early times, it energizes the massless and hence undelayed z/g beam, which is 
detected with full strength. Therefore, while for no mixing the v^- gave the larger fraction 
of the signal, for the case with mixing it is the that assume the more dominant role, 
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Figure 7.4: The ratio R(t) of the total charged current to neutral current event rate 
in SNO versus time. The solid line is for massless neutrinos, the short dashed line for 
neutrinos with complete flavor conversion but no delay, the long dashed line for neutrinos 
with only delay and no flavor conversion and the dotted line is for neutrinos with both 
delay and complete flavor conversion. Also shown are the ±1(7 statistical errors for delay 
with and without mixing in the 1**, 4*'* and the 7*'* time bins. 



and so even though the Uj- arrive more delayed compared to the case without mixing, the 
delay effect is diluted due to the enhancement of the Ug fraction and the depletion of the 
Ut- fraction of the neutral current events. We have also checked that although it may seem 
that the curve with delay and mixing can be simulated by another curve with delay alone 
but with smaller mass, the actual shape of the two curves would still be different. This 
difference in shape though may not be statistically signiflcant and hence one may not be 
able to see the effect of mixing in the time delay of the neutrinos just by looking at the 
time response of the neutral current event rate in the present water Cerenkov detectors. 
We therefore look for various other variables which can be studied to compliment this. 
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One such variable which carries information about both the neutrino mass and their 
mixing is R(t), the ratio of charged to neutral current event rate as a function of time. 
In fig. we give the ratio R(t) of the total charged current to the neutral current event 
rate in D2O in SNO as a function of time. Plotted are the ratios (i) without mass, (ii) 
with only mixing, (iii) with delay but zero mixing and (iv) with delay and flavor mixing. 
The differences in the behavior of R(t) for the four different cases are clearly visible. For 
no mass R(t)=0.3 and since the time dependence of both the charged current and neutral 
current reaction rates are the same, their ratio is constant in time. As the presence of 
mixing enhances the charged current signal keeping the neutral current events unaltered, 
R(t) goes up to 0.61 for only mixing, remaining constant in time, again due to the same 
reason. With the introduction of delay the ratio becomes a function of time as the neutral 
current reaction now has an extra time dependence coming from the mass. At early times 
as the Pr get delayed the neutral current event rate drops increasing R(t). These delayed 
z/t-s arrive later and hence R(t) falls at large times. This feature can be seen for both the 
curves with and without mixing. The curve for only delay starts at R(t)=0.52 at t=0 s 
and falls to about R(t)=0.26 at t=10 s. For the delay with mixing case the corresponding 
values of R(t) are 0.83 and 0.51 at t=0 and 10 s respectively. The important point is 
that the curves with and without mixing are clearly distinguishable and should allow 
one to differentiate between the two cases of only delay and delay with neutrino flavor 
conversion. 

In order to substantiate our claim that the two scenarios of only delay and delay with 
mixing are distinguishable in SNO, we divide the time into bins of size 1 second. The 
number of events in each bin is then used to estimate the ±lcr statistical error in the ratio 



R(t) in each bin and these are then plotted in flg. 1_A for the typical time bin numbers 
1, 4 and 7. From the flgure we see that the two cases of delay, with and without mixing, 
are certainly statistically distinguishable in SNO for the flrst 6 seconds. 

We next focus our attention on M"^(t), the neutral current ra-th moments of the 



neutrino energy distributions |Q observed at the detector, deflned as 
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Figure 7.5: The 1st normalized energy moment of the neutral current events in SNO 
m7'^(^) versus time. The solid line corresponds to the case of massless neutrinos, the long 
dashed line to neutrinos with only mass but no mixing, while the short dashed line gives 
the event rate for neutrinos with mass as well as flavor conversion. Also shown are the 
±1(7 statistical errors for delay with and without mixing in the l'^*, 6^^ and the ll*'^ time 
bins. 



while the corresponding normalized moments are given by 



Ml (t) = " 7.11 



We have shown the behavior of the 1st normalized moment (t) in fig. |7]^ as a function 
of time in SNO. For massless neutrinos, the M^^ has a value 40.97, constant in time, as 
this is again a ratio and hence the time dependence gets canceled out as in the case 
of R(t). For the case where the u^- is massive and hence delayed, it assumes a time 
dependence. Since the delay oc l/E"^ and since the neutrinos are produced at the source 
with an energy distribution, hence at each instant the lower energy Ur will be delayed 
more than the higher energy Ur- Therefore M^'^(t), which gives the energy centroid of the 
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neutral current event distribution in D2O, starts from a low value 38.76 at t=0 s as all 
the Vt- are delayed, rises sharply as the higher energy neutrinos arrive first and then falls 
slowly as the lower energy delayed Vr start arriving. If the Vr are allowed to mix with the 
Vei then they are de-energized and the above mentioned effect is further enhanced. To 
make an estimate of whether SNO would be able to distinguish the three cases discussed 
above, we compute the ±lcr statistical errors in the f** normalized moment for the two 
scenarios of delay, with and without mixing, and show them for the 1**, and 11*'' bins. 
We see that the errors involved are large enough to completely wash out the differences 
between the energy moments with and without neutrino mass and mixing. Hence the 
normalized energy moments fail to probe neutrino mass and mixing as at early times we 
don't see much difference between the different cases considered, while at late times the 
number of events become very small so that the error in M^'^(t) becomes huge, increasing 
the error in M^'^(t). 

The variable that can be a useful probe for differentiating the case for delay with 
mixing from the case for delay without mixing is the ratio of the unnormalized moment 
of the charged to neutral current events 

r it) - ^ (7 12) 

We present in fig. \l.6[ for SNO, the r„(t) vs. time plots (for n=l) for the cases of (a) 
massless neutrinos (b) with mixing but no delay (c) with delay but no mixing and (d) 
with delay as well as mixing. Since this is a ratio, the supernova flux uncertainties get 
canceled out to a large extent and since the unnormalized moments have smaller statistical 
errors, this is a better variable than the normalized moments to observe the signatures of 
neutrino mixing. In the flgure we have shown the ±la statistical errors in ri{t) for the 
two cases of delay alone and delay with mixing, for the f**, 8*^ and 15*^ bins in time, and 
the two cases are clearly distinguishable in SNO for early as well as late times. Note that 
ri{t) is different from the ratio R(t) as it gives information about the ratio of the energy 
centroids of the charged current and neutral current distributions as a function of time, 
while the latter gives only the ratio of the number of events as a function of time. 
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Figure 7.6: The variation of ri{t) with time in SNO. The sohd hne is for masslcss 
neutrinos, the short dashed hne for neutrinos with complete flavor conversion but no 
delay, the long dashed line for neutrinos with only delay and no flavor conversion and the 
dotted line is for neutrinos with both delay and complete flavor conversion. Also shown 
are the ±1(7 statistical errors for delay with and without mixing in the 1**, S*'* and the 
15*'* time bins. 

The advantage of using ratios is that, they are not only sensitive to the mass and 
mixing parameters but are also almost insensitive to the details of supernova models. 
Since they are a ratio they are almost independent of the luminosity and depend only on 
some function of the ratio of neutrino temperatures. All the calculations presented so far 
are for flxed neutrino temperatures. In order to show that the time dependence of the 
neutrino temperatures does not alter our conclusions much, we present our analysis with 
time dependent neutrino temperatures. We take 

= 0.16 log i + 3.58, Tp^ = 1.63 log i + 5.15, 7;^ = 2.24 log i + 6.93 (7.13) 
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Figure 7.7: The ratio R(t) in SNO for the two cases of fixed and time dependent 
neutrino temperatures. The sohd hne and the long dashed hne give R(t) for the cases of 
fixed temperatures and varying temperatures respectively for only delay, while the short 
dashed line and the dotted line give the corresponding R(t) for delay with mixing. We 
have also given the ±lo" statistical errors in the 1*** and the 4*'* time bin, for the both the 
curves for fixed and time dependent temperatures. 



These forms for the neutrino temperatures follow from fits to the results of the numerical 



supernova model given in Totani et al. ||T4| which we had used in the previous section. 



In fig. 7/7 we compare the ratio R(t) for the cases of delay and delay with mixing for 
the two cases of fixed temperatures and the time dependent temperatures. It is clear 
from the figure that that the time dependence of the neutrino temperatures does not have 
much effect on the time dependence of the ratio of the charged current to neutral current 
rates. In fact the two curves corresponding to fixed and time dependent temperatures, fall 
within ±1(7 statistical errorbars for both the cases of only delay and delay with mixing. 
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7.3 Summary 



In conclusion, we have shown that with the model of Totani et al. even with vacuum 
oscillations we obtain appreciable enhancement in the expected 1/ signal in SNO and SK 
even though the number of neutrinos arriving at the detector from the supernova goes 
down. In contrast to the case where we have MSW resonance in the supernova, with 
vacuum oscillations we get enhancement for both 1/^ as well as Pg events. If we have a 
galactic supernova event in the near future and if we get a distortion in the neutrino 
spectrum and an enhancement in the signal, for both i/g as well as Pg then that would 
indicate vacuum neutrino oscillations. 

We have shown that even though neutrino flavor mixing cannot alter the total neutral 
current signal in the detector - the neutral current interaction being flavor blind, it can 
have a non-trivial impact on the delay of massive neutrinos, which alters the neutral 
current event rate as a function of time. The neutral current event rate though does not 
depend on the neutrino conversion probability. In order to study the effect of neutrino 
mass and mixing we have suggested various variables. Of the different variables that we 
have presented here, the ratio of the charged to neutral current event rate R(t), can show 
the effect of mixing during the first few seconds, while the charged to neutral current ratio 
of the energy moments arc useful diagnostic tools for all times. These variables arc not 
just sensitive to flavor mixing and time delay, they are also insensitive to supernova model 
uncertainties and hence are excellent tools to study the effect of flavor mixing on the time 
delay of massive supernova neutrinos. Though we have considered a mass spectrum for 
the neutrinos where only the 1/3 have a measurable delay but the model considered is one 
of many and one can easily extend the above formalism to include more general classes 
of neutrino models. 
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CHAPTER 8 



Conclusions 



8.1 Summary 

In this thesis we have explored the signatures of neutrino mass and mixing in the 
solar and atmospheric neutrino experiments. We have derived bounds on the mixing 
parameters through elaborate analyses. We have studied the effect of the neutrino 
mass and mixing on supernova neutrino detection. 

In chapter 2 we have presented the basic aspects of neutrino oscillations, both in 
vacuum and in matter. 

In chapter 3 we have discussed the experimental status of the solar neutrino problem, 
the atmospheric neutrino anomaly and the accelerator /reactor neutrino experiments. We 
performed comprehensive statistical analysis of the solar and atmospheric neutrino data in 
terms of two flavor oscillations, identified the best-fit solutions and presented the allowed 
area in the neutrino parameter space. For the solar neutrino problem we have presented 
our results in Table p.3| . The LMA MSW solution emerges as the best-fit solution with 
Xmin — 33.42 which for 39 degrees of freedom is allowed at 72.18%. The other solutions 
with large mixing angles also give good fit. The CL. contours for the global analysis 
of the rates and rates-|-spectrum data are displayed in fig. K7\ where we note that after 



the inclusion on the SNO CC results the SMA solution fails to appear even at 3a. The 
t'e — t'steriio optiou is largely disfavored. For the atmospheric neutrino case we defined two 
methods of analysis and discussed their merits and demerits. We performed the bin- 
by-bin analysis of the 1144 day SK contained events in terms of pure — 1/^ oscillations 
and presented the results. The 90% and 99% CL. allowed zones in the parameter space 
are displayed in fig. p.lO . We gave an outline of the experimental bounds on the mixing 



parameters from the most stringent accelerator /reactor data. 
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In chapter 4 we continued our discussions on the solar neutrino problem and explored 
the viability of the energy independent scenario in explaining the global solar data. We 
identified regions of the parameter space where the survival probability is within 10% of 
the energy independent survival probability and called these the quasi-energy independent 
regions. Allowing for modest renormalizations for the CI data and the BPBOO flux 
we studied the comparative flt for the energy independent solution vis a vis the MSW 
solutions. We showed that these renormalizations enlarge the allowed large mixing angle 
MSW regions and most of these enlarged regions overlap with the quasi-energy indepen- 
dent region. We briefly commented on the potential of some of the future experiments in 
distinguishing the energy independent scenario from the MSW solutions. 

In chapter 5 we did a three-generation oscillation analysis of the 1144 day (SK) atmo- 
spheric neutrino data going beyond the one mass scale dominance (OMSD) approximation. 
We fixed A12 = A13 {Alsnd) in the range eV^ as allowed by the results from LSND and 
other accelerator and reactor experiments on neutrino oscillation and kept A23 (A^^m) 
and the three mixing angles as free parameters. We incorporated the matter effects, indi- 
cated some new allowed regions with small A23 (< 10~^ eV^) and sin^ 2623 close to and 
discussed the differences with the two-generation and OMSD pictures. In our scenario, 
the oscillation probabilities for the accelerator and reactor neutrinos involve only two of 
the mixing angles $12 and ^13 and one mass scale. But the atmospheric neutrino oscillation 
is in general governed by both mass scales and all the three mixing angles. The higher 
mass scale gives rise to Am^ independent average oscillations for atmospheric neutrinos 
and does not enter the analysis as an independent parameter. The A23 and the three 
mixing angles on the other hand appear as independent parameters in the analysis and 
the best-fit values of these are determined from an analysis of a) the SK data, b) the SK 
and CHOOZ data. The allowed values of the mixing angles 612 and ^13 from the above 
analysis are compared with the constraints from all accelerator and reactor experiments 
including the latest results from LSND and KARMEN2. Implications for future long 
baseline experiments are discussed. 

In chapter 6 we did a detailed x^-analysis of the 848 day SK atmospheric neutrino data 
under the assumptions of — Uj. oscillation and neutrino decay. For the latter we took 
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the most general case of neutrinos with non-zero mixing and considered the possibihties 
of the unstable component in i/^ decaying to a state with which it mixes (scenario (a)) 
and to a sterile state with which it docs not mix (scenario (b)). In the first case Am^ 
(mass squared difference between the two mass states which mix) has to be > 0.1 eV^ 
from constraints on K decays while for the second case Am^ can be unconstrained. For 
case (a) Am^ does not enter the x^-analysis while in case (b) it enters the x^-analysis as 
an independent parameter. In scenario (a) there is Am^ averaged oscillation in addition 
to decay and this gets ruled out at 100% by the SK data. Scenario (b) on the other hand 
gives a reasonably good fit to the data for Am^ ~ 0.001 eV'^. We discussed the possibihty 
of differentiating this latter scenario from the pure i/^ — Vj- scenario with more statistics 
on Vt appearance in SK. 

Neutrinos and antineutrinos of all three flavors are emitted during the post bounce 
phase of a core collapse supernova with t'/i/^'r(^'^/^'r) having average energies more than 
that of Ve{Ve). These neutrinos can be detected by the earth bound detector hke the 
SK and SNO. In chapter 7 we made realistic predictions for the observed signal in the 
detector. We studied the effect of three flavor oscillations in vacuum on the resultant 
neutrino spectra, in the framework of two different mass spectrum and presented our 
results. We showed that even though neutrino oscillations result in a depletion in the 
number of ^-nd coming from the supernova, the actual signals at the detectors are 
appreciable enhanced. In particular we found a huge enhancement in the ^^O charged 
current rate in the water Ccrcnkov detectors due to oscillations. Wc next considered a 
mass spectrum where one has detectable delay in the time of flight of the massive Vj. and 
studied its effect. We showed that even though the neutral current interaction is flavor 
blind, and hence neutrino flavor mixing cannot alter the total neutral current signal in the 
detector, it can have a non-trivial impact on the delay of massive neutrinos and alters the 
neutral current event rate as a function of time. We have suggested various variables of 
the neutral and charged current events that can be used to study this effect. In particular 
the ratio of charged to neutral current events can be used at early times while the ratio 
of the energy moments for the charged to the neutral current events can form useful 
diagnostic tools even at late times to study neutrino mass and mixing. 
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8.2 Future Prospects 



A pertinent question at this point is the potential of the future experiments in building 
upon the current information about the mass and mixing scenario. They should also look 
towards confirming the presence of neutrino oscillations in the deficit of the solar and the 
atmospheric neutrinos. Though neutrino oscillations is the favored solution for both the 
problems, there is no conclusive evidence in its favor. For the atmospheric problem one 
should be able to see the periodicity, or in other words the dips of the oscillations in the 
resultant neutrino beam. Likewise for the solar neutrinos the future experiments should 
be designed to give smoking gun evidence for the oscillation scenario and should have the 
strength to distinguish between the various allowed solutions. 

The MONOLITH experiment in Gran Sasso, Italy is a magnetized tracking Calorime- 
ter which will detect atmospheric neutrinos. This has the sensitivity in L/E to detect 
unambiguously the dip in survival probability predicted by the z/^ — Pr oscillation hypoth- 
esis. The other planned experiment which has the potential to observe the oscillations in 
the atmospheric neutrino signal is the ring imaging Cerenkov detector AQUA RICH also 
in Italy. 

Among the solar neutrino experiments the most promising upcoming detector is the 
BOREXINO in Gran Sasso, Italy. The BOREXINO will detect the monoenergetic Be 
neutrinos via — e scattering. This detector will be able to put strong bounds on the 
already disfavored SMA solution which predicts a very small survival probability for the 
"^Be neutrinos. BOREXINO also expects to see large day-night asymmetry in the LOW 
region and can put constraints there. LENS also in Gran Sasso will be able to detect 
pp, ''Be , pep and ^B neutrinos and hopes to compliment the results of BOREXINO. The 
other low energy detectors are HERON, HELLAZ and the Lithium- Berrylium detector. 

The Kamland in Kamioka, Japan will detect both solar neutrinos as well as from 
ten different reactors all across Japan with various baselines within 350 km. Using these 
known fluxes from the different reactors with known baselines, it is expected to scan the 
entire allowed LMA zone, which is currently the favored solution to the solar neutrino 
problem. The other exciting experiment is MINOS, which with a baseline of 730 km from 
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Fermilab to Soudan will use the disappearance channel to explore almost the entire 
area allowed by the atmospheric neutrino data. The other experiment probing this region 
of the neutrino parameter space is K2K in Japan. It has already given positive signal for 
oscillations. 

The MiniBooNE experiment at Fermilab will probe the oscillations of muon neutrinos 
to electron neutrinos and it is expected that within one year of data taking it will have 
the sensitivity to either confirm or refute the LSND claim of having seen oscillations. If 
MiniBooNE does see oscillations then there are plans to upgrade it to BooNE which will 
have better sensitivity and hence will determine the mass and mixing parameters more 
accurately. 

In conclusion, the large mixing angle solutions, the LMA, LOW and the vacuum os- 
cillation solutions are the favored solution to the solar neutrino problem with the LMA 
giving the best- fit. The — v^. oscillations give an excellent description of the SK at- 
mospheric neutrino data, which has confirmed the existence of neutrino mass. However 
one still needs smoking gun evidence to establish neutrino oscillations, particularly for 
the solar neutrinos. The future experiments should be able to throw more light on this 
question and are also expected to put severe constraints on the allowed values of the mass 
and mixing parameters. 
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